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PREFACE 


In our previous book, Mathematical Biophysics: Physico- 
mathematical Foundations of Biology, we attempted to present 
a system of matkematical biology based on physical concepts. 
While we tried, wherever it was possible, to establish contact 
between mathematical conclusions and experimental facts, 
still the theoretical foundations occupied a much more promi- 
nent place than their actual applications. 

During the two years which have elapsed since the comple- 
tion of the manuscript of that book, considerable progress has 
been made in this field, especially in the direction of applica- 
tions of the mathematical theory to various observations. 
This progress has largely been due to the use of the new ap- 
proximation method, very briefly outlined in the Appendix of 
the former book. In the present book it is used as the founda- 
tion for all the developments of chapters i-vi. Many prob- 
lems, which offered very little hope for a solution in the near 
future by means of the standard methods used before, have 
found their solutions by means of the new approximation 
method. The latter also considerably simplifies the mathe- 
matics used and wiU, we hope, make the present book acces- 
sible to a much wider circle of readers. 

Even more than the previous book, the present one owes a 
very large amount of most important material to the work of 
others. Without the important researches of A. S. Household- 
er, H. D. Landahl, Robert R. Williamson, Alvin M. Weinberg, 
and Gale Young this book would never have been possible. 
Throughout the volume the reader will find numerous refer- 
ences to the work of the above-mentioned members of the 
Chicago group of mathematical biophysicists. 

Yet, in spite of the progress made, the author is perfectly 
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well aware that he is not presenting anything like a finished 
product. This book is to be regarded merely as an early mile- 
stone on the road to further progress. 

The year 1939 has witnessed another important progress in 
the field of mathematical biology, namely, the founding of a 
journal devoted exclusively to publications in that field. The 
Bulletin of Mathematical Biophysics is published quarterly by 
the Psychometric Corporation imder the editorship of the 
author of this volmne and is entering now upon its second year 
of existence. The editorial and publication offices are located 
at the University of Chicago, 5822 Drexel Avenue, Chicago, 
Illinois. Practically all the original publications in mathe- 
matical biophysics appear in the Bulletin, which thus con- 
stitutes a record of research in that field. 

The author’s wife, Mrs. Emily Rashevsky, to whom this 
book is dedicated, has been responsible for the preparation of 
the whole manuscript. For the painstaking and tedious work 
of checking all formulas, both in the manuscript and in proofs, 
the author is indebted to Dr. Alston S. Householder, Mr. 
Herbert D. Landahl, Dr. Alvin M. Weinberg, and Mr. Gale 
Young. Thanks are due Mr. Robert R. Williamson for the 
preparation of the Index, and Mr. Robert Dubin for the prep- 
aration of a nmnber of original drawmgs. Furthermore, the 
author is under obligation to Professor N. J. Berrill for per- 
mission to reproduce Figures 2 and 3 from Growth, i, 213, 
1937; to Professors A. V. Hill and D. Y. Solandt for permission 
to reproduce Figures 29 and 30 from Proceedings of the Royal 
Society of London, B, 120, 400, 1936, and 121, 113, 1936; to 
Professor C. J. van der KJaauw for permission to reproduce 
Figures 12, 13, 14, and 26 from Acta biotheoretica. A, 4, Part II, 
1938; and to the Psychometric Corporation for permission to 
reproduce Figures 39, 41-46, 48, 49, and 53-58 from different 
papers published in Psychometrika. He also wishes to express 
his thanks to Professor Robert Chambers for permission to 
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use for measurements a part of his micromotion film of divid- 
ing blastomeres. 

Finally, acknowledgment must be made to the University 
of Chicago Press and particularly to Miss Mary D. Alexander, 
of the Editorial Department, for unfailing co-operation during 
the process of manufacturing the book. This co-operation 
must have been rendered especially difficult by the author’s 
tendency to make minor improvements and last-minute alter- 
ations at the time of proofreading. 

N. Rashevsky 


Chicago, Iliinois 
January 2, 1940 
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CHAPTER I 


DIFFUSION PHENOMENA IN CELLS 

In biology, as in other natural sciences, we find that, in 
spite of the tremendous variety of phenomena observed, cer- 
tain uniformities become apparent upon closer examination. 
The existence of such uniformities, of such a unity in diversity, 
is one of the prerequisites for the possibility of existence of 
any science. Among other uniformities found in biology, one 
of the most conspicuous is the mode of propagation of life on 
earth. While the methods of multiplication of organisms are 
almost infinite in their variety, they aU reduce, in the ultimate 
analysis, to cell division. Again, the details of cell division 
vary greatly from case to case. But everywhere tht fundamen- 
tal aspects of the phenomenon is the same: a cell, at a certain 
stage of its life, spontaneously breaks up, divides into two or 
more cells, which grow further as individual living units. It 
is in keeping with the spirit of all exact sciences to look for a 
simple common cause of such a general phenomenon. To effect 
such a division, some force of unkn own origin must exist. The 
very commonness of the cell division suggests that there must 
be some common cause underlying all cases of cell division 
and responsible for the force that produces it, though addi- 
tional factors may strongly modify the manifestations of cell 
division from case to case. However, if we assume such a 
common cause, it is only natural to connect it with some other 
phenomenon, one at least just as common to all living cells. 
Moreover, since the phenomena of spontaneous division are 
almost exclusively the prerogative of the living, we must 
therefore, to find their cause, look for something that is es- 
pecially characteristic of the living. Electrical phenomena. 
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which have sometimes been advocated for the explanation of 
cell division, do not meet that last criterion. They are uni- 
versal, found equally frequently in the living and the non- 
living. This, however, does not imply that their secondary 
effects upon cell division may not be of prime importance. 

One phenomenon, common to all growing and dividing 
cells, suggests itself particularly. In every cell there is always 
present a flow of different substances. Metabolism is one of 
the things common to aU living organisms — and almost ex- 
clusively common to them. Various foodstuffs flow into the 
cell; different waste products flow outward. The principal 
agency of this transport of substances is diffusion. When, for 
instance, a substance is produced inside a cell, its concentra- 
tion there exceeds the concentration in the surrounding me- 
dium. Because of the thermal agitation of its molecules, the 
substance is transported from regions of higher concentrations 
to those of lower concentrations. This transport would be very 
rapid if it were not for the resistance offered by the medium, 
in which the substance is dissolved, to the movement of the 
molecules flowing through it. But if the medium exerts a force 
of resistance opposing the flow, then the flowing molecules, by 
virtue of Newton’s Third Law, exert a force on the medium, 
in the direction of flow. Hence, in every living cell, by virtue 
of its ever present metabolism, there always is a more or less 
complex system of forces acting upon different parts of the 
ceU. And it behooves us to investigate whether the magnitude 
and direction of these forces may not be sufficient to produce 
eventually a division of a cell. 

Since these forces are due to diffusion flow, we must, as a 
preliminary study, investigate various general properties of 
such a flow. The direct mathematical manner of studying 
such a problem is to investigate the solutions of the differen- 
tial equations of diffusion, well known in physics. Here, how- 
ever, we meet at once with a serious difficulty. Given the size 
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and shape of a cell, the rates of production or consumption of 
a given substance, the diffusion coefficients, which character- 
ize the resistance of the medium to the flow of these sub- 
stances, and some other physical constants, we can determine 
the resulting flows in principle. But actually, except for the 
very simplest cases, the mathematical difficulties attending 
the solutions of the differential equations are, at present, in- 
superable. 

Moreover, even granting that we could solve exactly the 
differential equations of diffusion for a large number, or even 
for all cases, this would stiU leave us with a rather serious 
handicap. The distribution of the diffusion flows in every in- 
dividual case depends, among other things, on the exact shape 
of the cell. A slight variation of the latter will modify the 
analytical expressions describing the distribution of concen- 
trations and flows. But, since there are no two cells perfectly 
alike, the exact solution of the problem for a given case would 
contain a tremendous amount of detail which is biologically 
insignificant because it applies only to the given case. If we 
wish to draw general conclusions, we shall have to drop most 
of the details and pick out only the essential features of the 
solution, which would hold for a large number of different but 
somewhat similar cases. But with a complicated mathemati- 
cal solution this cannot be done unless we actually make a 
study of a very large number of cases and compare the re- 
sults. This certainly would not be availing ourselves of the 
powers of generalization offered by mathematical analysis. 

Two ways out of this difficulty may be suggested. One is 
to solve the problem exactly for some simple case and to trust 
that the conclusions so reached will hold, with some degree 
of approximation, for those actual and more complex cases 
which do not deviate too much from the ideal simple case 
chosen. This way has been followed in mathematical biology 
hitherto. Inasmuch as the sphere is almost the only shape for 
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which the general differential equations of diffusion can be 
solved exactly, we have limited ourselves to the theoretical 
study of spherical cells. In spite of such a limitation, some 
rather interesting and important results have been obtained 
(see our Mathematical Biophysics [Chicago : University of Chi- 
cago Press, 1938]; hereafter referred to as “MB"). While the 
greatest majority of cells do not possess a spherical shape, still 
a rather large class of cells are almost spherical, and we would 
naturally expect that for such cells our mathematical deduc- 
tions will hold approximately. But the inadequacy, in the long 
run, of such a procedure is obvious, especially where we are 
dealing with the cell division, during which the shape of a cell 
is anything but spherical. We could, by comparing the initial 
stage of a large spherical cell with the end-stage of two spheri- 
cal half-sized cells, demonstrate that the energy relations are 
such that the end-stage is mechanically more stable than the 
initial stage. From that we would argue that at a given size 
the large cell may spontaneously divide in two. But nothing 
could be said about the most important intermediate stages 
and about the mechanism proper of division. 

The second method of approach, which has been outlined 
in the Appendix of our previous book {MB), attacks the diffi- 
culty in a radically different way. We do not attempt to solve 
the mathematical problems exactly for a shape specified in 
detail, but rather we confine ourselves to finding mathemati- 
cal relations between such gross features as are common to all 
cells of a given type, in spite of a difference in detail. We may 
have two cells, with a given approximate ratio of length to 
width. The biologist would describe them as belonging mor- 
phologically to the same class. Yet the detailed shape of their 
surfaces wdU differ enormously from a mathematical point of 
view. It would be quite illusory, and practically unnecessary, 
to take into consideration such small, biologically irrelevant 
details. But, inasmuch as two such cells behave biologically 
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in a similar way, we must look for suck relations as depend 
only on their gross shape and structure. It seems to be more 
promising, from a practical point of view, to deal with rela- 
tions involving only orders of magnitudes, and not with exact 
expressions. 

Let us consider a cell having an over-all approximate length 
2^1 and an over-all width 2^2 (Fig. i). Its internal structure 
may be inhomogeneous; but for the time being we shaU neg- 
lect those inhomogeneities, for in detail they will differ from 
cell to cell. Consider first, for simplicity, that a substance is 



produced inside such a cell, at an average rate q gm cm""^ sec""^, 
a rate which does not depend on the concentration of the sub- 
stance. The actual rate of production will undoubtedly vary 
from point to point within the cell. But the details of such 
variations again constitute individual characteristics. What- 
ever we measure and observe are only average values, and 
only for them may we find some general regularities. 

Similarly, let us consider the average concentration c of 
the substance inside the cell. We shall express concentration 
in gm cm~ 3 . Again neglecting possible detailed individual 
variations, all we can say is that, since the substance diffuses 
outward and flows from center to periphery, therefore the 
actual concentration will be higher around the center and will 
drop toward the periphery. The average value c will be found 
approximately somewhere midway between periphery and cen- 
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ter. Along the periphery the concentration will also not be 
constant. Neglecting detailed variations, let us denote the 
average concentration inside the cell at the periphery along 
the “ends” AB and A'B' (Fig. i) by Ci, and the average 
peripheral concentration along the “sides” A A' and BB' by 
dTj. Then the average drop of concentration per unit length is 



in the longitudinal direction and 



(i) 


(2) 


in the transversal direction. 

If Bi denotes the coefficient of diffusion for the produced 
substance inside the cell, then, according to the fundamental 
law of diffusion, the average flows of the substance in the longi- 
tudinal and transversal directions are equal, respectively, to 

C ^ Cl 

2Di gm sec”^ 

fi 

and (3) 

2Di gm cm ^ sec ^ , 

fz 

When a cell produces a substance which diffuses outward, 
then in the immediate neighborhood of the cell there will be 
an excess of that substance in the external medium. If, as will 
generally be the case, the external medium, which we con- 
sider as extending infinitely, contains in the absence of the 
cell the substance in question in a concentration Co, then in 
the vicinity of the cell the concentration of that substance will 
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be greater than Co. For an oblong cell this external concentra- 
tion near the cell wiU, in general, be different along the “sides” 
A A' and BB' and along the “ends” AB and A'B' (Fig. i). 
Let us denote this concentration at the “ends” by d and at 
the “sides” by d. 

If the cell is bounded by a membrane of some sort, then the 
flow of the substance through this membrane is approximately 
proportional to the difference of concentrations on both sides 
of the membrane, as is known from observations on membrane 
permeabilities. If h denotes the permeability of the mem- 
brane, then the flow per second per square centimeter of the 
membrane is, for the “ends” and “sides,” equal, respective- 
ly, to 


h{ci — d) gm cm~“ sec“' 

and 

Hcz ” ci) gm cm“® sec*"^ . ^ 


( 4 ) 


From expressions (3) and (4) it follows that the dimension of 
the diffusion coefl&cient is lengthytimCj while that of permea- 
bility is length/time. In absolute units, therefore, the diffu- 
sion coefficient is expressed in cm"" sec”^, while the permeability 
is expressed in cm sec"™^ (MB, p. ii). 

The excess of the substance outside the cell becomes less 
and less as we move farther away from the cell. Practically, 
it is appreciable only over a finite range of distance from the 
boundary. We shall denote this distance by d. The distance 6 
is of the order of magnitude of the over-all size of the cell. 
The justification of this assertion will become apparent sub- 
sequently from comparison of our approximate equations with 
some known exact ones. The average drop of concentration 
per unit length outside the cell in its neighborhood is then 
equal, in the longitudinal direction, to (c^ -- Co)/S and, in 
the transversal direction, to (c^ — Co)/d. Denoting by De the 
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diffusion coefficient in the external medium, we have, for the 
corresponding flows, 


and 


— — gm cm““ sec"^ 


~ gm cm~“ sec ^ 
5 


(s) 


The amounts given by expressions (3) arrive per second to 
each square centimeter of the membrane. They must be equal 
to the corresponding amounts passing through every square 
centimeter of the membrane, as given by the expression (4). 
These, in their turn, must be equal to the corresponding 
amounts leaving every square centimeter on the outside of 
the membrane and given by expressions (5). Hence, writing 
down those required equalities, we obtain the following four 
equations: 

2jDi(c — Ci) = rik(ci — c'l) , (6) 

zDiic — Ca) = TzhiCz ~ ^a) > ^7) 

^{c-c.) = ^{ci-Co), (S) 

fi 0 

^ (c - Ca) = ~ (ci - Co) . (9) 

ri 0 


The total area of the “end” surfaces, AJB and A'B', is of the 
order of magnitude of Therefore, the total amount of 
substance leaving through the ends is, because of expressions 
(4) and (6), equal to 

42^M£TliE^gmsec- 


(10) 



DIFFUSION PHENOMENA IN CELLS 


9 


The area of the A A' and BB' being of the order of 

magnitude of 4Trxr2, the total amount leaving through the 
^‘sides’’ is equal to 

^TcDiTxic “ Ca) gm sec“^ . (ii) 

The volume of the cell is approximately If we approxi- 

mate it by a cylinder, we would get 2'Krjrl, which is of the 
same form except for the coefficient. For shapes with rela- 
tively rounded-up ends the first expression offers a somewhat 
better approximation. However, since we are after general 
types of relations rather than after exact numerical values, 
which differ from cell to cell anyway, it does not matter 
which expression we choose. This point shall be further il- 
lustrated below, when we discuss some results. For definite- 
ness we choose the first expression for the volume. Then the 
total amount of substance produced in the cell per second is 
given by 

l-Wag gm sec-'S (12) 

where q is the rate of production per cubic centimeter. The 
total amount of substance present inside the cell is equal to 
iTtrxriG gm, and its rate of change with respect to time is 
equalto 

^ gm sec-' (13) 

This total rate of change is equal to the total amount produced 
per second less the total amount leaving the ceU per second. 
Hence, expression (13) must be equal to expression (12) less 
the sum of expressions (10) and (ii). Writing down this re- 
quirement, shortening both sides by and factoring out 

Di, we find 
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Eliminating ci from equations (6) and (8), we obtain, after 
very simple rearrangements, 

2{DiDe 4 “ dDih)c + rJlDeCo / -v 

2DiJD, + 2Wih + nhDe ’ 

Similarly, eliminating from equations (7) and (9), we find 

2 {DiDe + Wih)c + rJiD^o , 

2 DJD,+ 2 Wih + uhDe ■ 

From (15) and (i6) we have 

2DiDe + 2hDJl + r,hDe ~ 


2 DJ), + 2 hDih + uhD, ~ 


Introducing expressions (17) into equation (14) and using the 
abbre-vdation 


ZJiDiD, 


(2Z)iA+ 2 Wjh^rJiD,) {2DJD ,+ 2 Wih-\-uhD,) 

2 { 2 DiDc\- 2 Wih-\-rihDc)ri+ {2DiD^-\-2Wih-\-r3hDe)ri ’ 


we find that equation (14) may be written thus: 

dc c — Co 

dt~^ A~ ■ 


Equation (19), integrated in the usual way, gives 


c = Co + A? - CAe-'/A , 
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where C is a constant of integration, determined by the initial 
conditions, that is, by the value which I has at a given fixed 
moment. Quite regardless of the value of C, and therefore re- 
gardless of the initial conditions, the last term of equation 
(20) decreases to zero with increasing t. Hence, after a suffi- 
cient lapse of time, c is always given by 

r = Co 4* Ag , (21) 

in which A is given by equation (18) and is independent of 
time. 

When we make ^1 = ^2= ro, which corresponds to a round- 
ed-up shape of the cell, expression (18) for A simplifies con- 
siderably. Making this simplification and introducing the re- 
sult into equation (21), we find 


= Co + 



, 2 ^ 


(22) 


We may compare this expression with the one obtained by 
solving exactly the differential equation of diffusion for a per- 
fect sphere. In this case the average value c is obtained by 
integrating the right-hand side of equation (22) of chapter ii 
of MB over the whole volume of the sphere and then dividing 
it by the volume. We thus find 


c = Co + 


Sh ^ isDi ^ s dJ ^ 


(23) 


We see that equations (22) and (23) are the same form if we 
make 5 = ro ~ ri ^ r2, which justifies the assertion about S 
made on page 7. The coefficient of the first term in paren- 
theses in equation (23) is somewhat larger than in equation 
(22), but the other two coefficients are smaller. 

Another interesting comparison is made by considering the 
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general case when rt > r^, and when both and h are very 
large, so that we may put in equation (i8) De = h= co . 
The expression for A again simplifies, and equation (21) now 
gives 


I rlrl 

$Di 2rl + r 




(24) 


Gale Young'* has computed the exaci expression for the aver- 
age concentration inside an ellipsoid of revolution with semi- 
axes r, and r. for the case De = h= =0 . His expression is 


C = Co + 


$Di 2r\ + 


72 ^. 


( 25 ) 


which again is identical with equation (24) except for a dif- 
ference in the numerical coefficient. 

As we have emphasized above, the individual variations 
from cell to cell are so large that it would be futile to compare 
exact numerical values. But the regularities that are found 
in the behavior of cells, in spite of these individual variations, 
suggest a search for relations of a general form which may re- 
main invariant from cell to cell despite variations of numerical 
values. Here we are dealing with precisely that kind of rela- 
tions. It will also be more clear now as to what we meant on 
page 9 by saying that it would not matter much whether 
we choose one or another approximation for the volume of the 
cell. A different choice in that case would merely modify the 
numerical value of the coefficients in equations (22) and (24) 
without changing their form. 

It is, of course, justifiable, to some extent, to make an ap- 
proximate mathematical study of oblong cells by assuming 
that their properties will be closely enough described by those 
of ellipsoids for which in some cases the mathematical prob- 

* All numbered references are given at the end of each chapter. 
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lem can be solved exactly. But the difference between an actu- 
al ceU and an ellipsoid is so great that we shall be able practi- 
cally to apply only the most general properties of the solution 
anyway, and the mathematical exactness becomes rather su- 
perfluous. If, moreover, we compare the derivation of our 
equation (24) with that given by Gale Youngs for equation 
(25), we shall see at once the labor- and time-saving advan- 
tages of the approximation method. Moreover, it maybe noted 
that, while equation (21), together with equation (18), gives 
us an approximate solution for the most general case, when 
Di, De, and h are finite, the ellipsoidal problem has, so far, been 
solved exactly only iox h= <0. 

We derived here the equation (21) for the stationary state 
as a particular case of the general equation (19). But we can 
just as well derive equation (21) directly, without integrating 
a differential equation, by limiting ourselves, from the outset, 
to stationary states. In such a case we require that the total 
amount of substance produced per second in the cell equals 
the total amount leaving the ceU per second. This requirement 
gives us, instead of equation (14), the following one: 

, = (, 6 ) 

Introducmg into this equation the expressions (17) and using 
again the abbreviation (18), we at once obtain equation (21). 
Thus, for a stationary state the whole problem reduces to the 
solution of five simple linear algebraic equations, namely, 
(6)-(9) and (26). For the general nonstationary case we have 
to integrate a very simple ordmary differential equation (19). 
The exact solution of a diffusion problem requires, in general, 
the handling of partial differential equations. 

If, instead of a substance produced, we consider a substance 
consumed by the cell and flowing into it, we obtain the same 
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expressions, except that we must substitute -g for 5 {MB, 
p. 9). Since A is always positive, equation (21) shows that 
for a produced substance there is always an excess of sub- 
stance inside of the cell, as compared with the outside, while 
for a consumed substance there is a deficiency. Since c cannot 
be negative, therefore for q<o, equation (21) has a meaning 
only for Co > A |g|, where ig| is the absolute value of q (cf. 
MB, chap. iii). 

Equation (20) shows that the smaller the A, the more rap- 
idly the exponential term vanishes and the more rapidly the 
stationary state is reached. The expression A decreases with 
increasing h, D,, and A and increases with increasing and r^. 
This can be seen directly either from expression (18) or, still 
easier, from the special cases leading to expressions (23) and 
(24). Thus, in a sense, the quantity A measures the resistance 
offered to the diffusion flow both by the external and internal 
mediums and by the membrane of the cell. It will therefore be 
appropriately called the total diffusion resistance of the cell. 
It depends both on the physical constants and on the size and 
shape of the cell. 

In the next chapter we shall learn to estimate the value 
of A for some metabolites and shall find it of the order of mag- 
nitude i-io seconds. This shows that the exponential term in 
equation (20) will decrease to about one-third of its initial 
value, which it has for t = 0, when / i to 10 sec. In a next 
similar interval it will decrease to one-ninth of its initial value, 
and so on. Unless the initial value of c at ^ = o is larger by 
several orders of magnitude than the stationary value as given 
by equation (21) — a case which is not likely to occur — this 
stationary value will be practically reached within a minute or 
less. Thus, with the foregoing value of A, a cell kept under 
constant conditions for a matter of hours will be practically 
in a stationary state. 

It is, however, very important to emphasize that a sta- 



DIFFUSION PHENOMENA IN CELLS 


15 


tionary state does not always exist. We have considered 
hitherto only cases in which the rate of production or con- 
sumption of a substance by the cell does not depend on the 
concentration of that substance within the cell. If the rate q 
is itself a function of c, so that we may write q = qic)i then, 
for the stationary state we still have five algebraic equations, 
namely, (6)-(9) and (26). But, putting in equation (26) q = 
q(c) makes this equation, in general, nonlinear in c; and under 
these conditions, for some types of functions q{c), the whole 
system of five equations may have no real roots, or the roots 
may be real but negative, which is physically impossible. This 
means that under such conditions stationary states do not 
exist. A study of the general nonstationary case in some such 
instances reveals that the concentration c increases indefi- 
nitely with time, regardless of initial conditions. Such is the 
case when a substance is produced inside the cell at a rate 
proportional to its concentration, so that q = kH, where is 
a constant. It is interesting that in this case a stationary state 
exists only for cells of sufficiently small size {MB, chap. iv). 
For a rounded-up cell — that is, for = Tz = ro — and for De = 
CO the critical size r*, above which no stationary state exists, 
is given^ by 



and by a somewhat more complicated expression for the gen- 
eral case of a finite De* However, it is interesting that, when 
q — — k'^c — ^in other words, when the substance is consumed at a 
rate proportional to its own concentration — stationary state 
does always exist and is approached regardless of initial con- 
ditions. 

A study of the foregoing case was made first by solving the 
differential equations of diffusion exactly for a sphere {MB, 
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chap. iv). Not only is the mathematics somewhat involved, 
but the expression for r* could not be obtained explicitly, be- 
ing given by a transcendental equation. The approximate 
method shows here also the advantage of both greater sim- 
plicity and generality. 

It is also possible that, while a stationary state exists, it is 
not approached asymptotically. The concentration c in this 
case oscillates periodically around the stationary value. For a 
very special case {Di = De= «>), such a situation was stud- 
ied by the author (MB, chap. vi). Alvin Weinberg^ worked 
out subsequently an exact solution for a more general case, 
namely, = co , -with Di and h finite. The mathematical 
procedure is again rather elaborate and lengthy. Such peri- 
odical solutions occur when we have at least two substances 
in the ceU, affecting the rates of formation of each other. 
Weinberg’s exact solution has been obtained only for the par- 
ticular case that the ratio h/Di is the same for both sub- 
stances. Using the approximation method, Weinberg'* finds in 
a rather simple way expressions of the same form as the exact 
ones, and without usmg the foregoing restriction as to the 
ratio h/Di. The comparison of the two papers of Weinberg is 
rather instructive. 

For use in subsequent chapters we shall derive, in conclusion, a 
few auxiliary equations. 

From equations (15) and (16) we have 

^ ^ 2 (DiPe -h 5Dih)hDe{ri — r^) (c — Cp) 

' ' l2DiDe+ 28 Dih+r,hDe)( 2 DiD,+ 2dDi/i-hr,hDe ) ' 

Introducing expressions (18) into (21) and the latter into (28) 
we find ’ 


^2 Cj — 


r£^ 

3 


^ 2(I>e+SA)(ri— fa) 

3(2DiDe-h2SD{h+rjhD,)ri+(2DiDe+2SDih+r,hDe)r:, 


(29) 
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For A = CO , equations (28) and (29) reduce, respectively, to 


and 


_ _ 28(^1 - , 

3 2 ( 25 Di + rxDe>, + ( 25 Di + r,De)r/ 


From equations (6), (7), and (17) we have 


, 2DtZ)e 

" 2DiDe + 25 Dih + rJiDe 

and ‘ 

, 2DiDe /_ X 

- C2 - 22)^, + 2 dDih + rJiDe • 


(32) 


For A = 00 this gives Ci — ci ~ C2 — ci - 0. Physically, this 
should be the case, since, when the membrane is infinitely perme- 
able, any finite concentration difference on both sides will result in 
an infinite flow. Hence, for any finite flow the concentration dif- 
ference in this case is zero. 
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CHAPTER II 


COUPLED REACTIONS: APPLICATIONS TO 
CELL RESPIRATION 

We shall now consider a somewhat more complex situation, 
in which several coupled chemical reactions take place in the 
cell simultaneously. Actually, a very large number of simul- 
taneous reactions are taking place in any cell, and we cannot 
hope to tackle the problem m all its complexity. We shall, 
however, investigate such relatively simple types of reactions 
as are suggested by actual observations. 

In a rather large class of cells we have roughly the following 
situation.^- ^ The cell consumes some type of sugar from the 
surrounding liquid, and it also consumes oxygen. The sugar, 
upon entering the cell, undergoes a series of breakdown proc- 
esses, during which a number of mtermediate products are 
formed. One of the most conspicuous of these mtermediate 
products is lactic acid, which frequently leaves the ceU in 
appreciable quantities. Within the cell, however, the lactic 
acid undergoes further changes, as a result of which it com- 
bines with o.xygen, giving carbon dioxide and water. It is also 
known that the formation of lactic acid from sugar is partly 
reversible, a part of the lactic acid being resynthesized into 
higher sugar. 

Following H. D. Landahl,^ we shall therefore investigate 
mathematically the following situation ; sugar flowing into the 
cell and there breaking up into lactic acid, part of the lactic 
acid resynthesizing back into sugar, part leaving the cell, and 
the rest combining with oxygen to produce carbon dioxide and 
water, both of which leave the cell. To make our problem defi- 

i8 



CELL RESPIRATION 


19 


nite, we shall assume that the sugar involved is glucose, 
CeHnOe. Furthermore, we shall assume that the breakdown 
of glucose into lactic acid, according to the equation 


C6H..06 = 2C3H6O3 , (i) 

goes on approximately at a rate proportional to its average 
concentration, Ci, in the cell, while the resynthesis of lactic 
acid is proportional to the square of the latter’s concentration, 
c^, As for the rate of oxidation of lactic acid, we shall make 
the plausible assumption that it is proportional to the 
product, ^2^3, of the average concentration C2 of lactic acid 
in the ceU and of the average oxygen concentration in 
the cell. The oxidation reaction takes place according to the 
formula 


C3H6O3 + 3O2 = 3CO2 + 3H2O . (2) 

For each gram-mol of lactic acid, 3 gm>mol of oxygen are re- 
quired. Hence, 90 gm of lactic acid require for their oxidation 
96 gm of oxygen; or, for each gram of oxygen consumed, 
0.94 gm of lactic acid is oxidized. 

These assumptions do not, of course, represent accurately 
the actual situation, but they may be used as JSxst approxima- 
tions. Subsequently we shall discuss other, more complex as- 
sumptions, which are closer to the true nature of cellular reac- 
tions. 

If, as we have done in the preceding chapter, we denote the 
rate of production of a substance by a positive quantity, then 
the rates of consumption are negative. Denoting the average 
rate of consumption of glucose by gm cm~^ sec“^, the 
average rate of production of lactic acid by ^2 gm sec"^, 
and the average rate of consumption of oxygen by — ^3 gm 
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cm ~3 sec“'; denoting, further, by a, h, and/ positive coeffi- 
cients of proportionality and putting n = 0.94, we have'* 


“f" ^^2 J 

(3) 


(4) 

?3 “ /^2^3 • 

(5) 


The amount of carbon dioxide produced per cubic centimeter 
per second is equal, according to equation (2), to = 
— 1.38J3 gm, while the amount of water produced per cubic 
centimeter per second is equal to = — o.sdga gm. Remem- 
ber that and q^ are negative, while q^ and gs are positive. 
The quantity q^ may be either positive or negative, depending 
on whether more lactic acid is produced from sugar than is 
oxidized, or vice versa. In the latter case, if there is a suffi- 
cient supply of lactic acid outside, it may flow into the cell, 
to be oxidized there. 

We shall confine ourselves to the three reactions represent- 
ed by equations (3), (4), and (5). For, when we have calcu- 
lated the values of Ci, Cj, and Cj and thus know q^, q^, and q^, 
we immediately have q^ and gj. For carbon dioxide the prob- 
lem thus reduces to one of a substance produced at a given 
rate, such as has been treated in chapter i. Water, being in- 
compressible, has practically the same concentration always, 
Cs = I gm cm~2; and its rate of production determines its 
velocity of flow from the cell. It may, at first, appear essential 
to modify oxir considerations on diffusion so as to take into 
account that the flowing water also carries part of the dis- 
solved substances with it. The latter are transported not only 
by diffusion but also by convection. However, a very general 
estimate {MB, chap, i) shows that, because of the usually very 
small concentrations of the different metabolites (io~3 — 
io ~4 gm cm- 3 ), the transport due to water flow is quite negli- 
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gible, as compared with the transport due to diffusion, and 
may be safely neglected. 

We shall consider here a stationary diffusion state, assum- 
ing that it exists under the conditions specified above. This 
assumption is borne out by finding that the equations have 
real positive solutions, as we shall see presently. 

Denoting by Coi the external constant concentration of the 
Ith. substance, we find, by the same line of argument as used 
in chapter i, that for each substance the following equation 
holds in the stationary state: 


h = Ooi + kiqi (6) 

The indices I have one of the values i, 2, 3 For Az we 


have expression (18) of chapter i, in which we shall write 
hi, Du, and D^i, referring to the permeabilities and diffusion 
coefficients of the different substances. Thus, hi denotes the 
permeability for glucose; h, that for lactic acid; etc. 

Putting in equation (6) / = i, 2, and 3 and substituting it 
into equations (3)~(s), we find 



II 

1 

+ Aig'i) + b{co2 + 

(7) 


qz = nq^ - 


( 8 ) 


?3 ” —/{Cq-z 

+ A 2 g 2 )(co 3 + . 

(9) 

Putting 




A 

^hKl, 

H = /A 2 A 3 , 


B 

Q.1)Cq2^2 j 

K = 1 H-ZcosAa , 


C 

= (zAi + I , 

AT * ^Cq^Al2 “ ildro ^'03 ) 

( 10 ) 

D 

acoi ” ^^01 

2 ) S = fCoiCo^ = 5^0^03 ) 


ic 

= /A 2 , 

0 ^ 

11 
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we may write equations (7), (8), and (9) thus: 

Aq% + Bq2 - Cqx - D = o , (li) 

+ 5-2 - rags = o > (^^2) 

Eq.q, -\-Kqi + Mq2 + S = o. (13) 


These are three equations for the three unknown quantities 
q^, q^, and q^. If we solve these equations, we can immediately 
find Ci, Cl, and C3 by using equation (6). But equations (11), 
(12), and (13) are, except for the different convention regard- 
ing the signs of the g’s, the same as equations (22), (23), and 
(24) of chapter v of MB, obtained previously by H. D. Lan- 
dahl by a different method, and their solution is the same as 
given there. While, however, the equations discussed in MB 
held only for spherical cells, the equations given here and ob- 
tained by the general approximation method hold for cells 
of any oblong shape. 

Equations (ii), (12), and (13) determine qi, qi, and q^ in 
terms of such constants as the diffusion coefficients, permea- 
bilities, and the external concentrations Coi, c^i, and C03. We 
shall confine ourselves here to the study of the average oxygen 
consumption q^ in its dependence on the external oxygen con- 
centration C03, as there are a number of experimental data 
with which our conclusions in this respect may be compared. 

Solving equations (ii), (12), and (13) for q^, we find that 
for C03 = <» , qi tends to a constant limiting value, qf, given 
by (MB, p. 48) : 


_ (aAi + i)coj +_aA.iCo- 
nAi{aAi i) 


Considering, instead of qi, the quantity y = qi/q* (that is, 
the relative oxygen consumption) and putting C03 = x, we 
find from equations (ii), (12), and (13) that, approximately, 

x = ^y+ 

I -y 


(IS) 



CELL RESPIRATION 


23 


where 


For 3; = ij rjc = while for 3; = o, = o. That is, if we 
plot y against we obtain a curve rising from zero and ap- 
proaching asymptotically the constant value i. Such should 
be the relation between the relative oxygen consumption y 
and the external oxygen concentration all other conditions 



Fig. 2. — The curve represents the theoretical equation (15). The 
points represent experimental data.s A3 = 2 sec; £>3 > 8 X io“7 cm^ 
sec“%* ^3 > 4.2 X 10“^ cm sec“*. 

being kept constant. To what extent equation (15) is verified 
experimentally is shown in Figures 2-7. 

Inasmuch as the approximation method used here applies 
to any shape of cell, it should apply also to pieces of tissue 
consisting of a large number of cells. The diffusion coefficient 
for such a piece of tissue is, of course, only a sort of average 
quantity, but the same holds for any cell that is physically 
heterogeneous. Figure 7 illustrates how the equations here 
developed apply to tissues. 




3.— The curve represents the theoretical equation (15). The 
points represent the experimental data. A3 = 4.3 X lo"'^ sec; A > 1.3 X 
10“^ cm» sec""^; > 3.6 X 10“^ cm sec’^ 



Fig. 4. — ^The curve represents the theoretical equation (15). The 
points represent experimental data ,7 A3 = 4 sec; Dj > 4 X io “7 cm^ 
sec~^; hi> 2 X io~4 cm sec~^ 
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In some of the experiments, namely, those with Arbacia 
eggs, the conditions were, however, not such as assumed in 
the derivation of equation (15). In the experiments with Ar- 
bacia eggs the latters were kept in saline solutions, practically 


P.S.TANG AND C.S FRENCH - CHLORELLA CELLS 



Fig. 5. — The curves represent the theoretical equation (15). The 
points represent experimental data.® Observations were made at different 
temperatures, as indicated on the graphs (degrees centigrade). A differ- 
ent scale for x is used for x > 0.75 X lo”® gm cm”^, in order to represent 
the whole range of data in the same figure. The constants calculated for 
the two different temperatures are as follows: For 10° C: f == 4 X 10”® 
gm cm“3; ^ = 1,5 X io “7 gm cm‘'3; A3 = 0.33 sec; > 1.4 X io~® 
cm^ sec“^; > 1.4 X io“4 cm sec~^. For i7?5 C: f = 3.7 X 10“® gm 
cm“3; = 4.2 X gm cm“ 3 ; A3 = 0.42 sec; £>3 > i.i X 10 ■’®cm=* sec“^; 

/Z3 > 1. 1 X cm sec“^ 

free from sugar. The eggs in this case oxidize their own re- 
serve of sugar stored within the cell. This reserve is apparent- 
ly large enough to maintain for a sufBcient time a constant 
supply and thus secure a quasi-stationary state. The theory 
for this case requires, however, only a very slight modifica- 
tion. All we have to do® is to substitute in equations (3) and 
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(4) a constant, a^, for the term acj. That constant, a^, repre- 
sents the constant rate of decomposition of the sugar reserve. 
However, this modification of the equations does not alter 
the final equation (15) (MB, chap. v). 


P. S.TANC and C.S, FRENCH - CHLORELUA CELLS 



Fig. 6. — Same as Fig. 3, with constants as follows; For 8° C; f = 
4 X io~* gm cm“; f = 1.5 X lo~^ gm cm Aj = 0.48 sec; Dj > i X 
io~® cm’ sec~’; hi > 0.9 X xo'"'* cm sec“’. For 20° C: f = 3.3 X io~® 
gm cm~ 3 ; f = 5 X icT’’ gm cm“X; A3 = 0.43 sec; Dj > i.i X io~* 
cm’ sec~q ^3 ± i.i X io~< cm sec~'. 

From the comparison of the experimental data with equa- 
tion (15) we obtain numerical values for f and Since in 
the experiment the quantity is given as the absolute maxi- 
mum rate of oxygen consumption per cubic centimeter, knowl- 
edge of I gives us the value of A3, the total diffusion resistance 
for oxygen. The value of A,, unfortunately, cannot be ob- 
tained from the knowdedge of f unless we know the coefifi- 
cient /. 

In experiments made with Arbacia eggs and with Chlorella 
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cells the suspension of the latter was thoroughly stirred in a 
respirometer. Such stirring uniformizes the concentrations of 
substances outside the cell and is equivalent in its effects to 
making A very large. In the expression for A3 we may there- 



Fig. 7. — Comparison of tlie theoretical equation (15), as represented 
by the curve, with experimental data, taken from a paper by P. S. Tang, 9 
Since the size of the fragments were not given, the constant A3 could not 
be computed. 


fore put, with a good approximation, = 00 . Moreover, 
since the Arbacia eggs are spherical, we have fi = ^2 = n* In 
this case the expression (18) of chapter i gives 




(17) 


While To is given experimentally, the knowledge of A3 does 
not yet give us either or Di^ separately. Since, however, 
both terms on the right-hand side of equation (17) are posi- 
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tive, therefore at least one of them, or both, must be of the 
order of magnitude of A3. Hence, we have 

L»i3>4- and (18) 

9A3 9A3 

where at least one or both of the equivalence (~) signs must 
hold. That is, if rl/gki, then certainly ~ aro/qAj; 

and, conversely, if A3 » 2V9-A-3. then rl/gk^. But 

both equivalence signs may hold simultaneously. 

Values of the different constants estimated in this way are 
given in the legends accompanying Figures 2-6. 

In principle A3 and Di^ could be obtained by experiment- 
ing on two different groups of cells of the same type but of dif- 
ferent size. This would supply us with two equations of the 
tj-pe of equation (17), one for each value of ro. The objection, 
however, could be raised that cells of different size may have 
different values of A3 and Di^. The equations developed here 
may, however, be applied not only to single cells but to slices 
of tissues of different size and shape. The expression for A3 
will, of course, now be of a more complex structure, though 
we could stm have Pe = 0= . One might hope to separate, in 
this way, A3 from Di^ However, in the case of a multicellular 
slice the inner meshwork of the membranes would add to the 
diffusion resistance due to Pij, and we would have to modify 
somewhat the equations to take this into account. 

Data by P. S. Tang and C. S. French® represented on Fig- 
ures 5 and 6 were taken for different temperatures. Determin- 
ing for each temperature the parameters f and we may 
study their variation with temperature, as is shown on Fig- 
ure 8. As we see, the quantity f increases with temperature. 
The exact variation cannot be ascertained from only four 
points. But, since f is proportional to the rate of reaction q* 
(cf. Eqs. [16]), we would exnect f to increase with tempera- 
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ture, since all chemical reactions do so. For J we should ex- ^ 
pect a decrease, since the reaction constant / would be ex- 
pected to increase. 

In many cell-kinetic problems, however, only the total dif- 
fusion resistance A of the cell is of consequence, and it is there- 
fore important to be able to determine these quantities, even 
if we cannot resolve them into their components. 



Fig. 8. — ^Variation of the quantities q^i (fuE line), ^ (broken line), and 
I" (alternate line) with temperature, as computed from Figs. 5 and 6. The 
scale for q^i is in 10“'^ gm cm“3 sec“^; that for in lo”? gm cm-^j and that 
for in 10”^ gm cm”^. 

The quantity q* is, as seen from equation (14), a function 
of the external concentration Cq2 of lactic acid and of the con- 
centration Coz of glucose. Unfortunately, there are no data 
available to check this relation. Should, however, such data 
be obtained and should equation (14) be found to be suffi- 
ciently accurate, this would immediately give us the values of 
Ai and A2, the total diffusion resistances for glucose and lactic 
acid. Other parameters may be determined by studying the 
consumption of glucose and the output of lactic acid in their 
dependence on the external concentration of glucose, lactic 
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acid, and os>"gen. The urgent need of further experimental 
work on cell respiration under conditions just mentioned is 
thus clearly indicated. 

If we deal with cells having a very thin membrane, such 
that its physical constitution is not likely to be very different 
from that of the rest of the cytoplasm, then the contribution 
of such a membrane to the total diffusion resistance is negligi- 
ble, and we may, with a good approximation, make in equa- 
tion (17) h= which gives us Di^ in terms of Aj. In many 
experimental determinations of cell permeability, especially 
in those using plasmolysis of plant ceils, it is the permeability 
of the relatively thick layer of cytoplasm plus surface mem- 
brane that is determined. What we actually get there is, there- 
fore, a quantity proportional to i/A. 

H. D. Landahl has also studied more complex situations. 
For instance, we may consider-t the case where the resynthesis 
of sugar occurs in two routes — one direct and proportional to 
Icl, and the other coupled to the energy liberated by the oxi- 
dation of part of the lactic acid and to the rate/caCj of oxygen 
consumption. As far as the dependence of oxygen consump- 
tion on ox\"gen pressure is concerned, we obtain equation (15). 

Some interesting results are obtained by considering the 
following type of reaction.-* 

a g (2) 

Hesose -f Phosphate :± He.tosediphosphate ± Lactic acid -|- Phosphate 
(i) (6) (7) + (6) 

Oxygen 

Carbon dioxide + Water 

(4) (s) 

In this scheme the number under or above each metabolite 
shows the index which refers to it in the following equations. 
Thus, Cj means the average concentration of hexosediphos- 
phate; q&, the average rate of production per cubic centimeter 
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of phosphate, etc. The letters a, g, and / next to the arrows 
are the coefficients of proportionality of the corresponding 
reactions, used in the equations which follow. Let us assume 
that the phosphate never leaves the cell, in any form, and 
that it neither accumulates nor is depleted, so that the 
total amount of phosphate per cubic centimeter is a constant, 
Cp. For simplicity v/e shall consider the reactions as not re- 
versible. Also, again the equations for carbon dioxide and 
water shall be omitted, for the same reasons as before. De- 
noting by Ml the molecular weight of the Zth substance and 
putting m = 2MtlM-^ and n = 2M2/6M3, we have for this 
case, instead of equations (3)-(s), the following set of equa- 
tions : 


= 


(19) 

~ 

-nfCiCi -f gcj , 

(20) 

?3 = 

fCiC} , 

(21) 

q 6 = 

— q^ = — OC1C6 -f gc, = 0 , 

(22) 

Cp “ 

C 6 -f mC'j . 

(23) 


Consider first the anaerobic case, that is, the case when 
there is no external supply of oxygen, so that C03 = = o. 

Then, solving equation (23) for Cj, substituting into equation 
(22), and rearranging, we find 


C6 


gCp 

amci + g ' 


(24) 


Substituting into equation (19) for the value given by 
equation (24) and remembering that Ci = Cai + Aiji, we ob- 
tain 


amA^ql + (amcn + agCpAi -f g)qz + agCoiCp — o . (25) 
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Putting 


„* - _I£2 ■ 

m ’ 




r' — . 

^ ~ m ’ 


r = , (26) 

am 


we may wnrite equation (25) thus: 

Cqz “ ^ 2 " 1 “ 


^2 
I — 


(27) 


Equation (27) gives a relation between the external concen- 
tration of sugar and the rate of sugar consumption. In ab- 
sence of oxygen the rate of output of lactic acid is proportional 
to the rate of sugar consumed. Therefore equation (27) also 
gives us a relation between rate of lactic-acid output and sugar 
concentration. Equation (27) is of the same form as equation 
(15), showing that for zero sugar concentration the lactic-acid 
output is zero and that the latter tends to a constant limiting 
value for large sugar concentrations. 

Data on this point are very meager; still, the following con- 
siderations may be of some interest. On Figure 9 data are 
plotted for lactic-acid output from experiments with glucose 
and fructose. If we use different scales for the abscissae, then 
the two sets of points seem to faU on the same curve, represent- 
ed by equation (27). Whether this is so or not, cannot be 
decided, because of a too small number of points. But assum- 
ing, for a moment, that this is so, we are led to the following 
rather interesting considerations. If the curves did coincide 
in the figure, then = 2of'i and tb = 2of 'j, where the sub- 
scripts refer to fructose and glucose and where the factor 20 
arises because there happens to be a factor of 20 between the 
two scales. Since the limiting values of the lactic-acid produc- 
tion is nearly the same in each,'" then by equation (27) we 
have Aif^ = 2oA.i^, or the total diffusion resistance for fruc- 
tose is twenty times that for glucose. If we assume, also, that 
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the secondary stage is, to a large extent, the same in both 
cases, so that Cp is nearly the same in each, then, since m is 
the same for both, g must be the same for both to satisfy 
equation (27). In other words, the second-stage reaction goes 
on at about the same rate in both. In this case one might 



Fig. 9. — The curve represents equation (27). The points represent 
experimental data. For explanation of double scale on the axis of ab- 
scissae see text; p. 32, 

then tentatively assume the secondary reactions to be largely 
equivalent. From equation (26) we have = 2oaf^, so that 
the rate of formation of the second product, hexosediphos- 
phate, from glucose is twenty times the rate of its formation 
from fructose. However, it could very well be that, at least 
in some instances, the two curves could not possibly be made 
to coincide by a linear change in the abscissa scale. In this 
case it would not be possible to compare the total diffusion 
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resistance or other constants directly, but it would be neces- 
sary to fit the curves empirically to evaluate and f The 
quantity Ai would then be given; but, in general, an addi- 
tional empirical relation would be required to evaluate a, g, 
and Cp, 

These considerations indicate the importance of accurate 
determinations of anaerobic lactic-acid production in its de- 
pendence on the concentration of different sugars. 

For the dependence of oxygen consumption on oxygen pres- 
sure the present case leads again to equation (15). The same 
equation is again obtained by considering a still more com- 
plex s^'stem,^ invohdng the presence of a catalyst necessary 
for the oxidation of lactic acid. All this shows that the oxy- 
gen-consumption-oxygen-pressure curves are determined by 
the gross features of the phenomena and do not depend much 
on intermediate details. The reason for the curves being such 
as shown in Figures 2-7 is physically this : If the concentra- 
tion of lactic acid in the cell were constant, then increasing in- 
definitely the oxygen concentration would increase also the 
oxygen consumption indefinitely. But, as the oxygen concen- 
tration, and hence the oxygen consumption, increases, so also 
does the consumption of lactic acid increase, and therefore the 
concentration of the latter in the ceU decreases. Since the oxy- 
gen consumption is proportional to the product of the concen- 
trations of lactic acid and of oxygen, therefore with increas- 
ing ox^^gen concentration and decreasing lactic-acid concen- 
tration the oxygen consumption increases less and less rapidly, 
approaching a limiting value, which is determined by the 
maximum output of lactic acid due to the breakdown of glu- 
cose. 

In the next chapter we shall use an important constant of the 
cell, which we call the “glycolytic coefficient,’’ This coefficient 
is defined in the following way: 3/3 is the ratio of molecules of sugar 
glycolyzed into lactic acid to those oxidized completely. If i mole- 
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cule of glucose is oxidized completely, then 3/3 are glycolyzed; and 
altogether 3/3 + i molecules of sugar are consumed. For the oxi- 
dation of I molecule of sugar, 6 molecules of oxygen are necessary, 
while 6 j 3 molecules of lactic acid and 6 molecules of carbon dioxide 
are produced. Hence, /3 is equal to the ratio of the number of mole- 
cules of lactic acid produced to the number of molecules of oxygen 
consumed. In other words, 


-g2 ^ 

Mz qz znqz ’ 


Mz 

^ = -77- = 0.94 . 
3^3 


(28) 


Since ^3 < o, /3 is positive when > o, that is, lactic acid leaves 
the cell. 

Solving equation (13) for using notation (10), and substitut- 
ing the value so obtained into equation (28), we find 

^ Znjkzqzicoz + k^qz) ^ 


Putting, again, y = q^/qf and using expression (16) for f and we 
have, from equation (15), 

c „3 + = c „3 - fy == 7:7^ . (30) 

Introducing equation (30) into equation (29) and using equations 
(16) and (14), we find 


/3 = 


(aki + i) 

3[{aki + i)coa + acoikz] 


Q/Cojkz I y 

jikx + I y 



(31) 


which gives us the glycolytic coefficient in terms of Coi, Coa, and c„s 
(through y, by means of equation [15]). For small values of exter- 
nal lactic-acid concentration Cozj equation (31) reduces to 
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and for large values of external oxygen concentration Cos and there- 
fore y I, equation (32), because of equation (30), reduces to 

^ ■” 3W/A2C03 ‘ 

For suiSciently large values of external oxygen concentration C03, 
the ghxolytic coefficient approaches zero, and may even become 
negative if ^02 is sufficiently large. The glycolytic coefficient is in- 
versely proportional to C03. For a constant and not too high value 
of Cos, 0 increases with decreasing total diffusion resistance, A2; and 
to a lesser extent /3 increases as the external glucose concentration, 
Coi, increases; as the total diffusion resistance to oxygen, A3, in- 
creases, and as the total diffusion resistance to glucose, Ai, de- 
creases. Also, jS decreases with / and, to a lesser extent, increases 
with a. 
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CHAPTER in 


DIFFUSION FORCES AS A POSSIBLE 
CAUSE OF CELL DIVISION 

As we have remarked in chapter i, a flow of diffusing sub- 
stance must exert a force on the medium in which the diffusion 
takes place, the direction of the force being the same as the 
direction of flow. The exact theory of these forces presents 
rather great difficulties because it involves a detailed knowl- 
edge of the moleciflar structure of liquids, especially of colloids 
and gels, in which the diffusion usually takes place. A pure 
solvent, like water, may behave differently, in this respect, 
from a complex gel. In order not to lose ourselves in generali- 
ties, it is well to make at once, at least in a preliminary way, 
some rough picture of the structure of the inside of the cell. 
We shall first adopt the simplest possible scheme, later on dis- 
cussing possible complications. 

Obviously, if the cell consisted merely of a sac, formed by 
the membrane and containing only the same solvent and 
solutes as are present outside, such a cell would not exhibit 
the simplest fundamental properties of actual cells. While, for 
instance, sugar and oxygen are present in sufficient concentra- 
tions, no appreciable oxidation occurs outside the cell. Other 
important reactions also take place only inside the cell, though 
the necessary components may be present outside. This must 
be due to the presence inside the ceU of various catalysts, to 
which the cell membrane is not permeable, so that they are 
confined within the cell. Most biological catalysts, or en- 
zymes, are known to be attached to various colloidal struc- 
tures.^ We may picture, therefore, the inside of the ceU as 
being composed of water, in which are suspended small col- 
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loidal particles. Water is the common medium outside and 
inside the cell. The cell proper consists, then, according to 
this picture, of an aggregate of colloidal particles, inclosed in a 
membrane. 

The drag forces, exerted by the flow of diffusing substances, 
will be acting both on the water and on the particles. We are 
interested especially in the forces acting on the particles, 
which represent, in a way, the “skeleton” of the cell. 

Considering that the diffusion flow always occurs as a re- 
sult of nonuniformities of concentrations, and that in a field 
of nonuniform concentration of a solute any particle, as well as 
any molecule, of the solvent is subject to an asymmetric bom- 
bardment by molecular collision of the nonuniformly distrib- 
uted diff usin g solute, the author derived {MB, chap, vii) the 
following expression for the force F exerted on a spherical 
particle of volume Vp: 

F = ^ aVp grad c . (i) 

2 M 

Here R denotes the gas constant per mol (R = 0.83 X 10® 
ergs degree"'); T, the absolute temperature; M, the molecular 
weight of the diffusing substance; and a, a constant, less than 
unity but of the order of unity, and determined by the molec- 
ular structure of the solvent. The symbol grad c (read “gra- 
dient of c”) means the change of the concentration c of the 
diffusing substance per unit length, in the direction of the 
greatest change. The minus sign (— ) shows that the force 
is directed from higher to smaller concentrations, that is, in 
the direction of flow. 

Gale Youngs pointed out certain inadequacies in the au- 
thor’s derivation of equation (i). He derives a similar equa- 
tion by considering, instead of the molecular bombardment, 
the viscous drag exerted by the flowing substance upon a 
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spherical particle. He obtains a formula very similar to equa- 
tion (i) except that the quantity —1RTY/2M is not multi- 
plied by a correction factor a but has a very smaU additive 
correction term. Inasmuch as both equations are derived un- 
der rather simplified assumptions, they cannot be expected to 
give more than a close order of magnitude; and no choice is 
indicated, at present, between them. It is, however, impor- 
tant that two entirely different approaches lead to essentially 
the same expressions. 

If we have iVj, particles per unit volume, then the force 
acting on a cubic centimeter of the solid cell constituents is 
approximately equal to 

/= grade. (2) 

But VpFp is nothing but the actual volume occupied by the 
particles contained in a cubic centimeter, a quantity which we 
shall denote by \i. We thus have for the force acting on the 
cell constituents per unit volume 

/ = -f ^ grad c . (3) 

Actually the colloidal particles constituting the body of the cell 
may not be spherical. From general considerations it appears 
plausible that for any shape of particles an equation similar to 
equation (i) will hold approximately, except for a different 
numerical coefficient^ (also MB, chap. vii). We may now im- 
agine the particles to be held together by highly viscous or 
even solid strands, thus forming a quasi-solid meshwork 
through which water and substances dissolved in the latter 
circulate freely. This conforms in general with accepted ideas 
on the structure of some jellies. In that case an equation of 
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the form of (2) will still represent the force exerted by the dif- 
fusion flow on the unit volume of the cell. The coefficient ix 
now denotes the relative true volume occupied by the solid 
structure. When we shall speak of the deformation of the cell, 
we shall refer to the deformation of this structure, consisting 
either of a jelly-like rigid network or of disconnected solid par- 
ticles. Water itself, which may either be stationary or circu- 
late, is considered just as a special metabolite. 

We may now proceed to calculate the mechanical effects of 
such forces on the cell. We must remember that each metabo- 
lite produces a force of its own and that the total force is the 
sum of the individual components due to each metabolite. Let 
us first investigate the effect of a single substance and then 
pass to the more complex case, which will enable us to com- 
pare with experimental data the equations derived. 

Let us again consider a substance produced in the cell at a 
rate q gm cm"“3 sec“^. The flow in this case is directed out- 
ward, and the forces are also directed outward. They have, 
therefore, the tendency to disrupt the cell. For a case of a 
perfect sphere, due to the complete symmetry of everything 
with respect to the center, the resultant force will be zero. But 
any deviation from the spherical shape will result in an asym- 
metry in the distribution of the forces ; the resultant effect will 
not be zero, and a deformation of some sort will occur {MB, 
chap. Lx). 

In the case of a cell producing a substance we, first of all, 
have a force, given by equation (3), acting on each element of 
volume. Moreover, at the membrane of the cell there is a 
pressure acting on the membrane because of the difference in 
the concentrations on both sides of the membrane. This dif- 
ference in concentrations results in a difference of osmotic 
pressures. The osmotic pressure of a solution of concentration 
c is equal to 



CELL DIVISION 


41 


The net pressure on the membrane is therefore the difference 
of the pressure due to the inside concentration, directed out- 
ward, and of the pressure due to the outside concentration, di- 
rected inward. For the “ends’^ and the ^^ides’^ of the cell 
(chap, i, Fig. i) this resulting pressme is equal, respectively, to 

RT 

= - cO dynes cm-= ; 

RT 

^ (cj — cO d3Ties cm~-^ . 


Whether we regard the protoplasm constituting the cell as 
a sol or as a semisolid gel, the cell must be considered as a body 
possessing an appreciable amount of viscosity. Any deforma- 
tion resulting in relative displacement of the different parts 
of the cell with respect to each other will set up viscous re- 
sistance forces which will be stronger the higher the rate of 
deformation. 

In order to calculate the deformation of the cell under the 
influence of the forces discussed above, we shall make use of 
a theorem by Betti, 3 combined with the laws of plastic fiow.^ 
If Iz is the length of the body in a given direction (say in the 
direction of the s;-axis), and if V is the volume of the body, 
S its surface, 77 its viscosity, X, F, and Z the components of 
the volume force per unit volume in the direction of the cor- 
responding coordinate axes, and X,,, F^, and Zv the compo- 
nents of pressure at the surface, then, according to Betties 
theorem, the average relative rate of change {ill^{dlzfdi) of 
Iz is given, for a body of any shape, by^ 


I ^ 


+ /J [zZ, - hiyY. + , 


(6) 
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where the first integral is extended over the whole volume, 
and the second over the whole surface of the body. 

Equation (3) shows that the volume force / can be derived 
from a potential (p = ^{RTp./M)c. In this case, as has been 
shown by Gale Young, ^ the volume integral in equation (6) 
can be easily transformed into a surface integral of the form 

jjc{s cos {v, s) - |[a: cos {v, x) + y cos (v, . 

Equations similar to (6) hold for all three directions. If we 
fiuid that the body elongates in two directions — that is, if for 
two directions the right-hand side of equation (6) is positive — 
then we shall fiind that it contracts in the third direction, the 
right-hand side of equation (6) being, for this direction, nega- 
tive. If the body possesses approximately axial symmetry, 
as we assume to be the case for a cell of the type considered in 
chapter i, then elongation in one direction results in a con- 
striction in the other two directions. 

Putting the s-axis in the direction of the largest dimension 
of the cell, we have k = n. Putting 

1*1 = /Jc{s cos {v, z) — ^[x cos (vj ac) + y cos (i^, y)] (7) 

and 

7. = [zZ, - lixX, + yY,)]dS ( 8 ) 

where F„ and are the components of the pressure given 
by equation (5), and using the above-mentioned transforma- 
tion of the volume integral, we have 

= 4--JL r 

U dt 2 Mr]V ^'^ ' 


(9) 
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We now proceed to calculate the quantities and They 
require a knowledge of the concentration at the surface of the 
ceU. Strictly speaking, when the shape of the cell changes 
with respect to time, we cannot take for the concentration in- 
side the cell at any given moment the values obtained for the 
stationary state of the cell having the same shape as the given 
cell has at this moment, for the latter were obtained in 
chapter i under the assumption that the shape of the cell does 
not vary. However, we have seen that, beginning with any 
initial state, the stationary state is practically reached within 
a couple of minutes. If the deformation of the cell proceeds 
very slowly, so that its shape does not change much within 
those few minutes, then we can apply with sufficient accuracy 
the expression derived for the stationary state. In many cases 
the time it takes a cell to divide is of the order of magnitude 
of hours. The time necessary to readjust the diffusion state 
to a different shape being much shorter, the system may be 
considered as a succession of stationary states. However, the 
more general case, when the rates of deformation are not slow 
enough, must also be studied. This will be done in chapter v. 

In calculating Ji we remark that on the “sides” of the cell 
cos {v, z) = o, while x cos (v, x) y cos {v, y) = Ti. Remem- 
bering, from chapter i, that the volume V of the cell is equal 
to 

V = iirr^rl , (lo) 

we find as a contribution to Ii for the “sides” 


— C2 X ir2 X znTTi X 2fi = - fVc, . (ii) 

In a similar way we find that the contribution of the “ends” 
to Ii is f Fci. Hence, 


h = |F(c, — C2) . (12) 
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In calculating I2 we have along the “sides” 

RT 

xXp "h yXv — J'afe cQ j 
while at the “ends” 

■n'p 

zZp = ^ fife - cO . 

Therefore, 

L = f Y ^ [(ci - 4 ) - (c2 - c',)] (13) 

Introducing equations (12) and (13) into equation (9) and 
introducing into that the expressions (28) and (32) of chapter 
i, we find 

I dn _ RT 

T-, dt 2 M 7 ] 

[3iJi.dh+ (sfi— 2)Di]hDiDe(ri—r2) (c— Cp) ' 

(2Z>iZ?e+ 2BDih-\- TihDe) (2DiDe+ 2 8Dih+ nhD^ ' 

For fj = Ti the deformation drx/dt vanishes, as should be the 
case. 

For q> 0, that is, when the substance is produced, 
c — Co > o. The quantity (3^ — 2) is usually negative. Since 
r, — > o, therefore for sufficiently small 

dri/dt < o for g > o. If, however, the quantities 8 
are sufficiently large, then the numerator becomes positive, 
and dfi/dt > o for g > o. In other words, omitting the rather 
unlikely case that ii> %, for cells of a sufficiently small size 
the effect of the diffusion forces generated by a produced sub- 
stance is such that any oblong shape tends to round up. But 
when the size of the ceU is sufficiently large, the effect of such 
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forces is to elongate the cell. Hence, a substance produced by 
the cell will always cause an elongation of the cell as soon as 
the latter grows sufficiently large. To insure positivity of 
we must have 

5 > 1 io"“^ cm . (15) 

If we consider a consumed substance, we have ? < o, and 
the situation is reversed. The diffusion forces tend to elon- 
gate the cell only when it is small enough. As we shall 
presently see, the surface tension of the cell always opposes 
the elongation. But the effect of surface tension decreases 
with increasing size of the cell Therefore, for q> that is, 
for a produced substance, for large sizes the diffusion forces 
will prevail, and the cell will elongate; for small sizes, however, 
it will not elongate. For 5 < o, that is, for a consumed sub- 
stance, for large sizes both diffusion forces and surface tension 
oppose elongation; for small sizes diffusion forces tend to 
elongate the cell, but the surface tension opposes it; for very 
small sizes the surface tension prevails. Only in a limited 
range of sizes can a cell elongate for q <0, and that only for a 
special choice of constants. These results are obtained by an 
exact treatment of the problem of stability of a spherical cell 
(MB, chap. ix). They can be obtained in a similar way by the 
use of the present approximate method. 

If h is very large, so that we can put in equation (14) 
}i ^ CO ^ then that equation reduces to 

i. ^ hDiDejri — r^{c — cd) / 

Tx dt 2M7} { 2 dDi + riDe)( 252 >i + 

In this case drx/dt is always positive for g > o and negative 
for g < o. For simplicity we shall confine ourselves to this 
case in the following discussion (cf. chap, ii, p. 30). 
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Making A = «> in equation (i8) of chapter i, introducing 
the expression for A into equation (21) of chapter i, and intro- 
ducing the expression for c — Co, so obtained, into equation 
(16), we find 

I dr, _ RTnuqti. (ri — . . 

fi dt aMtj 2{2dDi -t- riDe)rt -f {25Di + 

Putting S = r^, we have 

1 ^ (ri — r^r^rl 

fi dt 2 Mt] 2 { 2 rJ)i -p rJO^rt 4 - {^Di -f- D^rl ' ^ 

Gale Young-f worked out exactly the mathematical prob- 
lem of the elongation of an ellipsoid of revolution under the 
influence of diffusion forces, for the case h = <0 , assuming 
that, as it elongates, it still remains exactly an ellipsoid. The 
equation which he obtained is rather similar to equation (18). 
An actual comparison shows that both equations give rather 
similar numerical results within a range of elongation that 
is practically important. For very large ratios the two 
equations differ widely. Although Young’s formula, as we 
have said, is obtained by an elaborate exact solution of the 
problem, no preference can be attached to it because it holds 
exactly only for a perfect ellipsoid, a situation which is not 
met biologically. It is important, however, that within a 
practically significant range of elongations Young’s exact 
results agree with our approximate ones, lending thus a 
further justification to the method used here. 

Hitherto we have considered only the effect of the diffusion 
forces. But any cell is subject also to forces produced by sur- 
face tension. As is well known, these forces always tend to 
make the cell spherical, because the sphere has the smallest 
surface for a given volume. Thus it can be seen at once that 
the surface-tension forces will oppose any elongation of the 
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cell. The quantitative calculation of the effects of surface 
tension on the rate of elongation may be made approximately 
in the following way. 

At any point on the surface the surface tension y results 
in a pressure, directed inward and equal to® 

where Ri and R 2 are the principal radii of curvature. For a 
cell Uke the one represented on Figure i, we have approxi- 
mately at the “ends,” or “poles” 



while on the “sides” we have approximately 


(20) 


Ri = ri Ri = ~ Z '^7 ■ 

We thus have at the ends a surface force equal to 


and at the sides a force 


27 

(22) 

u 

-{ k + k }- 

(23) 


Those forces contribute to the expression /j, given by equation 
(8). We now have, on the “sides,” 

»X. + ,F.--rr.(i + i) 



and, at the ^^ends,^^ 
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By an argument similar to that which led to equation (13) 
we find the contribution of the surface tension to the quantity 
I 2 to be of the form — iVy(rj: — ^2)7^2. Therefore, accord- 
ing to equation (9), the contribution of the surface tension to 
the relative rate of elongation is 


X 

2r} Va 


(24) 


Expression (24) can be obtained also by a different argument 
{MB, Appendix). The same expression is obtained exactly for 
the case when the shape of the cell is an ellipsoid of revolu- 
tion.^ 

Adding expression (24) to the right-hand side of equation 
(18), we now obtain the following complete equation of 
elongation: 

I dfi _ RTqii (fi — r^r^rl 

fi dt 2M7} 2 { 2 r 2 Di + r^De)rz + { 2 Di + De)rl 


27] rjr2 * J 

By adding expression (24) to the more general equation 
(14) we find a corresponding more general and more compli- 
cated equation for the total relative rate of elongation. For 
the present we shall limit ourselves to the case of very large 
permeability, so that we can put A = 00 in equation (14), 
which leads to equation (18) and through it to equation (25X 
In order that elongation should take place at all, we must 
have dr^/dt > o. If an originally spheroidal cell should begin 
to elongate spontaneously, the coefi&cient of — ra must be 
positive for = ra = Under these conditions the slightest 
accidental deviation from the perfect spherical shape will be- 
come more and more enhanced. Putting = 7*2 = r 0 in the 
right-hand side of equation (25), we find that the expression 
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thus obtained is negative for small values of and positive 
for large ones. This means that a spheroidal ceU, producing a 
substance at a constant rate q gm cm~3 sec“*, wiU not elongate 
if its size ro is less than a critical value r*. If, because of 
growth, it reaches a size r* or exceeds it, it begins to elongate 
spontaneously according to equation (25). The value r* of ro, 
at which that expression just changes sign, is obtained by 
putting it equal to zero and solving with respect to r*. We 
thus obtain an expression for the critical size of a spheroidal 
cell, above which it elongates spontaneously. The expression 
thus obtained is 

r* = 3 UMi 2 Di + D,)y 
\ RTu^q 

A ceU does not necessarily need to increase in size in order to 
elongate. Its size ro may remain constant; but the different 
quantities, such as Di, De, 7, and q, which determine the 
critical radius r*, may vary. For a particular set of values of 
these quantities we may have To < r*, and the ceU does not 
elongate. For another set of values ro may become greater 
than r*, and the ceU begins to elongate. 

When the elongation begins, ri increases while ro decreases. 
If the growth of the ceU is sufficiently slow, so that the volume 
of the cell does not vary appreciably during the process of 
elongation, then we may express ro in terms of r^ and of the 
volume V from equation (10). Introducing the value of r, so 
obtained into equation (25), we obtain an equation containing 
ri only. When r-^ becomes very large and ro very small, so 
that ^2, then, as is readily seen from equation (25), the 
first term of the right-hand side varies like i/ri, while the 
second term varies like (i/n) Vri . Hence, denoting by A 
and B two constants, for large values of r^, dri/dt varies like 
A — Br\l^. But this expression is negative for sufficiently 
large values of r^- This shows that for sufficiently large elonga- 
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tions the right-hand side of equation (25) becomes equal to 
zero* For the value of ft, at which this happens, and for the 
corresponding value of the diffusion forces tending to 
elongate the cell are just equal and opposite to the surface- 
tension forces. The elongation therefore, once begun, pro- 
ceeds only to a certain finite value. It will be shown in chapter 
iv that, in some cases, when for a spheroidal cell the critical 
size r* is exceeded even by the smallest amount, the elonga- 
tion which sets in proceeds to a finite extent, determined by 
the foregoing considerations. That is, for = Ta = ro> 
the quantity dft/dt is positive and remains positive until the 
ratio fi/Va exceeds a definite value.^ 



Once the limiting elongation is reached, dr^/dt becomes 
zero. But Betti’s formula, used in derivation of equation (25), 
gives us the average rate of elongation for the cell as a whole. 
It may be readily seen that the elongation of the different 
parts of the cell will be different. In fact, at the equator the 
cell is pulled apart in two opposite directions by forces acting 
on every element of volume of each half of the cell (Fig. 10). 
On the other hand, any cross-section of the cell near the end 
is acted upon only by the force applied to each element of 
volume of the relatively smaller part of the ceU. Thus we 
should expect the middle of the ceU to elongate more rapidly 
than the ends. But, since the elongation is accompanied by a 
constriction in the perpendicular direction, the cell constricts 
more rapidly in the middle than at the ends and thus becomes 
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dumbbell-shaped. Even when the average elongation has 
stopped after reaching its limiting value, the cell may con- 
tinue to elongate and to constrict in the middle, at the same 
time shortening and thickening somewhat at the ends. The 
actual distribution of the elongations of different parts of the 
cell depends on the exact shape of the latter. For the case of 
an ellipsoid of revolution the problem has been treated by 
Gale Young.® The theory of the process of constriction in its 
early stages is still wanting. For sufficiently advanced stages 
an approximate treatment has been given by Gale Young.® 
When the constriction has proceeded sufficiently far, so 
that the cell has approximately the shape represented in 



Figure ii, we may estimate the further elongation of the cell 
by considering it as two spheres, connected by a neck and 
repelling each other. The repulsion is due to the fact that 
each sphere is producing a substance which flows outward, 
a force resulting which acts on each element of the other 
sphere in a direction away from the center of the first. 

The force per unit volume at any point outside the sphere 
and due to that sphere is proportional to the total rate of pro- 
duction, that is, to q/2, and is inversely proportional to the 
square of the distance from the center of the sphere to that 
point and to the external diffusion coefficient. The force must 
be multiplied by the volume of the other sphere. Therefore, 
the whole force is directly proportional to qrf and inversely 
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proportional to a‘Dc. The approximate expression for this 
force F, derived by Gale Yoimg,® is 




(27) 


This force may be considered as applied to the total surface 
of each end of the neck, so that per square centimeter 
there is a force F/tt^ directed along the line of centers out- 
ward. Besides, there is a lateral surface tension force tending 
to hold each end onto the neck, equal to 2rrj, which is again 
equivalent to a force 


dynes cni“^ , (28) 

applied to each square centimeter of the ends of the neck 
directed along the line of centers inward. On the lateral sur- 
face of the neck we have a capillary pressure of about 

— - dynes cm~“ , (29) 

directed inward. As the spheres tend to pull apart, the nar- 
rowest region on the neck constricts inward, much like the 
“necking” of an iron rod in a testing-machine. Introducing 
the above-mentioned forces into Betti’s equation and applying 
the latter to the “end” and “sides” of the neck alone, and 
ignoring the internal diffusion forces inside the neck, which 
in this case are rather small compared with the force F, we 
find, for the relative rate of change of the length I of the neck, 

1 dl _ X F — Try 
I dt ^TT) 


(30) 
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For a viscous incompressible body the relative lateral con- 
striction, due to elongation, is half the relative elongation, so 
that we have 

1 ~ = —1 — 
r dt~ 2ldf 


Putting 

p = 2- ■ o = 

6v ’ ^ s^vDer^'^M ’ 


(32) 


and introducing equation (31) into equation (30), we find 



(33) 


If Q/r^ P, which, with plausible values of the constants, 
is frequently the case, then 



(34) 


In any case, the condition Q > Pr is necessary for a con- 
striction to proceed at all. If this inequality is not satisfied, 
the cell will not constrict. Integrating equation (34) and de- 
noting by / the value of r at ^ = o, we find 

r" = - 2Qt . (35) 


This equation shows that r becomes zero; in other words, 
the two half-cells separate completely at 



(36) 


All the formulae derived in this chapter were obtained on 
the assumption that a single substance is produced in the cell. 
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Actually, however, a very large number of substances are 
produced and consumed in every ceU. The flow of the pro- 
duced substances results in forces which tend to divide a cell. 
The flow of the consumed substances gives rise to forces in 
the opposite direction. In order to know whether the ceU will 
divide or not, we must know the direction of the resultant of 
aU these forces. With the very large number of metabolites 
involved in any ceU, it is out of the question to calculate this 
resultant force exactly. However, some approximate estimate 
can be made. 

From the discussions of chapter i — in particular from equa- 
tion (23) — it foUows that for A = co , the average gradient of 
concentration in a spheroidal ceU is approximately propor- 
tional to qrJDi. Since the average force per unit volume is 
RTfi/M times the gradient, it foUows, therefore, that this 
force is proportional to gjMDi. The same thing foUows for an 
oblong ceU from equation (24) of chapter i. If we have a 
large number of metaboUtes in the ceU, those with a larger 
qjMDi wiU contribute most to the mechanical forces. As it 
happens, a group of reactions in the cell usuaUy stands out 
because of its particularly high rate q. This group involves 
the respiratory reactions discussed in chapter ii. For those 
reactions (consumption of sugar and oxygen, production of 
lactic acid and carbon dioxide) the values of q are about io“® 
gm cm“5 sec“’^ (MB, chap. x). The other important reactions, 
such as are involved in protein breakdown, are of the order 
of io “9 gm cm~s sec“^ The molecular weights of the sub- 
stances involved in the respiratory reactions are relatively 
low, while the molecular weights of the different amino acids 
are higher. This makes q/M larger for the respiratory reac- 
tions. To some extent this is upset by the circumstance that 
we would, in general, see the heavier molecules diffuse more 
slowly and thus have a smaUer Z?». However, this, as shall be 
discussed presently, is not necessarily the case. At any rate. 
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because of the large difference in q, it seems justifiable to dis- 
cuss the case, as a first and very rough approximation, as if 
only the respiratory reactions were contributing to the 
mechanical forces. Using the same notations as in chapter ii, 
denoting by q^ the rate of production of carbon dioxide, and 
omitting possible mechanical effects of flowing water, for 
reasons discussed in MB, chapter x, we find that the sum of 
the forces, owing to the four above-mentioned reactions of the 
respiratory group, is proportional to the quantity 

-2- = - + -2?- j ^ 

The quantity (37) can be expressed in terms of qJM-i, 
that is, the molar rate of the consumption of oxygen, and of 
the glycolytic coefficient / 3 , introduced in chapter ii (p. 34). 
From the definition of i 3 it follows that 


iL 




_2i = —3±. a ■ h. = 2l 

M.. ^ ’ M, M, • 


(38) 


Introducing this into equation (37) and putting 


we find 



(39) 


(40) 


In a previous discussion {MB, chap, x) we made the as- 
sumption that the diffusion coefficients are inversely propor- 
tional to the cubic roots of the molecular weights. In this case 
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it is found that A > o; B < o. Thus, for small values of /3 
the expression q/MD, is negative, showing that the forces 
due to consumed substances prevail and that the resultant 
force is directed inward. The assumption D ^ is^ how- 
ever, not always justified. Some indirect evidence^ indicates 
that the diffusion coefficient for lactic acid is very much 
smaller than that for oxygen, although the size of the molecule 
is not so much larger. The following explanation of this may 
be suggested. The diffusion of different metabolites takes 
place in the cytoplasm, which is a complex colloidal medium. 
Some metabolite molecules, especially such as those of lactic 
acid, which have an appreciable dipole moment, may be 
mostly adsorbed to much larger coUoidal particles. Such a 
molecule is then transported in the diffusion flow not as a 
free molecule but partly as an adsorption compound: Col- 
loidal micella-|-Metabolite molecule. The diffusion coefficient 
in this case will be determined largely by the size of the col- 
loidal micella and may therefore be very small. Two mole- 
cules may be of the same size, but one may be almost com- 
pletely adsorbed and the other may be in a free state. Thus 
they may have very different diffusion coefficients. A pre- 
liminary study of such a picture has been made by John M. 
Reiner.® 

Since lactic acid is more polar than glucose, it is to be ex- 
pected that Bt < Di, so that ^ > o, since < o. In that 
case the resulting force increases in the positive (outward) 
direction with increasing /3. For B < o the force is directed 
outward only when jS > —B/A. For B > o the force is al- 
ways directed outward but increases in intensity with increas- 
ing /3. Thus a ceU should divide the more readily and the 
process of division should be the more rapid, the greater the 
glycolytic coefficient jS. This, in a general way, suggests a 
connection with the weU-known findings of 0. Warburg® that 
abnormally rapidly multiplying tumor cells have a much 
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larger glycolytic coefl&cient than normally multiplying cells. 
While for the latter ^ is of the order of o.i, for the former it 
may be as large as 4. Inasmuch as Warburg^s findings refer 
also to the rate of growth, we shall discuss this important 
question more in detail in chapter iv. 

If an equation would involve g, My and Di only in the com- 
bination qJMDi, then all that we would have to do is to sub- 
stitute for that combination the expression + B. How- 
ever, the diffusion coefficients also enter separately in a num- 
ber of formulae. Inasmuch as we are dealing here with very 
rough approximations, we may substitute for M, Di, and D^ 
some suitably chosen average values, M and De, 

and put for q an expression of the form 

q^ AP + B . (41) 

The coefficients A and B are approximately of the order of 
AMD and BMD, respectively. Since A and B are both pro- 
portional to ^3, as seen from equation (39), the average q in 
our equations will also be proportional to q^, the rate of 
oxygen consumption; and in general both will be of the same 
order of magnitude. Except for the case where B < 0 and jS is 
just about equal to —BjAy the value of q will also be of the 
same order of magnitude as that of q^. We shall call the 
quantity q the “effective rate of metabolism.’’ 

With these approximations we may proceed to a compari- 
son of different equations, derived in this chapter, with ex- 
perimental data. First of all, we shall discuss equation (26). 
The work of Newton Harvey, Kenneth Cole, and others 
indicates^® that 7 ~ i erg cm”"^. For room temperature 
RT ~ 2.4 X ergs. Taking, as possible values for De and 
Di, io“7 cm^ sec“^ (cf. chap, ii), ix = o.i, g~ 10“^ gm cm ”3 
sec"^ and M ~ 100, we find the value of r* to be of the 
order of io“3 cm. These values are of the same order of mag- 
nitude as the average values of actual cells. It should be noted 
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that the relatively large variations in the values of D, q, y, 
and other constants entering into equation (26) result in rela- 
tively much smaller variations of r*. This throws some light 
on the understanding of the circumstance that exceptionally 
large or exceptionally small ceUs are found relatively seldom 
in spite of the tremendous variations in the physicochemical 
constitution of different cells. 

Equation (26) may he compared with observations in a still 
different way. Since we must introduce into it for q the 
quantity q, defined by equation (41), and since, as we have 
seen, q is proportional to the oxygen-consumption rate, there- 
fore, other parameters being kept constant, the size r* of the 
ceU should vary as the inverse cubic root of the oxygen con- 
sumption. If we could find a number of cells, differing from 
each other only with respect to size and to oxygen consump- 
tion, and if we would plot the former against the latter on a 
logarithmic scale, we should es^ject that the points will fall 
on a straight line with the slope — It is, however, impossible 
to find such an ideal case. Actually, other parameters char- 
acterizing the cell, such as the diffusion coefficients, surface 
tension, etc., will vary from case to case. Since there is no 
expected correlation between these constants, we should now 
expect that, in plotting cell sizes against their oxygen con- 
sumption on a logarithmic scale, the points will not fall on any 
line at all but will cluster along a straight line of slope — J. 
For a proper check of this expectation we should have at least 
a few dozen points, so that they would form a definite cluster. 
Unfortunately, the data available are very meager, and only 
very few experimental points are available. These, however, 
show the proper trend, as seen from Figure 12. The legend 
to this figure gives the sources of the data. For cells like 
those of the heart tissue, which are very oblong, the average 
size r* was computed from their volume V, defining r* by 
Ittt* = V. This is justified by equation (26), inasmuch as we 
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might have plotted the volume against q and would then 
expect a simple inverse proportionality. More exact measure- 
ment on average cell sizes and on their oxygen-consumption 



Fig. 12. — Radii fo of different cells plotted against q, the rate of oxygen 
consumption in gr cm“3 sec”^ on the logarithmic scale. The straight line 
has a slope — Values for q were taken from W. H, Bayliss, Principles 
of General Physiology (4th ed.), p. 612 (London: Longmans, Green & Co., 
1927). Values for n were taken from Stohr’s Lehrhwh der Eistologie, ed. 
W. von MoUendorf (22d ed.) pp. 84, 220, 283 (Jena: G. Fischer, 1930). 
For heart cells, an average Tq, corresponding to the average volume of the 
cell, is used. For values for Arhacia eggs and luminous bacteria, cf. 
Landahl .3 

rates are urgently needed and would be of immense value for 
the development of the theory. 

The comparison of equation (25) with experimental data 
presents similar dififfculties because of a practically complete 
absence of necessary relevant observations. 

K. Belar"^ has recorded the length of the dividing nucleus 
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in spermatocytes of a snail for different times. A dividing 
nucleus undergoes a large number of structural changes; and 
phenomena such as movements of chromosomes, appearance 
of spindle fibers and others, which are not included in the 
present theory, take place in it. Therefore, it might appear 
rather doubtful as to whether there is any sense in comparing 
our equations with such data as Belaf’s. On the other hand, 
however, as we have emphasized in chapter i, the approxima- 
tion method, by introducing average quantities and neglecting 
local structural variations, leads us to relations which should 
be independent of those local structural details. We would 
therefore expect that the gross features of a dividing nucleus 
and of a dividing cell as a whole would be represented ap- 
proximately by the same equations. Thus it is of interest to 
compare equation (25) with Belaf’s data. 

One thing, however, must be kept in mind when making 
such a comparison. Equation (25) is derived with the as- 
sumption that the viscous resistance and the surface tension 
are the only two forces opposing the diffusion forces. In a 
nucleus it appears highly probable that the elastic forces of 
the stretched spindle fibers are also an important contributing 
factor. These elastic forces can easily be taken into account 
{MB, p. 313) and result in the addition of an extra term to the 
right-hand side of equation (25). Although our knowledge of 
nuclear structure in mitoses is exceedingly meager, we may 
expect that these elastic forces will manifest themselves only 
in the earlier stages of the nuclear elongation, previous to the 
moment of complete separation of the daughter-chromosomes 
in the equatorial plate. According to a picture of the mecha- 
nism of mitosis suggested by the author {MB, chap, xiii), the 
elastic spindle fibers are connected at their ends to the chromo- 
somes, on one hand, and to the body of the cell near the poles, 
on the other. As the nucleus elongates, these fibers are 
stretched as long as the chromosomes remain in the equatorial 



CELL DIVISION 


6i 


plate. As soon as they separate, the fibers contract freely, 
moving the chromosomes to the poles. At this phase, how- 
ever, when the fibers contract freely, the additional term in 



Fig. 13. — ^Length, ri, of a dividing nucleus plotted against time, t. 
Points represent observations by K. BSlaf.^^ The curve to the right of the 
arrow represents the theoretical equation derived by G. Young .4 A 
graph of equation (44) is practically the same. At the point indicated by 
the arrow, the chromosomes separate. See test, p. 61. 

equation (25), due to elastic forces, vanishes. Thus, on the 
basis of this picture or a similar one we may expect a discon- 
tinuity in the rate of elongation, just at the moment when the 
chromosomes separate. Up to that moment the elongation 
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curve should be described by an equation containing the addi- 
tional elastic term. Beyond that moment the elongation 
should be described by equation (25). 



Fig. 14. — Same as Fig. 13 but for another set of data, by K. BSlaf” 

In considering the whole cell, such a discontinuity is not 
to be expected, owing to an absence of rigid connections of the 
nucleus with the ceU. In fact, the dividing spindle may be 
moved inside the ceU relatively easily. 

Although these conclusions are by no means binding, it 
is of interest to note that Belaf’s data actually show a dis- 
continuity, as indicated by the arrow on Figures 13 and 14. 
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The inaccuracy of Belaf ’s data makes any definite conclusions 
impossible. We have already mentioned above (p. 46) that 
Gale Young-t has derived an equation for the elongation of a 
cell, assuming that the latter has, during the elongation, the 
shape of an ellipsoid of revolution and that his equation is 
practically equivalent to our equation (18). By adding the 
effect of the surface tension, he therefore arrives at an equa- 
tion equivalent to our equation (25). He has compared his 
equation with the experimental points lying to the right of the 
discontinuity in Figures 13 and 14. The graph of equation 
(25) is practically the same as that shown in Figures 13 and 
14. 

For an actual comparison of equation (25) with any data, 
that equation is better written in a different form. By re- 
arranging the denominator of the first terms of the right-hand 
side, equation (25) can be written thus: 


hill — (ri — 

dt ~ 2M7] 2 De{rl -1- |f^) + 

_ y_ fi - r, 
27} rjTi 

Introducing the abbreviations 



a = 


3f . 

4T ’ 



_ RTixqa _ 
4Mr,D, ’ 



(43) 


multiplying equation (42) by rl / (ti — fj), and rearranging, we 
finally obtain 


ri 


Ti ^ 

— Ti dt 


I -}- 2(1 “b b) “3 ”f" b 


Va 


. (44) 


Since a can be directly measured, the foregoing equation has 
three parameters: K, L, and b. 
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Figures 15 and 17 illustrate an application of equation 
(44) to the case of the division of a cell as a whole. The points 
on Figures 15 and 17 represent the total length of dividing 
isolated Arhacia blastomeres, taken from a micromotion film 
made by Professor Robert Chambers. The broken lines repre- 
sent the theoretical curves, obtained by graphically integrat- 
ing equation (44). For each point, n was measured directly, 



Fig. 15. — ^Length of a dividing isolated blastomere from the second 
cleavage of an Arhacia egg, plotted against time. The points are taken 
from a micromotion film by Professor Robert Chambers. The full line is 
the empirical curve drawn thi*ough the points. The broken line represents 
the theoretical equation (44). Notice that the origin of coordinates in this 
figure corresponds to a length of 3 (arbitraiy units). The relative error 
of the theoretical curve is therefore rather small. 

and Ti was computed from the volmne, which was determined 
from the initial radius of the cell, just prior to the beginning 
of the elongation, by using equation (lo). Figures i6 and 
18 represent a similar comparison of equation (3 5), which 
gives the width of the constriction as a function of time. That 
equation, by the method of its derivation, should hold only 
when the constriction has proceeded sufficiently. We should 
therefore expect greater discrepancies between theory and 
experiment for smaller times. This is actually the case, as seen 
from Figirres i6 and i8. 





Fig. 1 6 . — ^Width of constriction of a dividing isolated blastomere of an 
Arbacia egg, plotted against time. The points were taken from the same 
film and for the same ceil as those of Fig. 15. The curve represents the 
theoretical equation (35). 



Fig. 17. — Same as Fig. 15 but for another cell, photographed on the 
same film. 
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In principle, one could attempt to calculate some of the 
constants of the cell from such curves as those discussed here. 
Unfortnaately, no exact record either of the magnification or 
of the speeds was made by Dr. Chambers in preparing the 
film, since the latter was not made with the present purpose in 
view. Therefore, we do not attempt here to calculate the abso- 
lute values of the constants K, L, and b. As a matter of fact. 



Fig. i8. — ^Same as Fig. i6 but for still another ceU, photographed on 
the same film. (Same ceU as used in Fig. 17.) 

the whole fihn was not taken at regular intervals. However, 
the circumstance that the points, measured from every 
twentieth frame, do fall rather well on a smooth line, as seen 
on Figures 15-18, indicates that the average speed was rela- 
tively constant. 

The foUowing point is of interest. A look at Figure 1 1 shows 
that, approximately, r'' = o.Sto and r'/ = i.6ro. The com- 
parison of equations (32) and (43) shows, then, that the con- 
stant Q in equation (35) is connected to the constant K in 
equation (44) by the relation 

Q n’ 


( 45 ) 
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The magnification of the fihn used is not known exactly; but, 
since for a whole Arbacia egg the value of the radius is ap- 
proximately 3.8 X io "3 cm, the radius of the two blastomeres 
after the first cleavage is about 3 X io “5 cm. Since the film 
used represented the second cleavage, the value of Tq can be 
taken as approximately 3 X 10“"^ cm. Introduced into equa- 
tion (45), this gives K/Q = 29 X lo^. The actual value of 
K/Q for the cell corresponding to Figures 15 and 16 
is found to be 24 X lo^, while for the cell corresponding to 
Figures 17 and 18 it is 18 X lo^. In view of the inaccuracy of 
the data, on one hand, and of the simplifications introduced 
in the derivation of equation (45), on the other hand, nothing 
more than an agreement of the orders of magnitude can be 
expected. 

While the examples given in Figures 13--18 can be con- 
sidered only as illustrations from which no definite con- 
clusions can be drawn, they make quite obvious the im- 
portance of further exact quantitative measurements of vari- 
ous aspects of cell division, made under properly controlled 
conditions. 

We have discussed in this chapter the effect of diffusion 
forces due to metabolism. As we see, those forces which by 
their nature are always present in any metabolizing cell, will 
produce a number of phenomena characteristic of cell divi- 
sion. It appears very likely that in these forces we have found 
the general cause of cell division, for we may say that, even if 
no other forces were present, cells would divide imder proper 
conditions, because the diffusion forces are always there. 
However, as remarked in chapter i, a number of other forces 
and other factors may strongly influence and modify the 
actual process of cell division. Only a thorough mathematical 
study of all possible factors will eventually enable us to decide 
their relative importance. Of particular interest are, of course, 
possible electrical forces. N. Rashevsky“ has pointed out 
that, when ions of opposite charges are metabolized in a cell. 
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then, because of the differences in their diffusion coefficients, 
the distribution of concentrations of the ions will be different 
and thus will give rise to the appearance of electric space 
charges. However, the exact treatment of this problem pre- 
sents great mathematical difficulties. Using the present ap- 
proximation method, Robert R. Williamson^s has solved the 
problem, confirming a surmise made previously {MB, p. 125) 
that, in general, the mechanical forces, due to such electric 
charges, are very small, as compared to the diffusion forces 
discussed here. 
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CHAPTER IV 
CELLULAR GROWTH 

In the previous chapter we have seen that a hi gh glycolytic 
coeihcient should, in general, increase the dividing forces in a 
cell and thus make the latter divide more readily. As has been 
pointed out, this agrees with the findings of 0. Warburg^ and 
others that cancer cells are characterized by an abnormally 
high glycolytic coefficient. But one of the conspicuous char- 
acteristics of tumors is their extremely rapid growth, and this 
is something that is not obviously connected with readier 
division. We shall therefore investigate in this chapter, some- 
what more in detail, the manner in which the glycolytic co- 
efficient may affect the rate of growth. And to do this we shall 
have to discuss the biophysical aspects of growth in general, 
a question which, so far, has been left outside the scope of our 
studies. 

The increase of the mass and of the size of a ceU is due to 
the formation of new substances, composing the body of the 
cell, from some simpler constituents, present in the surround- 
ing medium and diffusing into the cell. It is therefore to be ex- 
pected that the variation of the mass and size of the cell with 
respect to time wiU depend largely on the type of physico- 
chemical reactions which we assume to be taking place within 
the cell. If these reactions are in some way connected with 
glycolysis, then the glycolytic coefficient would naturally play 
a part in determining the rate of growth. It might, therefore, 
be of interest to study such types of physicochemical re- 
actions, and this would lead us into biochemistry rather than 
into biophysics. While such a study might prove to be of 
great importance and should be made, at present we shall 
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limit ourselves to a different aspect of the problem, an aspect 
which emphasizes more strongly the biophysical side and 
which shows possible indirect effects of glycolysis on the rate 
of growth of tissues. Although, unfortunately, no direct 
quantitative experimental data are available for comparison 
with the equations which we shall derive, we shall give some 
theoretically predicted curves which may stimulate the ex- 
perimenter to carry out some measurements. Such measure- 
ments, by either confirming or refuting our expectations, will 
serv^e to decide whether effects of the type here discussed are 
actually occurring in some cells. 

The simplest picture which we may make of the mecha- 
nism of formation of new cell material is that a number of sub- 
stances penetrate by diffusion into the cell and combine there 
to form more complex molecules, constituting the cell body, 
the rate of consumption of every substance being propor- 
tional to the products of the concentrations of cell substances. 
To be more specific, we may, for instance, visualize each sub- 
stance as an amino acid. These amino acids combine to form 
complex proteins. Actually, of course, the same amino acids 
may combine in different ways to form quite different pro- 
teins. But, for simplicity, we shall at present disregard these 
complications. They suggest, however, definite new types of 
problems, and the discussions of this chapter wiU perforce 
remain incomplete until those more complicated cases will be 
studied by a similar mathematical technique. 

Let there be n such substances, their average concentra- 
tions inside the cell being Ci, Ca, . . . . , Cn, and their concen- 
trations outside the cell being Coi, C02, . . . . , Con- To simplify 
matters, without much loss of generality, we shall consider a 
rounded-up, quasi-spherical ceU; and we shall assume all ex- 
ternal diffusion coefficients De as very large, putting • 

The inner diffusion coefficient for the Ith. substance shall be 
denoted by D^i, For each of the substances ci we shall have 
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an equation of the form of equation (21) of chapter i, in 
which, however, we have, since the substances are consumed^ 

qi — — klCiC2 . ... Cn, (i) 

where ki is a positive coefficient. 

Let us put 

kiCiC2C^ .... • ... Cfi ^ ki , (2) 

The left-hand side of equation (2) contains the product of all 
^^s except Cl. We now may write equation (i) thus: 

- kci . (3) 

As before, we may again apply the equation for a stationary 
state if the rate of growth of the cell is sufficiently small, so 
that the cell does not alter its size appreciably during the time 
it takes for the diffusion state to approach stationariness. 
This may be done even though the diffusing substances con- 
sidered here contribute directly to growth. 

Equation (21) of chapter i, re-written in the notations 
used now, reads 

Cl = Col + Mqi . ( 4 ) 



Let us consider the case, made plausible in chapter ii (p. 
30), that h is large and that therefore the first term in the 
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right-hand side of equation (6) may be omitted. If the mole- 
cules of the substances are rather large and strongly polar 
(chap, iii, p. diffusion coefficients will be rather 

small — say about io“^ cm^ sec“^. If the rates of reactions qi 
are of the order of gm cmr^ sec“^, and the concentrations 
Cl are of the order of 10“^ gm cm“3^ then, as follows from 
equation (3), ki must be of the order of 10“® sec“^ With 
ro^3 X io"~3 cm, the expression rlki/gDur^ 10 > i, and 
equation (5) reduces to 


Cl = 


gDiiCoi 

rlh 


(7) 


The assumptions made here are not essential and are used to 
simplify somewhat the final expressions. 

Putting 


V DiJ)» Din = D ■, 

02 * • • • Cofi — ^0 5 [ 

kn = k , 


( 8 ) 


we have from equations (2), (7), and (8) 


. . . 



(9) 


Let the substance into which the above-mentioned n sub- 
stances are transformed, and which forms the body of the 
cell, break down at a constant rate gm cm"”^ sec"”^ Since 
the rate of formation of the body substance is equal to the 
sum of the rates of consumption of the n component sub- 
stances, the rate of production of the body substance is equal 
to 


Cl ... .Cn ^ ki= kci .... Cn gm cm“3 sec"^ 


(10) 
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Denoting by V the volume of the cell, we have for a con- 
stant density, 5 , of the latter, for the rate of change of its 
mass 5 V the following differential equation: 

= kc,....CnV-q,V. (ii) 

Introducing (9) into (ii) and remembering that V =f7rfg, 
we find 

dro ^Dcok Toqb , . 

o-jT = - — f —■ (12) 

dt rji 2, 

Equation (12), integrated, gives,’ denoting by C a constant of 
integration, 

rl = . ( 13 ) 

rl(jh 

This shows that tends asymptotically to a limiting value 
- jDdok ( 14 ) 

The graph of equation (13) showing the variation of Tq with 
respect to time is given in Figure 19. 

If, besides the reactions discussed here, respiratory and 
other reactions take place in a cell, with the result that above 
a certain size the ceU wiU spontaneously divide in two, as 
discussed in the previous chapter, then two things may hap- 
pen: either the critical radius r^, above which division occurs, 
is larger than h, or it is smaller. In the first case the cell will 
reach asymptotically the size Tq and wiU practically cease 
growing when this size is approached sufficiently closely. In 
the second case it will divide before it approaches fo. As we 
have seen in chapter iii, the stronger the glycolytic coefficient, 
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the stronger the dividing forces and, everything else being 
equal, the smaller the r* . The curve of Figure 19 shows that 
the smaller the cell the more rapid its rate of growth. When 
the cell dhddes at a size r*, the size of the two half-cells is ap- 
proximately o.Sr* (MB, p. 71), and they grow up again to 
r*. Hence, the rate of growth of the ceU is always represented 
by that part of the curve of Figure 19 which lies between the 



Fig. 19 

ordinates r* and o.Sr*- Since the slope of the curve decreases 
with increasing t, therefore, as is readily seen from Figure 19, 
the smaller the rf the greater the average rate of growth be- 
tween o.8r^ and r*. Therefore, the smaller the r* the more 
frequently a cell will divide. Analytically the expression for 
the interval between two successive divisions is obtained by 
considering as the initial value of fo at i = o in equation (13) 
the value o.Szo and determining appropriately the integra- 
tion constant C. Then we obtain the time at which the cell 
divides by making u = r*, substituting for To its expression 
( 13), and solving with respect to t. 
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This gives“ 


I 

At = — log ■ 
2 qb ® 


0.64^6?-*^^ 

gDcok 

gDcok 


(is) 


The quantity At measures the time during which the num- 
ber of cells in a tissue is doubled, and thus At can be accurate- 
ly measured. Its inverse i/At can properly be called the rate 
of growth. 

For small values of r*, when the second terms of the 
numerator and denominator are much smaller than unity, 
expression (15) simplifies considerably. By expanding the 
logarithm, by expressing in terms of | 3 , and by denoting the 
diffusion coefficients for the effective respiratory metabolism 
by Di and De, as in chapter iii, p. 57, we find 


o.iZSk(M:i2Di + be)Y^^ 
D-a \RTM + B)) 


(16) 


Equation (16) shows that, other conditions being con- 
stant, the rate of growth i/At increases as the two-thirds 
power of the glycolytic coefficient /S. According to this theory, 
in its simplified form, which leads to equation (16), a tissue 
grows the faster the greater the glycolytic coefficient j 3 . But 
at the same time the average size of cells composing the tissue 
should be noticeably smaller. Such, however, is not generally 
the case for tumor tissues. This conclusion, however, does not 
follow from the more general equation (15), when kqirplgDcak 
is very close to unity. When it is exactly equal to i, then At 
becomes infinite, and the rate of growth zero, because in this 
case = fo, and the critical size is reached only after an 
infinite time. For qirpk/kgDco > i the interval is imaginary. 




Fig. 20.-See legend on bottom of facing page 
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When r* is still larger, so that the numerator becomes nega- 
tive, then the fraction in the log of equation (15) is positive 
but less than i, so that At is negative. But when is slightly 

smaller than ^^gDcokfkqi,, then a very small variation in rf 
will produce a very large variation of At. Since r* varies as 
and since q is linear in i? (chap, iii), a large increase in 
jS may result in a very small decrease of r* but in a very large 
increase of i/At. In Figure 20 the full line shows the varia- 
tion oil/ At with respect to calculated from equation (15), 
into which expression (26) of chapter iii is introduced for r* . 
The broken line shows the variation of 

Glycolysis may, however, influence the rate of growth in a 
more direct way than that considered above. When jS is large, 
the net effective metabolic flow is directed outward and re- 
sults in a pressure on the cell membrane directed outward. 
This pressure is given by equations (5) of chapter iii. For a 

Fig. 20. — The theoretically expected variation of the rate of growth 
of a tissue and of the cell size with the glycolytic coefficient j 3 . The full 
line represents the rate of growth, according to equation (15), with the 
constants: 35/2^5 = lo'* sec; D =io“* cni=* sec“^; Co = io~3 gm cm'"^; 
qi, = io”s gm cm“3 sec““^; k/k = i. In expressing r J in terms of q and ^ by 
means of equations (26) and (41) of chapter iii, the following values were 
used: ^ = 0.018 X 10“^ gm cm”^ sec; B = 0.738 X 10“^ gm cm'^ 
sec”"^; $M{2Di -|- D^y/RTix = 2 X gm sec~^ The broken line rep- 
resents the variation of the rate of growth with the glycolytic coefficient 
iS, based on equation (23). In addition to the same constants as before, 
the following were used: h = io““ cm sec“^; x = S X 10“^ cm^ dyn~^ 
sec“^; M = 100; Po = 9.51 dyn cm"^. The alternate line represents the 
variation of rj, which is the same in both cases. The graph does not pre- 
tend to any accuracy and is merely used as illustration. A rather high 
value for h has been used in order to justify the use of equations (26) of 
chapter iii, since the latter holds only for large values of h. Otherwise, 
the constants were chosen so as to give the same range of values for 1/ At 
in both cases, and thus illustrate the steeper ascent of the curve cor- 
responding to equation (23). See text, p. 80. 
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rounded-up cell, when ri = fa = ro, the differences Ci — c' 
and Ci — c'i both become equal to 2qral^h, as is seen from 
equations (32) of chapter i, q being defined by equation (41) 
of chapter iii. Thus, the pressure in this case is equal to 


2 RT grp 
9 M h 


(17) 


where h again refers to an average value for respiratory 
metabolites. The pressure, represented by expression (17), 
may be partially compensated by a static osmotic pressure pp, 
owing to different nonmetabolized substances dissolved in 
the ceU and in the surrounding medium — substances for 
which the membrane may be impermeable. The total pres- 
sure is then equal to 


2 RT qfp 
9 M k 


(18) 


Let a reversible reaction take place in the ceU, such as 

Ap -t . . . . An Pi -h Pa + .... Pm , (19) 

where Ai represents some amino acids and Pi represents 
proteins. If this reaction requires a special enzyme, which is 
present only inside the cell, then it will not take place out- 
side the ceU. If the membrane is permeable to Ai but not to 
Pi, then the concentrations of Ai will be the same inside the 
cell and outside; and if the cell size remains constant, the 
amounts of Pi wiE not change, everything remaining static. 
However, let the ceU size increase very slowly as a result of 
some causes. Then the concentrations of A i and Pi wiE de- 
crease. But, because of the supply from outside, the concen- 
tration of Ai wEl be practicaEy immediately re-established. 
Thus the concentration of A i wEl remain the same and that of 
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Pi will decrease. As a result of this, some Pi will be formed ac- 
cording to equation (19). The cause of the increase may be the 
pressure given by expression (18), which tends to expand the 
cell. According to this picture, a cell does not grow in size be- 
cause of increase in mass, but rather increases in mass be- 
cause of growth in size. 

In the first approximation the rate of increase in size will 
be proportional to the expression (18) and will be given by 


drp 

dt 


X 


2 ^ grp 
9 M h 


(20) 


where x is principally determined by the resistance of the 
membrane to the flow of water, which, because of its incom- 
pressibility, must flow into the cell whenever the latter ex- 
pands. Equation (20), integrated, gives 

rp= 9. . (21) 

2 Klq 

For jS just above the critical value /5o = —BjA^ when 
q + B is positive and very small, such a cell might 

grow to tremendous dimensions, though the rate of growth 
would be extremely slow. For normal jS it will grow to a 
normal size r* and then divide. Again, as before, let us put 
t = o, when Tq = o.8r^. We then obtain^ 


_ gMhpo 

2RTq 



. (22) 


The interval At between two successive divisions is determined 
from 


gMhpo 

2RTq 



gMhp^ 

2 RTqj 


J2RT-x2iQWfi)At — r* 

& / 0 > 
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or, solving it for A/, 

_ gMh gMhpa — 2 RTqrt 
iRT xq ® gMhpo — i. 6 RTqr^ 


(23) 


For very small values of po, equation (23) reduces to 


At 


gMh 

2RTxq 


log 1.25 


Mh 
RTxq ' 


(24) 


Again expressing the critical size r* in terms of the glyco- 
lytic coefficient / 3 , by means of equations (26) and (41) of 
chapter iii, we obtain from equation (24) the interval At in 
terms of / 3 . Such a relation, obtained from equation (24), is 
shown in Figure 20 by the broken line. 

Unfortunately, no quantitative data are available for com- 
parison with the calculated curves. We see from Figure 20 
that a variation of 18 from o.i to 4 results in a tremendous in- 
crease of the rate of growth. It is known empirically^ that 
values of ^ about i correspond to benign tumors, while values 
such as 3 or 4 correspond to de&iitely malignant, very much 
more rapidly growing tumors. The importance of exact 
quantitative e.xperimental studies of the rate of growth of 
different tumors, together with simultaneous determinations 
of the glycol3rtic coefficients of the same tumors, becomes ap- 
parent from these considerations. By comparing the results 
of such experiments with the two curves of Figure 20, either 
we may be able to decide in favor of the one or of the other 
of the equations here derived, or else we may get definite 
indications for the necessity of considering more complicated 
mechanisms involving the effects of glycolysis on the rate of 
growth.* 

A careful, well-controUed quantitative statistical study of 
cell sizes in tumors, as compared to normal tissues, is also 
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urgently needed. From Figure 20 we see that the change in 
cell size of only a small percentage may correspond to a tre- 
mendous increase in the rate of growth. It is usually ac- 
cepted that tumor cells do not show any appreciable differ- 
ences in size from normal cells. But this is based on rather 
crude general observations. A difference of a small per- 
centage will escape notice unless a thorough statistical study 
of size distributions is made and mean values and standard 
deviations computed and compared. 
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CHAPTER V 


CELLULAR FORMS AND MOVEMENTS 

In chapter iii we have studied the deformation of the cell 
under the influence of diffusion forces for the case where the 
rate of deformation is very slow. In this case we may con- 
sider the diSusion field at any moment as quasi-stationary. 
Since in some cases the ceU division lasts only lo minutes or 
less, it is important to discard the restriction of slow rate of 
deformation and to attack the problem in its more general 
aspects. 

In this case the variation of the concentration c in the cell 
is given by equation (19) of chapter i. However, the quantity 
A, which is a function of the ceU shape, now itself is variable, 
its variation bemg determined by the diffusion forces accord- 
ing to the discussion of chapter iii. These diffusion forces are 
themselves determined by the values of c, Ct, c^, etc. 

For the sake of simplicity, we shaU confine ourselves to the 
case of a very large permeability, h, though the more general 
case is treated in the same manner. We also consider, again 
for sake of definiteness, the production of a substance; that is, 
we take 5 > 0. 

Making h= <x> and 5 = in equation (18) of chapter i, 
introducing the expression for A into equation (19), and 
introducing a new variable 

c = c - Co, (i) 

we may write equation (19) of chapter i in the form 
— = — 3 -Pi-Pe rl{De + 2 Di) 4 - ijnDe -h 

dt ® (riD,+ 2riDi){De-h 2Di) ^ 

82 


(2) 
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The rate of elongation of the cell is given by equation (16) 
of chapter iii, which in that form holds for the general non- 
stationary case. Putting in that equation h = r, and adding 
to the right-hand side the surface-tension term (chap, iii, 
Eq. [24]), we obtain 


r-i it sMtj 

~ f'l) _ 7 rt — r, 

^ i2r,Di+r,D,){2Di + De)^ 2r, ' 


We shall now discuss the simultaneous differential equations 
(2) and (3). These equations are nonlinear, and their exact 
solution presents very great difficulties. We shaU, however, 
study the general properties of their solutions by means of a 
semigraphical method which is rather simple and useful. 

The requirement 


dc ^ 

di-°’ 


(4) 


introduced into (2), gives 

= r^rlq (r^De+ 2r2Di){De+ ^Dj) 

^ ^ sDiDe 2rj (riDe + 2rJDi) + rliJDe + 2Di) * 


The equality sign in (5) gives the equation of a curve for 
which dc/dt = o. Since 


^ iZl. 


V = iirrl ; u 


47rre 

the equation (5) may be written 
Vq(De+2Dd 


4'ir ] 


1/3 


c = 


47rDJ)e 


Deft + 2D, 


X 



iDefl + 4 ^ 


\ 


3F(D. -f 2Di) 

4 wri 


( 6 ) 


(7) 
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For verj' large values of r-, we have 


Vq(J),-\-2Di) I 
STrDiDe ' 


( 8 ) 


For very small values of fi we have 


( 9 ) 


From (8) and (9) we see that c = o f or ri = o and for = 00 . 
The value of c is everywhere positive. From physical con- 



siderations it is cle ar that c should have a maximum for 
f Q = '^3F/47r, that is, when the cell is spherical. For 
De = <x> this is actually the case, as can easily be seen. 



CELLULAR FORMS AND MOVEMENTS 85 


For finite De (7) gives a maximum at a point slightly different 
from Tq. This is due to the approximation introduced in 
chapter in by making S = fz, which is not quite exact* 

Let us consider Tj and c as rectangular coordinates in a 
plane (Fig. 21). Equation (7) then represents a line, with a 
maximum at approximately ri = n (Fig. 21, full line). For 
all points of the plane which lie below that line, dcjdt > o. 
In other words, if we choose any pair of values and c, that 
is, a given elongation and a given concentration, such that the 
representative point is below the line, then for such a pair of 
values the concentration c ^ c — Co will increase. The re- 
verse holds for all points above the line. Here we have 
dc/dt < o. 

Because of equation (3) the requirement 



VII 

0 

(10) 

leads to 

_ + Z).) / 2 Di De\ 

zRTtxDJ), 

(ii) 

for 

rt> ft, 


and to 

_ M^{ 2 Di + A) / 2 Di I>e\ 
SRTiiDiD, \r,^rj 

(12) 

for 




^1 < f 2 . 


Because of (6), equations (ii) and (12) may be written thus: 
. My{2Di-\- [2Di , De r ^ r \ 
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and 

< -^T(2.Pi + A) 2Z?i , ^ /— 

sRTfjJ)iDe r, 


for < r2. (14) 


For n = ^2 equation (3) gives always drj:/dt = o. The 
equality sign in (13) and (14) gives us c as a function of r^. 
This function is infinite for = o and varies for small values 
of like i/fi; while for larger values of Tx it varies as Vfj. 
(Fig. 21, broken line). It has a minimum for 




Fig. 22 


The two lines (7) and (13) either may not intersect at all 
(Fig, 21) or may intersect as in Figure 22 or in Figure 23. 
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The situation shown in Figure 21 will occur when, for a 
given To and other constants, the value of q is sufl&ciently 
small, resulting in too low a concentration c ^ c — Cq Sit the 
maximum, which corresponds to a sphere. In this case the 
two lines represented by equations (7) and (13), together with 
the line AB, which corresponds to fa = fo, divide the quad- 



rant of the fj^c plane into six regions. In region I, dc/di < o, 
while > o. If we consider a cell, characterized by a 

given elongation fj and a given concentration c = c — Co, and 
if to these particular values of fj and c the corresponding 
point lies in region I, then fi will increase while c will decrease. 
We may say that at a certain moment any configuration point 
in the r^c plane moves in the direction indicated by the ar- 
rows. In region II we have dc/dt < o and dr^/dt < o; the 
direction of the movement of a configuration point is also 
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indicated by the arrow. In region III, > o and 

dcja < o; while in region IV, drrjdt < o, icidi < o. A 
similar survey of the other regions leads us to a distribution 
of pathways along which the configuration point moves ac- 
cording to equations (2) and (3), as shown in Figure 21. For 
any point on we have dr^ldt = o. Therefore, above the 
heavy line any such point will move downward; below it, up- 
ward. All paths converge to the point 0 , which corresponds 
to = ra == ro, that is, to a sphere. This is quite clear 
physically, since for very small values of g the sphere is a 
stable shape of equilibrium. 

A similar analysis of the case represented in Figure 22 
leads to a different distribution of the paths. The point 0 is 
now a point of unstable equilibrium. There are now two 
points, Oi and O2, of stable equilibrium, toward which the 
paths converge. The point Oi corresponds to < ro, or 
n < fa — in other words, to a flattened shape. The point O2 
corresponds to fi > fo or fj > u — that is to say, to an elon- 
gated shape. This happens when, owing to a sufficiently 
large g, the sphere becomes imstable for a given fo, that is, 
when fo becomes greater than the critical radius f * given by 
equation (26) of chapter iii. In this case the metabolic forces 
begin to elongate the cell. Either this elongation may pro- 
ceed in one direction only, leading to oblong shapes with an 
eventual division in two; or it may proceed in two transverse 
directions, leading to flattened shapes, like, for instance, red 
blood corpuscles. 

In the case represented in Figure 23 the spheroidal shape, 
corresponding to f^ = fo, is in a relatively stable equilibrium. 
It requires a finite disturbance to make the cell elongate to 
the point O2. 

Around the point O2 the path of the configuration point 
may be a closed line. This means that the cell may periodical- 
ly contract and expand around a configuration of equi- 
librium.^ Let us now keep aU other parameters in equation 
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(7) constant but increase V or, wbat amounts to the same, 
increase fo = (3F/47r)^^3. Xhe point 0 in Figure 23, cor- 
responding to the maximum of the hea\y line, will shift to the 
right and will also move upward, since with increasing V the 
maximum concentration c increases. Finally, 0 will fall on the 
broken line, so that 0 and Oi will coincide. In other words, 
the sphere will become a configuration of imstable equi- 
librium, and any infinitesimal elongation will bring the cell 
into configuration O2. The value of n for which this happens 
is nothing but the value given by equation (26) of chapter 
iii. The foregoing reasoning proves the assertion made in 
chapter iii, page 50, that, in some cases when the spheroid 
shape becomes unstable, the deformation proceeds at once to 
a fimite value. 

A similar situation will occur for the case in which the 
maximum of the heavy line is to the right of the mini-mum of 
the broken line. In this case, however, as soon as instability 
for spheroidal shapes will set in, the cell will flatten out by a 
finite amoimt and assume the configuration corresponding to 
the point Oi of Figure 22, 

For very slow deformations, such as studied in chapter iii, 
the configuration point of the cell in the (fi, c) plane moves 
along the heavy line, which gives at each point the value of 
c = c — Co, corresponding to a quasi-stationary value for a 
given elongation, fi. In general the configuration point moves 
along one of the pathways indicated by the arrows. 

Considering, for instance, the case where Di De, so 
that De may be neglected, as compared with Di, when oc- 
curring together as a sum or a difference, and considering, 
moreover, the case where the point O2 is very much to the 
right of 0 , we may use,^ instead of equation (7) the approxi- 
mate equation (8). The latter now simplifies to 


Vq I 

fi ’ 


(16) 
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and equation (13) now simplifies to 


c 


2My 

2,RTiJ.De 


2Di ■ 

f, 


A . /- 


Vr 


(17) 


By equating the right-hand sides of equations (16) and (17) 
we find the abscissa of the point Oa. Making use also of equa- 
tion (26) of chapter iii, and putting 


V* = 



we find approximately^ 

fi — To 


{fr 


(18) 


(19) 


In a similar way we find the value Cz of c corresponding to 
the point Oz- 

If the condition discussed in chapter iii, page 53, is satis- 
fied, then the cell not only will elongate but also will divide. 
Otherwise it remains in an elongated nonspherical shape 
indefinitely. 

Since by means of equation (26) of chapter iii, rf is ex- 
pressed in terms of the effective metabolic rate g, equation 
(19) gives us the elongation of a cell in terms of the effective 
metabolic rate. 

We have pointed out before {MB, chap, ix) that, because 
of the presence of metabolic forces, we should expect the 
possibility of existence of nonspherical shapes of equilibrium 
in liquid cells. We also pointed out that, inasmuch as the 
existence of such shapes is tied up to the presence of the 
metabolic forces, we should expect that a cessation of metabo- 
lism — in other words, the death of the cell — would result in a 
rounding-up of the cells unless there are some rigid structures 
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present within the cell. This agrees with the general observa- 
tions that a number of cells which have an oblong shape dur- 



Fig. 24 , — The theoretically expected variation of the length- to- width 
ratio fz/fz of a unicellular organism with varying rate of metabolism q, 
w^hen the latter is influenced by external agencies, for the following val- 
ues of constants: Dq = 10“^ cm^ sec"^; == lo"^ cm=* sec~^; 7 = i erg 
cm“^; = o.i; 2RT /M = 5 X lo^ erg gm'^. The case corresponds to that 

of Fig. 22- The discontinuity is due to the fact that, for that particular 
value of q, the configuration Ot on Fig. 22 becomes unstable and the cell 
assumes the flattened shape, corresponding to 0 i of Fig. 22. If q is again 
increased, the cell flattens out further (broken line). We have, thus, a 
hysteresis effect. Such an effect is not present for the choice of constants, 
which correspond to Fig. 23. 

ing life round up upon death, as discussed in MB, chapter Lx. 
But, as has been pointed out, we should expect not only a 
complete roxmding-up of the cell upon death but also a 




92 


MATHEMATICAL BIOLOGY 


gradual rounding-up during life if the cellular metabolism is 
impaired by some drugs (MB, p. 103). Although this seems 
to be the case, no definite quantitative observations are avail- 
able, however. The discussion of this chapter leads us to a 
definite prediction concerning the dependence of the elonga- 
tion r'/ of a cell as a function of the effective metabolic rate q. 
The use of equation (19) in this case is not justified, for that 
equation holds only for a rather narrow range of values of 
r('. We can, however, easily construct rC as a function of q 
graphically, by plotting the lines given by equations (7) and 
(13) for different values of q, by using a given set of values for 
the other constants, and by determining graphically the 
abscissa of the point O2 for every value of q. The result of 
such a determination is shown in Figure 24. 

Since the effective q is itself a function of the other com- 
ponents of metabolism, we may thus study the effects of 
those components on the shape of the ceU. The urgent 
necessity of stud3dng the effects of various drugs and other 
environmental conditions affecting metabolism on the shape 
of protozoa becomes evident. 

As we have said, around the point O2 the cell may exhibit 
periodical contractions and expansions. This case has been 
studied approximately analytically,* and the bearing of it on 
a possible biophysical theory of such phenomena as con- 
tractility of tissues and amoeboid movements has been dis- 
cussed. 
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PROTOPLASMIC STREAMING 

In all cases hitherto discussed the diffusion forces were 
directed either outward or inward, tending either to disrupt 
or to compress the cell. It is conceivable that in some cases 
a sufficiently as5Tnmetrical distribution of these forces inside 
the cell will produce not merely elongations or contractions 
of the cell but continuous displacements of different parts of 
the ceU with respect to each other, and thus produce internal 
movements in the ceU. Such internal movements are, as a 
matter of fact, frequently observed in actual cells. They form 
a wide class of phenomena, known as “protoplasmic stream- 
ing.”^ While in plant cells these protoplasmic streamings oc- 
cur rather continuously during the normal life of the ceU, in 
some animal cells they appear in different forms during the 
process of ceU division. 

As mentioned above, some sufficiently as3mmetric dis- 
tribution of the forces must be present in order to induce such 
movements. Such an asymmetric distribution is, however, 
quite likely to occur whenever the cell exhibits gross inhomo- 
geneities of its structure, which is more a rule than an excep- 
tion. Thus, for instance, in most cells the nucleus does not 
occupy a central position but lies eccentrically. Quite regard- 
less of any detailed assumptions we make about the chemical 
and physiological role of the nucleus, if we consider that it is 
a seat of some important reactions, we see that an asym- 
metry of location of the nucleus will result in an as)mmetry 
of the reaction rates at various parts of the cell. In other 
words, instead of having an average constant rate of reaction 
q per unit volume, we now have the situation in which the rate 
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of reaction for some metabolites will be much larger near one 
end of the cell than near the other. If this is the case, then 
the concentrations of the corresponding metabolite will not 
be symmetric with respect to the equator of the cell. Con- 
sidering again, for sake of definiteness, the production of one 
substance, we can readily see that, if in the region AB (Fig. 
25) the rate of production is higher than in the region A'B\ 
then the concentration near the end AB will also be greater 
than the concentration near the end A'B\ Instead of having 



Fig. 25 


a maximum value near the center of the cell and dropping 
down toward loth ends, the concentration now will have a 
maximum somewhere near that indicated in Figure 25 by the 
shaded area and wiU drop toward the end A^B'. A similar 
thing will happen if the rate of reaction q is ever3rwhere the 
same but the permeability of the cell membrane at AJ5 is 
much less than the permeability at A'B\ Such a difference in 
permeability may also be caused by an asymmetric location of 
the nucleus if the latter produces substances which affect the 
permeability. There will therefore be a force directed from 
left to right, as indicated by arrows. The force represented 
by the upper arrow tends to set the interior of the cell into 
clockwise motion; the force represented by the lower arrow 



PROTOPLASMIC STREAMING 


95 


tends to set it into counterclockwise motion. In order that 
the interior of the cell may be maintained in a steady motion, 
either clockwise or counterclockwise, continuous work must 
be done by the driving forces, because, owing to the \dscosity 
of the cytoplasm, any continuous motion is connected with a 
dissipation of energy. This implies that the work done by the 
forces in carrying a particle of the cell interior around a closed 
path would be finite. As we have seen (chap, hi), the diffusion 
forces may be considered as being derived from a potential 
which is proportional to the concentration. In that case, how- 
ever, the integral of the force along any closed circuit is al- 
ways zero. The circumstance that now we may have a dif- 
fusion in a moving medium does not alter this situation. The 
motion of the medium will affect the distribution of the con- 
centrations, making the gradients different, but the force still 
remains derived from a potential. Thus, at first sight it 
would appear that the diffusion forces cannot produce any 
permanent circulatory motion of the cytoplasm. A closer 
examination, however,^ shows that this is not necessarily the 
case. The foregoing argument holds if we consider that the 
whole of the diffusion force is sustained by the liquid medium. 
If, however, there are some immobile rigid structures present, 
then part of the force will be sustained by those, and that 
part will not contribute to the motion of the liquid. Either 
these solid structures may be constituted by a jelly-like 
meshwork, through which the liquid may be “squeezed 
through,’’ or the rigid walls of the cell may sustain part of the 
force, which will be transmitted to them by the viscosity of 
the liquid. 

Suppose, now, that the physical constitution of the rigid 
framework or of the walls is different in the top and in the 
bottom part of Figure 25, so that in the two regions different 
fractions of the total force are sustained by the liquid. Then, 
while the two total forces, represented by the arrows, will be 
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equal, their parts sustained by the liquid will not be equal. 
The total integral of that part of the force along a closed cir- 
cuit will not be zero, and streaming will result. 

Gale Young has developed a theory of protoplasmic stream- 
ing along those lines. ^ Referring to Young’s original paper for 
a more exact treatment, we shall give here merely an approxi- 
mate outline of the mathematical argument, which leads to 
somewhat similar results. 

We may estimate the order of magnitude of the upper 
limit of the velocity of streaming in the following manner. 
Suppose, for definiteness, that the asymmetry of the ceU, 
represented schematically by Figure 25,. is such that a greater 
portion of the force represented by the upper arrow is trans- 
mitted to the liquid, so that that force prevails and the liquid 
is set in clockwise circulation. The strongest driving force 
would be obtained if the force represented by the upper arrow 
would be completely transmitted to the liquid, while the force 
represented by the lower arrow would be completely sus- 
tained by the solid structure. In that case the total driving 
force would be represented by the upper arrow. According 
to equation (3) of chapter iii, this force is of the order of 
magnitude of {RT/M) X grad c per unit volume. Denoting 
the concentrations at . 4 J 5 and at by Cj and C2, respec- 
tively, we find that the average value of grad c is of the order 
of (^i — cP)/ 2^1. Hence, the driving force per unit volume is of 
the order of 

RT{cj — C2) . . 

2Mn 

Considering, for simplicity, the case where the asymmetry of 
concentration at the two ends is due to a very low permeabil- 
ity at AB and a very high one at A'B', we can put C2 = Co, 
where Cq is the external concentration. 

We shall now express Cj — 62 = — Co in terms of the re- 
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action rate q. Besides the transport by diffusion, a transport 
of substance by convection takes place. If v denotes the aver- 
age velocity of streaming, and if c is the average concen- 
tration in the cell, then the convection transport is equal to 
VC gr cm""' sec“^ (MB, p. 13). Remembering that the velocity 
has a different direction in the upper and in the lower part of 
the cell, and that the cross-section of each part is approxi- 
mately (tt/ 2)rl, we find that the net transport from left to 
right is 



VC -{■ Di 




2 

2^^ 


= Di TTfl 

2rr 


(2) 


This net transport is of the order of magnitude of the total 
production rate ^irrjrlq. Strictly speaking, we must also take 
into accoimt the loss of substance through the ^^sides” AA^ 
and To neglect this introduces a rather large error but 
does not affect the order of magnitude, which is all that we 
are now concerned with. Equating expression (2) to the 
total rate of production and solving for Cj, — Co, we find 


('l Cq 


§ 2!l 

sDr 


(3) 


Combining this with expression (i), we find the driving force 
per unit volume to be of the order of 


4RT 

iMDi ■ 


(4) 


The total force is equal to expression (4) times approximately 
half the volume of the cell, that is, times When there 

is a steady streaming, then the total force must be equal to 
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the viscous resistance. Denoting by rj the viscosity, and re- 
membering that the average “gradient of velocity’' is 2z;/r2 
(because over the distance 2^2 the velocity changes from +7; to 
—v), we find the resistance per unit area to be 22/77/^2. The 
total resistance is obtained by multiplying this by the total 
surface, which is of the order of 4rxr2, Now, by equating the 
total resistance to the total force, we obtain a linear equation 
for the determination of v. Solving it, we find 


TrRTqTxfl 

gMDiri 


(s) 


Taking q ^ io“® gm cm""^ sec”^, rj ~ ~ cm, Di ~ 

lo"'^ sec~b a,nd 17 i, that is, a value a hundred times 
that of water, we find 2; ~ i cm sec^h The actual velocities^ 
are of the order of 10“^ — lo”^ cm sec“h This is not surpris- 
ing, for we have given here the upper limit of the velocity, 
under the assumption that the difference of the two oppositely 
acting forces (Fig. 2 5) is simply equal to one of them, the other 
force being zero. The fact that the actual velocities are so 
much smaller than the upper limit shows that a very small 
asymmetry, resulting in a very small difference in the amounts 
of the forces sustained by the solid structure, will be suf- 
ficient to produce streaming. 

Another interesting consequence of equation (5) is the fol- 
lowing: The velocity v depends both explicitly and implicitly 
on the temperature. The temperature-depending quantities 
on the right-hand side of equation (5) are g, Di, and rj. Inas- 
much as Di is approximately^ inversely proportional to tj, 
being, besides, proportional to T, the product Dir) is ap- 
proximately proportional to T. Therefore, approximately, 
the velocity v varies with temperature, as does q. If in an 
actual cell we would have to deal only with a single reaction, 
we would expect the rate of that reaction to increase ex- 
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ponentiaUy with temperature, according to the Arrhenius 
equation, that is, to vary like being a constant). As 

we have seen in chapter iii, g is actually a composite quantity, 
and therefore the relation g can hold only approjd- 



Fig. 26. — Logari thm s of velocity of protoplasmic streaming plotted 
against the inverse of absolute temperature. Experimental data from 
C. Oppenheimer and L. Pincussen (1927), Tabulae Uologicae, 4» 465 
(Berlin: W. Junk). The relationship is approximately linear, as required 
by the theory. 

mately. To the extent of this approximation we should, how- 
ever, have V ~ or 


log D = const — ~ . (6) 

In other words, if we plot on semilogarithmic paper the loga- 
rithm of the velocity of streaming against the inverse tempera- 
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CHAPTER VII 

EXCITATION IN PERIPHERAL NERVES 

In this chapter we shall discuss another rather general 
property of living cells, namely, their irritability. VTiile, to 
some degree, irritability is found in almost all cells, the nerve 
cells exhibit it to a particularly high degree. A great many 
different physical and chemical factors, when acting on a 
limited region of an irritable cell, produce at that region a 
peculiar physicochemical state characterized by a local change 
in the metabolic rates, an increased general permeability, and 
an electronegativity relative to the nonaffected regions of the 
cell. This physicochemical state has the property of spreading 
to adjacent regions which are not directly affected by the 
external local physical or chemical factors. It thus propagates 
all over the cell. In the case of a nerve fiber, which is a highly 
elongated part of a cell body, this physicochemical change, 
which we call “excitation,” propagates along the fiber, travel- 
ing over appreciable distances. 

Another peculiarity of the excitatory state is that it usually 
develops suddenly, after the intensity of the external physical 
or chemical factors which induce it exceeds a rather definite 
threshold value. 

The problem which confronts the mathematical bio- 
physicist in this case is to derive the existence and properties 
of excitation from other basic biophysical properties of the 
cell, such as metabolic activity. In view of the great general- 
ity of excitatory phenomena, such an attempt is highly justi- 
fi.ed. Unfortunately, until now only a few general indications 
are available as to how phenomena of excitation may be re- 
duced to other more general ones. The actual solution of the 
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problem is still wanting. This, however, should not detain us 
from a purely formal mathematical study of the properties of 
excitation, though we may not know what the true nature of 
excitation is. In the early development of physics, phenomena 
of heat, optics, and electromagnetism were treated as sepa- 
rate, rather unconnected, branches. Now we know that 
optical phenomena are merely a special case of electro- 
magnetism. But even long before the physicists knew that, 
they were quite able to develop the mathematical theory of 
optical phenomena, though the true nature of light was un- 
known to them. And, when this true nature was revealed by 
the work of Maxwell, all the studies on mathematical optics 
were automatically incorporated into the general electro- 
magnetic theory of light. 

We shall follow the same procedure with respect to the 
excitation phenomena, pending their hoped-for reduction to 
more basic phenomena of cell biophysics. 

Under physiological conditions different t^ypes of nerves 
require for their excitation different physical or chemical fac- 
tors. Thus, the taste nerve fibers are excited by chemical sub- 
stances. The acoustic nerve fibers and those of tactile sense 
respond to mechanical pressure changes, etc. 

Under artificial conditions of the experiment we find, how- 
ever, that there is one, so to say, tmiversal factor for excita- 
tion, namely, the electric current. The protoplasm of the 
nerve cell through which any exciting current passes is an 
electrolyte. Therefore, the most obvious and immediate ef- 
fect that the passage of a current may have upon the proto- 
plasm is a displacement of various ions in one or another 
direction, depending on the sign of their electric charge. Such 
a displacement will result in local variations of concentrations 
of different ions. On the other hand, the classical work of 
J. Loeb and others has shown that the excitation state may 
be induced in many cells by treating them with solutions of 
ions, and that a certain threshold concentration of ions in the 
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solution must be exceeded in order to induce excitation. Con- 
siderations similar to these led W. Nernst/ in 1899, to assume 
that the excitation by electric currents is due to such changes 
in the concentrations of ions which, when they exceed a 
threshold value, induce excitation. 

Suppose we let an electric current of some kind pass 
through a nerve. Then, if we can set up the differential equa- 
tions which describe the motion of ions and the change of 
their concentrations under the influence of such a current, 
and if we can integrate these equations, we should be able to 
predict how long such a current would have to pass through 
the nerve in order to excite it. 

The differential equations involved in such considerations 
are partial differential equations of a form whose integration 
sometimes presents very great difficulties. Setting up and 
solving a three-dimensional problem with approximately 
axial symmetry, as we have it in a nerve fiber, presents almost 
insuperable difficulties. Therefore, beginning with Nernst, 
different investigators'* have resorted to exact solutions of 
highly oversimplified cases which had hardly any properties 
in common with the real cases. Yet some general results, 
agreeing with observations, were obtained. The situation was, 
to some extent, similar to the exact treatment involving the 
diffusion problem for a perfectly spherical ceU, a problem dis- 
cussed in mathematical biophysics, except that the unsatis- 
factory aspects of the procedure were even more apparent in 
the case of nerve excitation. 

-We have seen how, in the case of mathematical biophysics 
of the cell, the introduction of a rather drastic but powerful 
approximation method reduces problems of partial differential 
equations to those of ordinary ones, and problems of ordinary 
differential equations to problems of algebraic equations. A 
procedure similar to that seems indicated, therefore, in the 
case of nerve excitation. 

If an electric current is flowing through the nerve fiber, then 



104 


MATHEMATICAL BIOLOGY 


the positive ions are driven to the cathode; the negative ones, 
to the anode. Let us first confine ourselves to the former. 
Although the actual concentration of the ions will vary from 
point to point in a complex way, the average concentration, c, 
will increase at the cathode approximately at a rate propor- 
tional to the flowing current, I. The coefScient of proportion- 
ality, Kj may be roughly evaluated from the average initial 
concentration the mobility of ions, and other electro- 
chemical properties of the protoplasm and cell membrane. 
This particular problem of determining K in terms of other 
constants will not be discussed here. 

Any excess of concentration will tend to decrease either 
because of a diffusion-like process oj because of a chemical 
reaction which involves the destruction or removal of the 
ions. The rate of decrease will be roughly proportional to the 
amount of excess, c — Cq, Denoting, here, the coefficient of 
proportionality by k, we thus arrive at the differential equa- 
tion 

f^ = KI-k{c-c.), (i) 

in which I is some function of the time t. 

If, as said above, we assume that excitation occurs when 
the concentration c exceeds a critical value then, by intro- 
ducing the new variables 

€ = r — Co; h — — Cq j 

we have 

with excitation occurring for 


( 2 ) 

( 3 ) 


€ ^ A. 


( 4 ) 
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Equations (3) and (4) form the basis of H. A. Blair’s theory 
of excitation, 3 which was proposed several years ago, although 
the argument given here and leading to its approach was not 
used. 

In the absence of any current /, the quantity e varies ac- 
cording to 

€ = 606-*' , (5) 

where €0 is the value which e has at the moment ^ = o. This 
shows that, when for a sufficiently long time the nerve has 
not been subject to any current, € = o and c == CojSls should 
be the case, physically. If, under these conditions, we now 
establish a sudden current of constant intensity I and count 
the time from the moment of the establishment of the current, 
then the quantity € varies according to 

e = ( 6 ) 

It tends asymptotically to the limiting value Kl/k; and the 
initial slope of the e, t line is equal to KL Both the limiting 
value and the initial slope increase with increasing current 
I {MBj p. 151). If, in order to produce excitation, the value 
of e must exceed the threshold h, then Kl/k must exceed h; 
or, in other words, the current I must exceed a critical value 
R = kh/K. This value R is technically known as the ^^rheo- 
base.” Once the rheobase is exceeded, the quantity € exceeds 
the threshold h after a time t, this time-interval decreasing 
as I increases. The relation between the intensity of the 
exciting constant current and the time which this current 
must take in order to produce excitation is obtained by 
equating the right-hand side of equation (6) to h and by solv- 
ing this with respect to L As has been shown by H. A. Blair, 
the relation thus obtained is in fair agreement with experi- 
mental data .3 
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One important conclusion is to be drawn from the compari- 
son of the theoretical results Tvith the experiments: the de- 
crease of e — c — Co cannot be due to any diffusion-like 
process. The order of magnitude of the coefficient k is about 
10* — io 3 sec~*. In other words, it must take about one- 
hundredth of a second or less to re-establish practically the 
initial concentration. As we have seen in chapter i, with 
plausible values for diffusion coefficients the time it takes for 
a diffusion process to reach a stationary state is of the order 
of a few seconds at least (cf. also MB, p. 12). Notice that 
equation (i) is essentially the same as equation (20) of chap- 
ter i, except that, instead of the accumulation of substance, 
owing to the rate of production q, we have an accumulation 
due to the current I. Hence, we should have k= i/A. But 
with any plausible values for the constants involved in A, we 
cannot obtain such high values of h, and that discrepancy is 
found by comparing Blair’s theory with experiments. And 
this holds also for the generalizations of Blair’s theory, to be 
discussed later in this chapter. The only conclusion is that, 
while the ions are accumulated by the current, their concen- 
tration decreases, owing to some chemical reaction whose rate is 
proportional to c — Co. The initial finite concentration Co may 
be due to the fact that there are two reactions going in oppo- 
site directions— one producing, the other destroying, the par- 
ticular type of ions. These two reactions would be in equi- 
librium for the particular value Co of the concentration. 

Instead of a suddenly estabhshed constant current, we may 
consider any other kind of current, such as an exponentially 
varjdng current or an alternating current. This amounts to 
making the quantity I in equation (3) a definite function of 
time, integrating the equation, and again setting e = A. In this 
way a number of other results in more or less good agreement 
with experiments were obtained by H. A. Blair^ (cf . also MB, 
chap. xvi). 
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There are, however, several difficulties of a fundamental 
nature which are not taken care of by H. A. Blair’s theory. 
One of the most important is the behavior of a nerve at the 
anode of an electric circuit. We have hitherto confined our- 
selves to the cathode, to which positive ions are driven by the 
current. When a nerve is made to be a part of an electric cir- 
cuit, and a current is suddenly established, excitation always 
occurs at the cathode. This indicates that the excitatory ions 
are positively charged. But the establishment of a constant 
current of indefinite duration usually results in a single excita- 
tion impulse, which may be observed either by a short twitch 
produced in a nerve-muscle preparation or by oscillographic 
recording of the propagation of the action potential ac- 
companying the excitatory process. After the first impulse the 
nerve is not excited, though a constant current is passing 
through it. But if such a constant current is suddenly inter- 
rupted, then another excitatory impulse is observed, this time 
originating at the anode of the circuit. This phenomenon 
offers very great difficulties to any theory based on the as- 
sumption that excitation is induced by an excess of concen- 
trations of some ions. If positive ions are excitatory, as seems 
evident, then at the anode their concentration decreases, and 
a sudden break of current results merely in the restoration of 
the initial concentration. One might think of some inertial 
effects which would make the concentration of the ions at the 
anode ^^swing” over the initial normal value. But the slowness 
of the diffusion processes makes the introduction of such a 
hypothesis physically quite untenable. Thus, until recently, 
the “break excitation” at the anode presented a stumbling- 
block to all theories. 

In 1933 the author^ suggested a natural generalization of 
H. A. Blair’s theory, which solves the difficulty in a simple 
way. It is known that the protoplasm of any cell contains not 
one, but many, kinds of ions. Moreover, not all positive ions 
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produce excitation. Some of them inhibit the excitatory ef- 
fects of other ions.^ This suggests the study of a case where 
we have not one, but at least two, t 3 ^es of ions — one excita- 
tory, the other inhibitory — both positively charged and there- 
fore driven to the cathode. Denoting the concentration of the 
former by 6 and of the latter by we assume that for each of 
them an equation of the form of the equation (i) holds'**' except 
that the constants involved are different for the two types of 


ions. Thus, we now have two equations: 


^^=KI-k(e- €0) 

(7) 

and 


MI - m{j - j\) , 

(8) 


where again Z", M, k, and m are constants and €o and jo are 
the initial normal concentrations of e and j. In line with our 
knowledge of the antagonism of ionic action, ^ we assume that 
excitation occurs when the ratio e/j exceeds a critical value h. 
As far as single nerve fibers are concerned, we may, without 
loss of generahty, put A = i.f We thus have, as a condition 
of excitation, 

« = y • (9) 

Since in the absence of any current the nerve is not excited, 
we also have 

€o<io. (10) 

* Note that the symbol e stands here for a quantity different from 
that in equation (2). 

t Cf. the following chapter, p. 132, concerning the cases when this 
restriction is not applicable. 



EXCITATION IN PERIPHER.\L NERVES 


109 


Furthermore {MB^ chap, rvdi), we make the following as- 
sumption about the constants iT, M, k, and m: 

k'^m] K^M) — . (ii) 

km ^ 

Figure 27 represents graphically what happens at the cath- 
ode and at the anode when a constant current, /, is suddenly 




e 
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Fig. 27 


closed or opened. Owing to inequalities (ii), the initial slope 
of the €-line is greater than that of the 7-line. The quantity € 
also approaches its limiting value more rapidly. As a result, € 
exceedsy for a short time, during which excitation takes place 
at the cathode. Upon opening the current, the quantity j 
at the cathode decreases less rapidly than the quantity e, be- 
cause olk> m. Thus, in this case e always remains below 7’, 
and no excitation occurs. 

At the anode the variation of e and 7 is represented by a 
pair of equations similar to equations (7) and (8) but with the 
sign of I reversed- Both e and 7' decrease and for sufficiently 
strong currents may become negative, which is physically 
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impossible if they represent concentrations. The solution of 
the difficulty lies in the fact that equations (7) and (8) are 
only approximate and hold only within a limited range of 
currents. For that matter, the empirically found simple laws 
of excitation also hold only within a limited range of exciting 
currents. When the current at the anode is broken, both e 
andj return to their initial values eo andjo, according to equa- 
tions which are of the form of equations (5). Because of 
k> My the quantity e rises more rapidly, and again for a 
short time it exceeds the quantity j, so that excitation occurs. 

Some other general facts, such as the failure of very slowly 
rising currents to excite, have been shown to follow naturally 
from this “two-factor’^ theory of excitation^ (ifB, chap. xvii). 
Other special cases, however, such as the time relations for 
excitation by alternating currents, were not studied in detail. 

In 1935 A. V. HlO.^ proposed a somewhat different theory, 
which, however, shows much similarity to the above-men- 
tioned one. A. V. Hill and his collaborators have subjected 
different consequences of that theory to thorough experimen- 
tal tests and have found them in good agreement with ex- 
perimental results. 7 ’ ® 

A. V. Hill prefers to talk not of “ion concentrations” but 
rather of two factors, without specifying their physical nature. 
This has the advantage of a greater generality. For that 
matter, nothing would be changed in the author’s theory if, in- 
stead of ionic concentrations, we would refer to € as the 
“excitatory factor” and to j as the “inhibitory factor.” We 
then would not even need to restrict ourselves to positive 
values of e and/, as has been done on pages 109-10. A. V. Hill 
calls j the “accommodation factor.” 

As differential equations for the variation of the two fac- 
tors, A. V. Hill chooses, in our notations, the following: 

de 
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and 

^ = if'(e-€o)-w'(j-_7„), (13) 


with relation (9) as condition of excitation. Moreover, Hill 
imposes a more restricted condition than (ii), n am ely, 


El 

k' 



(14) 


Equation (12) is identical with our equation (7). However, 
equation (13) states that the rate of increase of the fac- 
tor j is proportional not to the current I directly but to 
the excess 6 — e,, of the excitatory factor, created by the 
current. 

The general situation for “make” and “break” of a con- 
stant current is represented in HiU’s theory by a diagram 
essentially similar to that of Figure 27. The only difference 
is that they-curve is not a simple exponential one but begins 
with a zero slope and has an inflection point, as shown on 
Figure 28. 

In view of the great similarity of the two theories, it seemed 
of interest to compare the detailed consequences of both. In 
attempting to do so, F, Offner® found that, as far as intensity- 
time relations for excitation with any kind of current are 
concerned, the two theories are entirely identical. Gale 
Young” proved subsequently that any other two-factor 
theory which is based on a pair of simultaneous differential 
equations, linear in the current I and in the factors e and j, 
and with real characteristic roots, can be brought, by a suit- 
able transformation, into the form given by equations (7) and 
(8). Thus, any experiments on excitation by electric currents 
which are in agreement with A. V. Hill’s equations (12) and 
(13) are also in agreement with our equations (7) and (8). In- 
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asmuch as equations (7) and (8) are independent, while equa- 
tions (12) and (13) are simultaneous, the mathematical 
analysis is somewhat simplified in the former case. As an il- 
lustration, we shall derive here from our equations, following 
A. Weinberg, the relations for excitation by alternating cur- 
rents and for excitation at the anode — relations which have 
been studied theoretically and experimentally by A. V. Hill 



and his collaborators. For a number of further details we 
must refer to the original paper of A. V. HiU and his school, 
covering a very extensive field. Although the fundamental 
idea of the two-factor theory had been suggested by the au- 
thor a few years previously, to A. V. Hill and his school be- 
longs aH the credit for actually elaborating the theory in de- 
tail and subjecting it to rigid experimental tests. 

Let the current I in equations (7) and (8) be of the form 

I~IoSmo)t. (15) 
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Introducing equation (15) into equations (7) and (8), and 
integrating them by elementary methods {MB, p. 157), we 
j&nd 


loK 


€ “ €0 = 


3 -Jo = 


^3 {k sin — 05 cos o)t) ; 


loM 

+ 05 - 


(16) 


{m sin 05 ^ — 0) cos o;^) 


Hence, 
j = 1 0 


Kk 


Mm 


hJ + + 05 ^ 

, Xo5 Jkfo; 

+ 05 ^ 


sm 05/ 


cos 05/ 


+ €0 — Jo • 


(17) 


By a procedure similar to that used before {MB, p. 158), we 
find that the amplitude or the maximum value of e— j is 
equal to 


{^-j\ 

=-f-V 


Kk 


Mm 


^2 ^2 ^2 _j_ ^2 


+ 


iro 5 

+ 05 ^ 


If 05 


+ 05 ^ 
+ Co — Jo 


(18) 


After elementary rearrangements, using equation (14), which 
is merely a restricted form of the last equation (ii), we find 




(K - M)co 
V + co^) (m* + 0)“) 


4 - eo - jo ■ (19) 


As we have seen in MB, page 167, equation (13), equations 
(7) and (8) lead to the following approximate expressions for 
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the rheobase Ro of the cathodic excitation by a constant cur- 
rent 


Ro 


(j\ — eo)k 

K-M ■ 


Hence, 


jo - 6o = 


R,{K - M) 
k 


(20) 


(21) 


Excitation by an alternating current will just occur when 
the amplitude h of the current is high enough to make 

(e — i)max == o . (22) 

Introducing equation (20) into equation (19) and the latter 
into equation (22), we find 



Equation (23) gives us the ratio of the liminal amplitude of an 
alternating sinusoidal current to the liminal value of a con- 
stant current in terms of the frequency co. A. V. Hill has de- 
rived a similar relation on the basis of his theory, expressed 
by equations (12) and (13). In Hill’s notations our k is de- 
noted by i/k, our m is denoted by i/X, and co is equal to 
2Tni, Moreover, Hill uses Iq for our Rc and I for our 7 o. With 
this change of notations equation (23) becomes identical 
with equation (i) of A. V. Hill, B. Katz, and D. Solandt.^ 
Figure 29 shows a comparison of the equation with experi- 
mental data, taken from a paper by these authors. 

We now shall derive from equations (7) and (8) the rela- 
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tion between the duration of a constant current and the 
minimum intensity necessary to produce an excitation at the 
anode upon breaking the current. If I represents the duration 
of the current I from the moment of the “make” to the mo- 



Fig. 29. — Excitation by alternating current of frequency co. Threshold 
current I plotted against log low. The curve is the theoretical one, de- 
rived from equation (23). Circles represent observed values. From A. V. 
Hill , B. Katz, and D. Y, Solandt, Froc. Roy. Soc, London, B, 121, 113, 
1936. 


ment of the ^‘break/^ then, according to equation (6) and the 
remark about behavior at the anode made on page 109, at 
t — i the quantities e — €0 and 7 —/o are given by 


e - €0 = (i - ; 

j- jo = (i - e-™"*) . 


(24) 
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From the moment t = 1 on, these values vary according to 
equations of the form of equation (5) with the right-hand sides 
of (24) as initial values. Hence, 


e — 60 = (i — ; 

;■ -io = -^ (i - 


As we have seen in MBj chapter vii, when excitation occurs 
just at threshold we have de/dl = dj'/dt Introducing equa- 
tions (25) into this requirement and solving the resulting 
equation for we find for the time of excitation after re- 
arrangements, again remembering equation (14), 


/ = 



log 


m € 


)—mt 


~kt 


(26) 


On the other hand, at the moment of excitation we have e = j. 
Introducing expression (26) for t into equations (25), and 
the latter equations into e = j, we find, after some rearrange- 
ments, 

^ ^ El (l - f ^ \ 

° k (j_ ^ ^-mi^k/im-k) \m) \k 

Since m <K k, {m/k) — i — i; and therefore equation 
(27) gives approximately 

k(Jo “ € 0 ) (l — 0-^iy/[i-(m/k)] \^fn 

Again making use of equation (21) and of inequalities (ii), we 
finally obtain 

J (i _ ! ]^y/[{k/m)-i] 

Rc (i — \m 


(28) 



EXCITATION IN PERIPHERAL NERVES 117 

Equation (28), except for differences in notations, is again 
identical with the corresponding equation deduced by A. V. 
Hill from his equations (12) and (13). Figure 30 shows a com- 
parison of that equation with experimental data. 

A problem closely connected with that of excitation is the 
problem of conduction of the nervous impulse. It has been 
rather generally accepted that the cause of the spread of 





Fig. 30. — “Break’’ excitation at end of constant current pulse. 1 rep- 
resents current threshold for excitation at pulse duration i. Points are 
experimental values. The curve is the theoretical one, derived from equa- 
tion (28). From D. Y. Solandt, Proc, Roy. Soc. London, B, 120, 400, 1936. 

excitation lies in the local electronegativity of the excited 
region. This electronegativity with respect to adjacent non- 
excited regions results in an electromotive force which sets up 
local electric currents. These currents act in such a way as to 
excite directly the regions immediately adjacent to the al- 
ready excited one. 

It follows, from this picture, that the mechanism of propa- 
gation of the excitation along the nerve will depend largely 
on the mechanism of excitation by electric currents. It also 
depends on the electrical properties of the nerve fiber. As- 
suming H. A. Blair^s theory of excitation and neglecting the 
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distributed capacity of the nerve fiber, the problem of conduc- 
tion can be treated exactly as has been shown by the author;” 
it requires the solution of a special type of functional equa- 
tion. For the two-factor theory of excitation in the author’s 
form, the problem still can be solved exactly,” though the 
exact expressions are rather cumbersome, and considerable 
simplifications are obtained by some approximations. Where- 
as Blair’s theory leads to a strictly constant velocity of con- 
duction, the two-factor theory requires that the rate of con- 
duction should vary very slightly in the neighborhood of the 
originally excited region, but within about a millimeter of 
this region, there wiU be practically a constant value. The 
experimental detection of such slight variations would present 
some difficulties. 

Alvin Weinberg'3 has studied the problem for the case of a 
nonnegligible capacity of the nerve. In this case only an ap- 
proximate solution has been obtained so far. 

All the above-mentioned theories lead to expressions for 
the velocity of propagation m terms of a number of other con- 
stants of the nerve fiber, such as resistance of myelin sheath, 
rheobase, etc. Most of these constants have, unfortunately, 
not yet been determined with any degree of accuracy. 
Therefore, while we may say that the expressions obtained 
agree in a general way with observations, a really quantita- 
tive and exact verification is still wanting. At any rate, these 
mathematical studies definitely suggest new experimental 
problems. 

A theory of nerve conduction, based on an approach some- 
what different from the above-mentioned ones, has been 
developed in an interesting paper by W. A. H. Rushton.”* 
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CHAPTER VIII 


SOME REMARKS ON THE PHYSICOMATHE- 
MATICAL ASPECTS OF CENTRAL 
EXCITATION AND INHIBITION 

In the third part of our book Mathematical Biophysics we 
attempted to develop a systematic abstract mathematical 
theory of the functions of the central nervous system. The 
fundamental idea was to postulate a few mathematically 
definite laws of interaction between two adjacent neuroele- 
ments and then to consider to what consequences such laws 
of interaction lead when applied to different geometrical ar- 
rangements of the interacting elements. In other words, the 
attempt was made to reduce the tremendous complexity of 
functions of the central nervous system to the complexity of 
its structure, keeping the fundamental dynamic processes as 
simple as possible. A different line of approach is also pos- 
sible, namely, to consider theoretically the structure of the 
central nervous system, with its tremendous number of neuro- 
elements, as quasi-homogeneous and to try to account for the 
enormous complexity of its functions by postulating cor- 
respondingly complex dynamical laws of interaction between 
the individual elements. 

It appears to us that the first point of view may be con- 
sidered as indirectly better Justified by actual observation. 
As far as experience goes, there does not seem to exist any too 
pronounced fundamental difference between the elementary 
physiological properties of the neuron of the frog and those of 
the neuron of man. One would therefore naturally tend to 
explain the enormous difference in the functions of the brain 
of the frog and that of man by the equally enormous differ- 
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ence in the complexity of the structure of their brains. Theo- 
retically, also, the first approach seems to be the more promis- 
ing one. In some earlier publications^ we have attempted to 
develop the second point of view. While some rather inter- 
esting general conclusions were thus obtained, nothing of a 
quantitative nature seems to be gained without a large num- 
ber of additional and rather disconnected assumptions. The 
interesting speculations of W. Kohler/ which implicitly as- 
sume the second point of view, also do not materially con- 
tribute to a quantitative mathematical theory of the brain. 
On the other hand, the first point of view, as developed in 
MBj not only has led to a rather natural systematic mathe- 
matical theory of the brain but has also shown its empirical 
usefulness by leading to quantitative relations that are found 
to be verified experimentally. We refer to H. D. LandahFs 
work on psychophysical discrimination and on reaction 
times, to A. S. Householder’s study of the sensory discrimina- 
tion, and to N. Rashevsky’s work on the aesthetic measure of 
geometrical patterns, to be discussed in subsequent chapters. 

While the primary object of the developments given in MB 
was to work out an abstract physicomathematical theory of 
the brain and to demonstrate the applicability, in principle, 
of the mathematics to this field, the above-mentioned work 
of Landahl, Householder, and Rashevsky does indicate that 
something more than a mere illustration is contained in our 
fundamental postulates. In choosing these postulates, we 
were guided by some general experimental evidence. On the 
other hand, however, the possible connection of these postu- 
lates with some other empirical evidence may not be so clear. 
We shall therefore analyze such a connection somewhat more 
elaborately. 

Taking a suggestion from the weU-known results of E. D. 
Adrian and his school, we first assume that under physio- 
logical conditions a constant peripheral stimulus of any nature 
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results in a sequence of impulses in each single afferent nerve 
fiber. The intensity I of each impulse is assumed to be con- 
stant^ thus preserving the “all-or-none” law, while the fre- 
quency V of the impulses is taken to be an increasing function 
of the intensity S of the stimulus. For small values of S the 
frequency v is considered as approximately linear: 

V = a(5 — h) , (i) 

where h is the threshold below which a stimulus does not 
elicit any excitation and where a is a constant. Equation (i) 
is understood to be only an approximation to some more gen- 
eral type of equation, in which v tends to a constant value 
with increasing 5 , as must be the case from considerations of 
the refractory state of the fiber {MB), Such a more general 
equation may be, for instance {MBj p. 218), 

X, = i [i — ^ 

u 

where € is the refractory period; or, as another possibility, 

5 

V = ah log = ; 5 « i , (2) 

x + *f 

which, for small values oiS/h, reduces to 

V = ah log ^ . 

Defining the intensity, E, of excitation of a fiber by 

£ = /., (3) 
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which, together with (i), gives 

E = aI{S - h) , (4) 

we now establish the fundamental equations of interaction be- 
tween adjacent neuroelements. 

The axon end of every neuron is assumed to produce two 
factors — one excitatory, e; and the other inhibitory, j. It 
must be emphasized that we do not consider these factors as 
substances. They may be such, but they may represent any 
kind of physicochemical state having the property of exciting 
or inhibiting the neuron upon which it acts. We assume that 
the quantity e— j affects any neuron in the same way as does 
a natural stimulus; that is, e andy, produced by the ith neuron 
at the synapse with a ^th neuron, result in an excitation of 
the latter with an intensity given, in analogy to (4), by 

Ek = aklki^ — j ~ hk) . (s) 


Excitation of the Ath neuron occurs only if e — > hk- 

The equations necessary to complete our set of postulates 
specify the production of e andj by a fiber. They are 


i-BE-ij. 


( 6 ) 


Equations (i)-(6) are sufiScient to describe quantitatively 
any situation arising from the interaction of any number of 
neurons arranged according to a definite geometrical pattern. 
Provided we can master the purely mathematical complexity 
of any such problem, we can proceed to the consideration of 
various forms of neuronic circuits and their functions. 
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Equations (6) are formally identical with equations (7) 
and (8) of the preceding chapter. When a stimulus S, of con- 
stant intensity, is suddenly established and kept indefinitely, 
then, according to equation (4) or to a corresponding, more 
general one, the result is the sudden establishment of a con- 
stant E. And following this, there is, as in chapter vii, an 
exponential increase of e and 7, both tending asymptotically 
to limiting values, corresponding to AE/a and to JBE/b. The 
larger the a or b, the more rapidly these limiting values are 
approached. 

We do not need to make any definite assumptions about 
the constants A, B, a, and b in equation (6). The following 
cases, however, present special interest. 

If 

A>B; a>b-, (7) 

then for a constant stimulus of indefinite duration the fiber 
always exerts an excitating action on the other fibers with 
which it synapses, with e being produced in excess of j. If, 
on the contrary, 

A <B a<b; -<t (8) 

a b 

then the fiber is an inhibitory one (MB). 

Whenever the constants satisfy the relations (7), so that 
the nerve is an exciting one, we shall attach to those constants 
the index e. W^e shall have, in this case, -Be, and When- 
ever the constants satisfy relations (8), we shall mark them 
Aj, Bj, Gj, and bj. A neuron producing only e shall be called 
“purely excitatory”; one producing only 7, “purely inhibi- 
tory.” While, in a sense, we thus assume the existence of 
special inhibitory fibers, an assumption which does not meet 
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at present with general recognition (see J. C. Eccles,^ R. 
Lorente de No,^ and H. S. Gasser-i), yet, at the same time, ac- 
cording to this theory, an inhibitory fiber diff ers from an 
excitatory one only quantitatively, the two being otherwise 
fundamentally identical in their general behavior. Other rela- 
tions between the constants A, B, a, and h lead, moreover, to 
situations in which one and the same fiber may act either as 
an excitatory or as an inhibitory one, depending on the time 
relations of the stimuli {MB). 

We may also have cases intermediate between those repre- 
sented by relations (7) and (8) ; namely, we may have 

A < B with a<h and — > w ; 

a 0 

A B 

A > B with a > b and — < — . 

a b 

The first case gives for a constant, continuous stimulus a 
brief inhibition, followed by a lasting excitation {MB, p. 
224). The second case gives, for the same conditions, a brief 
excitation, followed by a lasting inhibition. 

While, as has been remarked, the above-mentioned postu- 
lates lead, in a number of cases, to consequences which are in 
quantitative agreement with experimental data, and while 
such an agreement may be considered as justifying them, to 
some extent, stiU, at first sight, these postulates do not fit 
quite well into the picture made by neurophysiologists of the 
basis of some electrophysiological and other experiments. 

First of all, the situation as represented by our equations 
(i)-(6) is undoubtedly highly oversimplified, as compared 
with the actual situation. The wealth of observations re- 
viewed, for example, by Fulton or by Eccles^ reveal a maze of 
complexity in the interaction of central neurons. While, ac- 
cording to our equations, the elementary process of excitation 
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and conduction at a sjoiapse appears to be a rather simple one, 
we know that it is actually decomposable into a rather large 
number of factors. 

It must, however, be kept in mind that, from the point of 
view of the theory of interaction of the very large number of 
neurons which constitute the central nervous system, the de- 
tails of the mechanism of interaction of two neurons may 
be quite irrelevant as long as we have a quantitative expres- 
sion which describes the net result of this interaction. And, 
while the details may be very complex, these net results may 
be described by relatively simple equations. An example from 
physics may best illustrate the point. In the theory of heat 
conduction it is assumed, on the basis of some generalizations 
of experimental observations, that the heat transport in a 
given direction is proportional to the gradient of temperature 
T in that direction. Thus, choosing the direction of flow Q as 
that of the ir-axis, we have 

e = -« , (i°) 

where k is the thermal conductivity. On this simple relation 
is based a mathematically most complete and practically very 
useful theory of heat conduction in bodies of various shapes. 
Yet the actual mechanism of heat transport is a very complex 
one, determined by the detail of the thermal motion of the 
molecules. The constant k is a constant only to a first ap- 
proximation, and even then it resolves itself into a number 
of factors, expressed in terms of gas-kinetic parameters. The 
complete and exact expression for the thermal conductivity k 
is not available even now, though we have a pretty good in- 
sight at present into the mechanism of heat transfer. At the 
time of Fourier, who founded the mathematical theory of 
heat conduction, practically nothing was known of the com- 
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plexities of the kinetic theory. Yet, without waiting for the 
detailed mechanism which underlies equation (10) to be 
worked out, the theory of heat flow has reached a high degree 
of theoretical and practical perfection. 

Similarly, the simple equations (5) and (6), which form 
the essence of our theory of central excitation and inhibition, 
may be used for a development of the theory without waiting 
for a detailed physiological interpretation of the underlying 
mechanism. Unquestionably, these equations are very crude 
and approximate, but so is equation (10) of the heat transfer. 
Yet for many practical purposes they may be accurate 
enough, as the above-mentioned work of Landahl, House- 
holder, and Rashevsky indicates. The factors e and 7 may 
each be of a rather composite nature, just as is the coefficient 
K in equation (10). A physicochemical interpretation of these 
factors must be hoped for in the future. At present we may 
justifiably proceed on the more formal level. From this point 
of view it may be proper to caU the developments of the third 
part of MB and such developments in subsequent publications 
“mathematical neurology” rather than “mathematical bio- 
physics of the central nervous system,” to emphasize the 
lack, at the present stage, of any physiological and physical 
interpretation. 

It has been emphasized recently, on many sides, that the 
inhibitory effects are more likely to be due to an interaction 
of specially arranged excitatory mechanisms than to specific 
inhibitory impulses. Some such mechanisms have been sug- 
gested, for instance, by H. S. Gasser."* The net result of the 
functioning of such a scheme or any other similar one is the 
creation of an inhibitory state, which gradually disappears. 
If the rate of establishment of this inhibitory state is ap- 
proximately proportional to the intensity of excitation of the 
neuron responsible for setting in action the whole mecha- 
nism, and if its decay, in the absence of disturbances, is ap- 
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proximately exponential, things will occur as if they were 
described by equations (5) and (6). The 7-factor in this case 
may stand for a very complex excitatory state of a rather 
complicated anatomical structure. All that matters is the end 
result. 

It must be emphasized that the element of “as if” enters 
into mry physical problem described by differential equations. 
A differential equation deals with infinitesimals at a mathe- 
matical point, and, while it may be suggested by experiments, 
it is actually beyond any direct experimental verification. In 
many cases the most important differential equations of 
physics have only a statistical meaning and are known to 
break down for very small regions. Such is the case, for in- 
stance, with the differential equations of hydrod3mamics and 
aerodynamics. Concepts like pressure, density, etc., which 
enter as continuous functions into these equations, have a 
meaning only for volumes which contain a very large number 
of molecules. For a volume of io~^* cm^ these concepts break 
down, and the differential equations do not hold at all. Yet, 
with very great accuracy, things actually happen as if the 
differential equations would hold down to infinitesimal 
regions. A similar interpretation is possible for our fundamen- 
tal postulates (i)-(6). In the future they must be reduced 
to the details of physicochemical mechanism, just as equation 
(10) is now being reduced to kinetic interpretations. There 
is no doubt that a more detailed interpretation of our postu- 
lates wiU result in their modification and improvement, and 
a study of possible modifications is definitely indicated. In 
the meantime, however, we can and should proceed with the 
development of the theory based on these postulates. 

From these general remarks we shall now pass to the dis- 
cussion of a few specific questions of the relation of our postu- 
lates to some current neurophysiological views. 

The intensity E of excitation of a fiber, as defined by (3), 
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can have only a statistical significance. When we are dealing 
with only one, or a very few, irregularly spaced impulses, such 
as are used in electrophysiological exjjeriments, the concept 
of frequency, v, loses its meaning. 

Considering continuous stimuli, in which case (3) can be 
used, we find (MB, chap, xxiv) that the synaptic delay be- 
tween iwo neurons decreases with increasing intensity of the 
stimulus. In the general case the synaptic delay is finite at 
threshold, decreasing to zero with increasing 5 . This latter 
decrease, however, will occur only if we consider equation (i) 
as holding for any value of S. If, for higher values of S, we 
adopt equation (2), then, as is readily seen from the argument 
in MB, increasing S decreases the synaptic delay only down to 
a finite value. Thus the synaptic delay varies within fixed 
limits; and by a proper choice of constants we may make 
those limits sufficiently narrow so as to be in agreement with 
the experimental findings of R. Lorente de No. 

Let us now look at the picture from a different angle. Sup- 
pose that the transmission at the synapse is due to the action 
current of the presynaptic impulse. Two possibilities must 
then be considered: (a) the intensity of a single “all-or-none” 
impulse is subliminal; (b) the intensity of a single impulse is 
superliminal. 

In case (a) a volley of impulses of frequency v will produce 
a pulsating excitation current at the S5mapse which will act 
much as an average equivalent direct current of constant in- 
tensity would. Considering, for simplicity, the shape of the 
action current of an impulse as rectangular, we can calculate 
the excitation time of the neuron at the synapse by using the 
equations developed in MB, pages 220-26. If we assume that 
the law of excitation of the neurons is the same as that of 
peripheral nerves, based on the HiU-Rashevsky two-factor 
theory, we shall find for the S3maptic delay formally the same 
expression as before. 
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In case (b), however, the synaptic delay between two neu- 
rons will be determined solely by the intensity of the indi- 
vidual impulse. It will be quite constant for a given pair 
of neurons and will not vary with either E or S. If, however, 
several fibers of different thresholds form synapses with the 
body of the same neuron, then increasing the intensity of the 
stimulus increases the number of excited fibers, and the body 
of the neuron is acted upon by the sum of several impulses, 
with a corresponding reduction of the synaptic delay. This 
seems to represent the case studied by R. Lorente de Nb.^ 
WHle his experiments give strong evidence for this interpre- 
tation, there is no reason known, as yet, to assume that case 
(a) may not also be present in some parts of the central 
nervous system. R. Lorente de No’s observations refer only 
to motoneurons and do not preclude a variability of the 
s}maptic delay between two neurons of some other types, es- 
pecially because the two cases differ only in the quantitative 
values of the parameters of the neurons. 

Considering, now, the inhibitory effects, we may well as- 
sume that the inhibition is merely the result of a refractori- 
ness produced by proper excitation. Let, for instance, a par- 
ticular case be characterized by a very large a (Eq. [i]), so 
that even for very weak stimuli, 5’, the frequency, v, will be 
so high that the successive impulses will fall within an inter- 
val of time smaller than the refractory period of another fiber 
with which the first one synapses. Then for average stimuli 
.S', applied to the first fiber, the second fiber will be main- 
tained continuously in a refractory state and thus may be 
considered as inhibited. It would be only natural to assume 
that, in general, the inhibitory state of the second fiber, as 
measured by the increase of the threshold, will be the stronger 
the higher the frequency v of the first, "‘inhibiting” fiber. At 
the same time, it is not likely that this state of inhibition, or, 
if we prefer, of “relative refractiveness,” will be established 
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at once, with the first impulse of the first fiber. If this state 
develops gradually at a rate approximately proportional to 
the frequency v of the first fiber, and if, when left to itself, 
it decreases exponentially with respect to time, then, mathe- 
matically speaking, this refractory state will have, approxi- 
mately, all the characteristics of our j-factor, as described by 
equation (5) and ( 6 b). It must be remarked that equation 
(5) may be written thus: 

Ek — "" (hic "by)] , 

with 

e> h +j (12) 

as a condition of excitation. In this form the equation merely 
says that the effect of j is to increase the threshold of the 
other neuron. 

If the first fiber s3mapses with a fiber which has a very 
short refractory period, then the former will excite the latter, 
producing in it a continuous volley of impulses of frequency v. 
Thus, according to this conception, the same fiber may act as 
an excitatory or as an inhibitory one on another fiber, de- 
pending on the relative values of the physical constants of 
the two fibers. 

Our equations (i)-(6) have been patterned very much 
along the lines suggested by the two-factor theory of periph- 
eral excitation. In our original form® we called the factors 
“excitatory” and “inhibitory.” In the slightly modified but 
essentially equivalent form of the theory proposed later by 
A. V. Hill® our inhibitory factor is called “accommodation.” 
The name does not seem to matter much, since the effects of 
the increase of accommodation are opposite to those resulting 
from the increase of excitation, and so accommodation is es- 
sentially an inhibitory factor. The experimental evidence ob- 
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tained by A. V. Hill and his school has put the two-factor 
theory on a firm empirical basis. It is, therefore, quite in 
keeping with the prevailing tendency to reduce central excita- 
tion and inhibition to laws of peripheral nerves to adopt here, 
also, some sort of a “two-factor” theory, as we have done. 

One might even go farther and identify the central e and 7 
with the peripheral ones. We may assume that e and J pro- 
duced at the end of the axon act directly on the adjacent 
neuron. In the theory of peripheral excitation we have con- 
sidered (cf. MB) as a condition of excitation 



and for simplicity have chosen the units so as to make h = 1, 
which reduces (12) to 


This is admissible as far as a single fiber is concerned. But, 
in considering several fi.bers, we must take into account the 
possibility that h varies from fiber to fiber, in which case we 
cannot consider it as always equal to 1. If a fiber with a 
given h synapses with another, which has an V > h, then 
while € = Ay in the first fiber produces excitation, it fails to 
do so in the second one, because there we must have 
€ = h'j > hj. Inasmuch as, in the two-factor theory of periph- 
eral excitation, e first exceeds j and then again drops below 
it, thus creating a state of inhibition which for properly 
chosen constants may be stronger than the initial (jo — 6o) 
(MB), several interesting possibilities are open here. 

* The h in (12) is not to be confused with the h in equations (i) and (5) 
(cf. MB, chaps, xvii and xxii). 
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CHAPTER rX 


MATHEMATICAL BIOPHYSICS OF SOME SIMPLE 
NEUROLOGICAL STRUCTURES: APPLICATIONS 
TO REACTION TIMES 

Having established in the previous chapter the fundamen- 
tal equations which govern the interaction between neurons, 
we shall now follow the procedure clearly indicated. We 
must systematically study different, more or less complex, 
geometrical arrangements of neurons and mathematically de- 
rive their properties. We shall start with the simplest possible 
cases and then gradually complicate the picture. 

Z JT 


Fig. 31 

The simplest possible structure is that of a neuron with a 
peripheral fiber synapsing with a neuron of higher order. As 
we have seen in the previous chapter, we do not need to enter 
into the considerations of the detailed structure of the synapse 
in order to apply our fundamental equations. We shall there- 
fore represent two synapsing neurons schematically, as shown 
in Figure 31, simply by the adjacent lines. In MB we indi- 
cated by arrows, placed at the synapse, the direction of trans- 
mission. Inasmuch as a forklike scheme is used in neuro- 
logical diagrams to represent schematically the branchings 
of the axon at the synapse and a point is used to represent 
schematically the body of the higher-order neuron, our former 
notation appeared somewhat confusing to the neurologists,' 
for in the classical neurological notations the nervous impulse 
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travels from the fork of one neuron to the point of the other, 
that is, in the direction opposite to that assumed in our nota- 
tions, since a fork can be looked upon as an inverted arrow. 
Therefore, we shaU now place the arrows which indicate the 
direction of transmission in the middle of the lines which rep- 
resent the neurons, thus avoiding confusion. It must be 
noted, however, that some neurologists use the same conven- 
tion as we used in MB.^ 

We may now consider various stimuli, S, which are differ- 
ent with regard to their variation with respect to time, as 
being applied to the afferent end of the neuron I (Fig. 31). 
We shall confine ourselves here to the consideration of a 
stimulus of constant intensity 5 i, applied suddenly and kept 
indefinitely. 

In this case neuron I will be subject to a suddenly estab- 
lished intensity of excitation connected to S-^ by one of 
the equations discussed in chapter viii. From the moment the 
excitation arrives at the S3mapse, the two factors e andj will 
vary according to 

€ = — (i - ; J = — (i - . (i) 

As soon as € — j reaches or exceeds the threshold ki of the 
neuron II (Fig. 31), the latter becomes excited, and the 
excitation is transmitted farther along it. The time which 
elapses between the moment of arrival of the excitation 
along the fiber I at the synapse and the moment when the 
neuron II becomes excited — ^in other words, the synaptic de- 
lay — is obtained by solving, with respect to t, the equation 
e — j = ha, after introducing into the latter the values for e 
and j from equations (i). The details of the calculations are 
given in MB, chapter xxii. The relations become particularly 
simple when B. is very small— -in other words, when neuron 
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I produces very little of the inhibitory factor j, so that the 
latter may be neglected, as compared with e. In this case the 
synaptic delay is given by {MB, p. 220, Eqs. [12] ff.) 



AeE, 


(2) 


The relation for the general case is more complicated, though 
similar in general, to the foregoing {MB, chap. xxii). 

Equation (2) shows that the stronger the intensity of exci- 
tation Et the shorter the synaptic delay 4. Since Ex increases 
with increasing intensity 5 i of the peripheral stimulus, there- 
fore the stronger the stimulus Sx the shorter should be the 
synaptic delay 4- To obtain a relation between 4 and S,, we 
must use one of the equations of the previous chapter. As- 
suming relation (2) of chapter viii to hold, we find, approxi- 
mately, for not too large values of Sx, 



log Sx log Jtx 

log Sx Aelxdihx log Hx * 


(3) 


in which the subscript i indicates that the constants refer to 
the neuron I. 

Putting 


Ex = log hx ; 




Q'shx 

■AgJ xdixHl 


we obtain equation (3) in the form 


4 


AW. JggA- gA_ 

a, log Sx- Ex- H,‘ 


(4) 


Using the approximate expression (4) of chapter viii, we 
find for 4 an expression of the form {MB, chap, xxii, Eq. [13]) 

< = — locr AeaiIi{Sx — hx) 

‘ a, ^ AeaxIx{Sx - hx) - aA ' 


(5) 
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Consider a chain of neurons transmitting the excitation 
successively to each other, and remember that the velocity of 
propagation along each individual link of the chain is inde- 



Fig. 32. — Reaction time for visual stimuli, plotted against stimulus 
intensity. The curves represent equation (4). The circles represent ex- 
perimental values by G. 0 . Berger and J. McK. Cattell for two different 
subjects. The following values for the constants in equation (4) were 
used. Subject B: i/ae = 0.413 sec; Hi ^ — 1.96; E2 = 0.57. Subject 
C: i/ae — 0.279 sec; Hi = 1.57; Hz — 0.57. From H. D, Landahl .3 

pendent of the intensity of excitation. In that case a qualita- 
tively similar relation between the intensity of the peripheral 
stimulus and the total time of transmission along the chain 
must hold, that is, the stronger the peripheral stimulus the 
shorter the total time of transmission along the chain. We 
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cannot expect equations (4) and (5), whidb were derived for 
a single synapse, to hold quantitatively in that case. However, 
in one particular case equations (4) or (5) may hold even for 
such a rViain of neurons — ^while not exactly, yet with good ap- 



Fig. 33. — Same as Fig. 32 for a set of auditory data by H. Pi^ron. The 
values of the constants are: = 0.304 sec; Ei ^ — 0.37; ila = 

0.23. From H. D. Landahl .3 

proximation. This will occur when all S3niapses but one in the 
chain have, for given conditions, approximately equal and 
very short delays, 4, while the one exceptional synapse (or per- 
haps a very few synapses) has a much larger synaptic delay. 
Then the total time of transmission along the chain will be 
controlled mainly by the one “long” synapse, and the relation 
between stimulus intensity and total transmission time will 
be of the form given either by equation (4) or by equation (5). 
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The total time of transmission is then obtained by adding to 
ts a constant time, ta, which represents the duration of the 
transmission along each fiber, delays at end-organs, etc. 

In view of these facts, it may be of interest to compare equa- 
tions (4) or (5), with a constant term added, to experimen- 



Fig. 34. — Reaction times for gustatory data, plotted against intensity 
of stimulus. The curve represents the theoretical equation (5). The 
circles are experimental values by H. Pi6ron. The values of the param- 
eters in equation (5) are, in this case: 1/% = 3.8 sec; hi = 0.6; 
aJhlAgo-ih — 0.2. From H. D. Landahl.J 


tal data for reaction times. Such a comparison has been made 
by H. D. Landahl,^ and the results are shown in Figures 32, 
33, and 34. Visual and auditory data (Figs. 32 and 33) are well 
represented by an equation of the form of equation (4), while 
the simpler equation (5) is sufficient to represent the gusta- 
tory data. For discussion of details we must refer to the 
original paper by H. D. Landahl.3 

A more complex neuronic structure is represented in 
Figure 35. Here a neuron 11 of second order is excited by two 
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neurons, I and III, of the first order. We may consider cases 
where the stimuli Sx and ^3 are applied at different times. Of 
particular interest is the case where ^3 is too weak to produce 
a response in neuron II and is always applied before Sx, Fol- 
lowing H. D. Landahl, we shall consider a special case, name- 
ly, where neuron I produces only the factor e, while neuron 
III produces both e and 7 but is of the excitatory type (chap, 
viii, expressions [7]). Later on we shall discuss a more general 
assumption regarding neuron L The constants a, B, and b 
are considered as different for the three neurons and are de- 
noted by corresponding subscripts. It is assumed that 
A^/a^ = B3/53, although A^>Bi and a^>b2,- 



Fig. 35 


If a constant stimulus Sx of indefinite duration is applied 
suddenly to the neuron J, then at the synapse the amount €i 
of € produced by that neuron is given by 


= 



(i — , 


( 6 ) 


the time tx being coimted from the moment of arrival of the 
excitation of I at the synapse. Denoting by ^3 the time from 
the moment of arrival of the excitation of neuron III at the 
synapse, we have, similarly, for the variation of the amounts 
€3 and 7*3 of the factors e and 7 produced by neuron III: 




(i — 


h 


— ±2 


B^E^ 


(i 


e^hh) . 


(7) 
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Excitation of neuron II occurs when €i + ej — = h^. Ac- 

cording to our assumption, however, €3 — remains always 
less than Let us denote by the time between the arrival 
of the excitation of the neuron I at the synapse and the initia- 
tion of excitation in neuron II. In other words, is the value 
of ti at the moment when neuron II becomes excited. Let 
be the value of at this moment. (Remember that ti and t; 
are counted from different origins, since St and are not 
applied simultaneously.) With the foregoing assumption, tj 
is obtained by introducing expressions (6) and (7) into 
ei -f 63 — 73 = hi, putting ti = and ^3 = ^3, and solving 
the resulting equation with respect to This gives 


ft = 




(i - 
03 


Ml 

h 



>(8) 


The total time U between the application of stimulus S, and 
the final reaction of the end-organ is obtained under similar 
assumptions, as before, by adding to a constant to. 

Of particular interest is the study of the relation between 
the reaction time U and the time l'^ for the case where the in- 
tensities of both stimuli are kept constant and where the 
conduction time along neuron I and III is very short, as com- 
pared with either or In other words, we consider that 
most of the constant time to is due to conduction on the 
efferent side and to delays at the end-organs. In that case, 
approximately represents the time between presentation of 
the stimulus and the beginning of the reaction less the 
time to. 

These considerations can be applied to some experiments 
on the effect of a warning or preparatory stimulus upon the 
length of the reaction time. In such experiments the subject 
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is to react to a given stimulus; but some time before that 
stimulus is applied, a different preparatory or warning stimu- 
lus is given. We may tentatively identify our stimulus Si with 
the stimulus to which the reaction takes place, while may 
be considered as the warning stimulus. In practice the inter- 
val iw between and is always much larger than the re- 
action time U. In other words, 

L = ti- i'l » tr = t'i + io. 


In that case we have approximately 


Putting 

; M = 




aih 




we obtain from expression (8) 


tr = i'lA- ta — U — “ log [M + J{e — e “sVi)] ^ (lo) 

dj 


which is the relation between the preparatory interval and 
the reaction time U. In Figure 36 are shown two sets of actual 
data, 3 compared with curves represented by equation (10). 

It is of interest to note that the value of i/^e is of the same 
order of magnitude for the following curves (Fig. 36) as for 
the visual data (Fig. 32) and the auditory data (Fig. 33), al- 
though there is no apparent connection between the curves. 
However, it should be noted here that the intensity-time 
curves are largely determined by the quantity Hi the 

value h, and the product of H, and i/(Z«. As long as Hj/oe 
remains constant, changes in Ha or in i/oe do not appreciably 
affect the shape of the curve. Then, making the value i/oe 
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the same for two different situations imposes a condition upon 
H2. If, then, any one is determined independently, the 



Fig. 36. — Reaction time plotted against the interval between warning 
and final stimulus. The curves represent equation (10). The points and 
circles represent observations by H, Woodrow. The two sets of data were 
obtained on the same subject but under different conditions. For the 
case represented by circles the subject did not know the lengths of the 
warning period, whereas in the other cases the subject was given practice 
with a particular warning period, after which his responses with that 
warning period were recorded. The values of the parameters are as fol- 
lows: Upper curve: to == 0.13 sec; i/ai = 0.2 sec; M = 0,638; / = 2; 
G3 = 0.132 sec”^; &3 = 0.12 sec“^. Lower curve: to “ 0.13 sec; i/gi = 
0.2 sec; M — 0.67; J = 0.41; — 1.28 sec"^; bs = 0.128 sec"^ From 

H. D. Landahl.3 

others are fixed. Thus we see that the curve of equation (4) 
is, apart from the constant term added, practically a two- 
parametric curve. 
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The peculiar minimum of the tr, ty, curves is due, in the 
present theory, to the fact that the constants of the neuron 
III satisfy the relations (7) of chapter viii. As has been 
shown in MB, page 227, in this case the difference €3 — 
first increases with then reaches a maximum, and finally 
decreases, tending asymptotically to a constant value posi- 
tive for Ajae > Bjie, or zero for Ay/ae = Be/he. Thus, for 
a certain value of a maximum amount of e, — is added 
at the s}mapse. Since the condition of synaptic transmission 



Fig. 37 


is €i -h 63 ~ 73 = h, the larger the €3 — j^, the smaller is the 
necessary amoimt of e, and the sooner will this necessary 
amount be reached according to equation (6). 

If we consider the more general case where neuron I also 
produces both e and j, things become more complicated. 
However, approximately the same relation will hold for the 
case where is suf&ciently larger than the threshold neces- 
sary to produce any reaction at all. In this general case the 
difference a, — of e and j produced by neuron I will vary 
with respect to t't, as shown in Figure 37. When and, there- 
fore, El are near threshold, then the excitation of neuron 11 
begins at the time t'l (Fig. 37, dotted line). When Si and Ei 
are large, then neuron 11 is excited at the time t'l (Fig. 37, full 
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line). But in this case the segment OA of the curve may be 
represented, with sufficient accuracy, by equation (6). 

From such simple neuronic structures as are considered 
here, we may pass to more complex ones. A particularly inter- 
esting structure is that consisting of two purely excitatory 
neurons forming a closed circuit, so that one neuron excites 
the other, and the latter in its turn excites the first. As has 
been shown in MB, chapter xxiv, and later on in more detail 
in a paper by A. S. Householder, 4 such a circuit possesses the 
property of remaining unexcited as long as any external 
stimuli do not exceed a certain threshold. When, however, 
an external stimulus exceeds a given threshold, then the cir- 
cuit is brought irreversibly into a permanently excited state, 
in which it remains even after the removal of the external 
stimulus. Applications of such structures to the theory of 
conditioned reflexes have been made.® 

A. S. Householder and H. D. LandahF have studied more 
complex circuits, involving neurons producing both e and j, 
so that they are not purely excitatory. In that case the perma- 
nent excitatory state of the whole system is not necessarily 
constant, as in the former case. The intensity of excitation 
may oscillate with respect to time. We thus have a system, 
based on our fundamental equations, which possesses a 
property of producing spontaneously, without external stimu- 
lation, a rhythmic excitation. Possible applications of such 
considerations to the theory of the spontaneous rhythmic 
activity of the brain^ may be indicated. 

It has been pointed out in MB (chap, xxii, p. 234) that, 
when we have a chain of neurons with different thresholds 
and other constants, the intensity of excitation of the indi- 
vidual neurons of the chain may either increase or decrease 
along the chain. LandahF has shown that a rather simple 
relation between thresholds of the individual neurons leads 
to exponentially increasing or decreasing excitation. House- 
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holder^ has studied more complex cases. If the end-link of a 
chain is an inhibitory neuron, the chain has a net inhibitory 
effect which may either increase or decrease with the length of 
the chain. 
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CHAPTER X 

DISCRIMINATION OF INTENSITIES 

Let a stimulus of intensity 5 i be applied to a sense organ. 
Under actual physiological conditions, perceptible stimuli, no 
matter how weak and how sharply localized, excite not one 
but a large number of peripheral fibers. A nerve innervating 
a given sense organ usually consists of a large number of 
fibers with different thresholds. A stimulus of a given in- 
tensity, ^i, excites only a fraction of those fibers, namely, 
those whose thresholds are less than Si. If we apply to the 
same sense organ a stronger stimulus of intensity, S2 > Si, 
then that second stimulus will excite all the fibers which were 
excited by the first one, plus an additional number of fibers, 
namely, those whose thresholds lie between Si and S3. An 
excessively strong stimulus, such that its intensity exceeds 
the highest threshold of the bundle of nerve fibers in the 
nerve trunk, wiU excite all the fibers. 

Whether each peripheral fiber is connected by a chain of 
neurons with a corresponding single fiber of an effector end- 
organ, or whether through branches and collaterals this 
peripheral fiber becomes connected to several fibers of the 
effector end-organ, thus producing a sort of “multiple re- 
sponse,” the foregoing considerations lead to the conclusion 
that a response due to a weaker stimulus is, so to say, always 
“contained” in the response for any stronger stimulus of the 
same type. In this simple scheme, whenever a stimulus. Si, 
produces a reaction, i?,, then a stronger stimulus, S2 > S,i, 
produces necessarily such a reaction, i?j, that it includes the 
reaction R,. However, the fact that we can actually dis- 
criminate between stimuli of different intensities shows that 
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the situation is not simple. We respond to a weaker stimulus 
by calling it “weaker,” and to a stronger stimulus by calling 
it “stronger.” The words “weaker” and “stronger” them- 
selves constitute certain motor reactions of the lungs, pharynx, 
and tongue to the corresponding stimuli; and neither of these 
reactions “includes” the other in the above-mentioned sense. 
A still better illustration is obtained by considering cases 
where we are taught to respond to a weaker stimulus in a 
way qualitatively di fferent from our response to a stronger 
one. Thus, a weak sound of a given pitch may be used as a 

C 

\ 

ne Si 

Fig. 38 

signal for beginning to read a book, while a stronger sound of 
the same pitch may be used as a signal for an entirely dif- 
ferent action. 

Hence, some neurological mechanism must exist which 
provides for the possibility of a qualitatvoely different response 
to quantitatively different stimuli. Such a mechanism can be 
provided by the existence of inhibitory fibers, intercalated in 
a proper way between excitatory ones. 

Let us consider the foUowing structure (Fig. 38) and let us 
confine ourselves to stationary states — that is, we shall con- 
sider constant stimuli of sufficiently long duration so that at 
aU synapses the factors e andj have practically reached their 
constant limiting values. In that case e is proportional to E 
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(cf. Eq. [i] of chap. ix). A sensory peripheral fiber or a chain 
of fibers eventually divides into several branches. Each of 
these branches synapses with a neuron ne of higher order. 
The thresholds h of the neurons ne are, in general, different 
and are distributed according to some distribution function, 
N(h). That is, if we take h as abscissa and plot the number 
N{h) of neurons ne having a threshold between h and h + dh, 
where dh is a very small fixed quantity, we shall obtain a 
curve of some sort, which represents N{h) as a function of h. 
If, as in MB, chapter xxii, we consider the simpler case where, 
whenever a fiber divides into several branches, the intensity 
E of excitation in each branch is the same as in the original 
fiber, then we have the following situation. 

For a given intensity S of the peripheral stimulus there 
corresponds a definite value £i of intensity of excitation in 
each branch {MB, chap. xxii). If the chain consists of only 
one fiber, then Ej is a linear function of 5. In general, how- 
ever, Ex is a more complicated function of S, depending on 
the constants of the individual members of the chain 


Ex = P{S ) . (i) 

The intensity of excitation E„e of each neuron ne is given 
by {MB, chap, xxii) 


Ene = ^{PEx - h) . (2) 

The total intensity of excitation of all neurons {ne) which 
have a threshold h < PEx is then given by 

E{h) = ^{PEx - h)N{h) . (3) 

If iV(o) = o, as should be expected from general considera- 
tions, then E{h) is zero for A = o and for h = PEx and is posi- 
tive within this interval. Hence, in that interval E{h) has 
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at least one maximum. Let us first consider the case where 
it has only one maximum, h = hm, hm being a fimction of E,. 

Let all neurons ne send fibers to a region of the brain, Ci, 
in which all neurons ne with the same threshold h form a 
synapse with the same neuron of third order ne^-, and let, 
moreover, each neuron ne of threshold A* send off a collateral 
which excites inhibitory fibers ni, leading to all the synapses 
Si corresponding to neurons of a different threshold hi. 
Consider, for simplicity, that all thresholds hi of ni are the 
same. Each S3niapse Si is excited by an amount of (e — j)t 
due to neurons ne and proportional to the quantity E{}i) in 
equation (3) (cf. chap, ix, p. 135). Moreover, each of these 
synapses receives from neurons ni a certain amount of 
(/ — e)2 > 0. Since each neuron ni is excited by a neuron 
ne, its intensity of excitation is stronger, the stronger the in- 
tensity of excitation of the corresponding neuron ne. The 
few S3mapses, j<, corresponding to such a value of h that E(A) 
is very high, will receive a large amount of (e — _/)i from the 
ne neurons and relatively lesser amounts (_/ — from the 
inhibitory fibers, coming from other less excited neurons ne. 
The total amount oie— j will therefore be positive and large 
enough, and the corresponding neurons ne^ will be excited. 
But the large number of synapses, corresponding to such 
values of h for which EQi) is much smaller will receive little 
(e — j)i from their ne neurons and a large amount of (7 — e), 
from the strongly excited neurons ne. As a result of this, 
only those synapses which correspond to sufficiently large 
values of E{h) will transmit excitation to neurons ne^. But 
such E{h) has a maximum for a value hm of h^, and therefore 
only those synapses Si will transmit excitation that cor- 
respond to values of h, that lie ia the neighborhood of hm, and 
that are therefore included between two fixed values h^ and hi 
(Fig. 39). Since, according to equation (3), Am is a function of 
El and hence also a function of S, the Ai and K are also func- 
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tions of S. If, by varying the intensity S of the peripheral 
stimulus, we vary £„ this wiU result in a variation of h^; and 
if, for a new value E[ of Ei, the corresponding will be suf- 
ficiently different from the A™, then entirely different groups 
of synapses Si will be excited by the stimulus S' t-ha-n by 
S. Thus, to any intensity S of the peripheral stimulus 
there is a corresponding excitation of a definite group of 
S3mapses Si. A stimulus of intensity 5 will produce a reaction 
R through a group Si of s3mapses, while a stimulus of a differ- 
ent intensity. S', will not produce R because it involves a 
totally different group, s'i of synapses. Each intensity 5 of 
the same stimulus has thus a representative individual group 

A, A 


Fig. 39 

of synapses in the nerve centers, and therefore each intensity 
may, in a way, be considered as a dif event stimulus-pattern. 

However, things will happen in this fashion only if the 
difference between the intensities S and S' is sufficiently 
large. If this is not the case, then hm and h'^ differ very little, 
and the corresponding intervals {h^, hi) and {h^, h'^ will 
partially overlap (Fig. 39). If, now, 5 produces reaction R, 
then S' also produces R through aU the synapses which cor- 
respond to the interval (Fig. 39). Since, because of 

equation (3), hx and h, as well as hm, are functions of Ex, or, 
what amounts to the same, functions of S, the minimum 
difference AS — S — S' for which the intervals {hx, h) and 
(h'x, h'j) do not overlap at aU is itself a fimction of S. The 
determination of this fxmction AS = U(S) is, in principle, a 
simple problem, which may, however, involve some compli- 
cated algebra. The function U(S) is determined by F( 5 ) in 
equation (i) and by Nih). Several problems suggest them- 
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selves at this stage. For instance, we may ask the following 
question. How should F{S) and Nih) be chosen in order that 
AS = U(S) would have a prescribed form — say that of 
Fechner’s Law? 

When (hi, h) and (h'l, hi) partially overlap, we may in- 
vestigate the relative intensity of R, as produced by S', com- 
pared to that produced by S. In the simplest case this rela- 
tive intensity will be given by the ratio rj of all neurons 
lying in the interval (h'l, Jh) (Fig. 39) to those lying in the 
interval (k,, hi). That is. 



If we make more complicated assumptions about the possible 
interaction of the neurons which lie on the efferent side of the 
synapses j,-, we shall obtain expressions different from (4). 
This leads to another interesting group of problems. 

The case where E(h) has several maxima in the interval 
(o, PEz) is more complicated but is treated in a similar way. 

A similar but slightly different situation is obtained by 
considering the case where the inhibitory fibers are not col- 
laterals of the ne neurons but are themselves excited at the 
S3mapses Si (Fig. 40). This case has been discussed in more 
detail in MB, chapter xxii, and has been used by A. S. 
Householder to develop a quantitative theory of discrimina- 
tion along the lines outlined above. We shall now follow 
Householder’s presentation.^ 

The quantity E(h) in equation (3) is a function of h and of 
the parameter S. Let us denote it by 4 >( 5 , h). Assuming that 
the relation (i) is a linear one, we can bring equation (3), by 
a suitable choice of units, into the form 

^(S, h) = iS- h)f(h) . 


(s) 
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Since e is proportional to E, we may, by a proper choice of 
units, make e = $. Let each synapse be connected with 
each other synapse Sh by an inhibitory fiber. Then the net 
amount e — _/ at s%, to be denoted by a{S, h), wiU be obtained 
by subtracting from #(5, h) the amount of j produced by all 
inhibitory fibers leading from other synapses to sj,. More- 
over, a-{S, h) is the effective stimulus acting upon all inhibi- 
tory fibers leading from sk. To set up the equation giving a, 
we make the following assumptions concerning the inhibitory 


Fig. 40 

fibers: {a) they are aU similar; (J) for each fiber the amount 
ofy produced is a linear function of a; and (c) the threshold of 
each fiber is negligible. 

Then the total amount of j produced by these fibers at any 
synapse Sh is 

I{S) = X/(r(5, h)dh , 

where X is a constant of proportionality measuring the activ- 
ity of the inhibitory fibers. The integration is to be extended 
over all values of h for which <r > o. Hence, 

aiS, h) = $(5, h) - X/(i(5, h)dh . 

We may suppose f{h) to be at least continuous and to 
vanish at h = 0 . Then h) vanishes sX h = o and at 
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h = S, for any given 5. Hence, #(5, h) has at least one 
Tnayimum in the interval. If we suppose it to have only one 
maximum, then for any 5 the graph of # will have just two 
points of ordinate I — say at A, and at — and for values of h 
between these two values c will be positive. Call this interval 
over which tr > o the “excited interval.” 

Putting 

I{S) = X (\is, h)dh , ( 6 ) 

we find 


and hence/ 


(r(5, h) = $(5, h) - I{S) ; 


I{S)^ 


X I '#(5, }i)ih 

Jht 

I + Xijh ~ ^i) 


(7) 


For determining h{S) and }h{S) we have, then. 


#(5, h) = $(5, h) = 


X (^^{S,h)dh 

Jht 

I + \{1h — hi) ’ 


( 8 ) 


since at these points o- vanishes. 

We now define the Weber ratio d(S) by the equation 

h.(S) = hiiS + S5 ) . ( 9 ) 

This is equivalent to saying that discrimination between the 
intensities S and 5(i + 5) is possible (in a suitable per- 
centage of trials) when the excited intervals are just distinct. 

By simple considerations A. S. Householder proves' that for 
very large values of S equation ( 9 ) is inconsistent with the 
foregoing assumption. In other words, the above-described 
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mechanism does not work when 5 is excessively large. The 
limits of possible variations of S are approximately defined by 
the value h* corresponding to the maximum otf(h). Neuro- 
physiologically, there is nothing unlikely in that result, in- 
asmuch as any observed regularities hold with sufficient 
exactness only within more or less limited ranges of stimulus 
intensities. 

Under those conditions we need consider only the ascend- 
ing branch of /(A). Since nothing is known about that func- 
tion, we may make the simplest assumption, namely, that 
over a wide range of values of S, the function f(h) is, with 
sufficient approximation, linear, so that 

#( 5 , h) = hiS - h) . (ro) 


In that case the graph of $( 5 , h) is an inverted parabola 
with a maximum at S/2. Hence, h and A, are equidistant 
from S/2. Define the “relative interval,” *, by 


Sx = S — 2hi= 2ht — S . 

(ii) 

Thus, 


^{S,h/) = $( 5 , h) = 7 

(12) 


(13) 

Then equation ( 8 ) becomes 


2'hSx^ + — 3 == 0 . 

(14) 

If we set 


2 \S 

u = 

(x 5 ) 
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we obtain 

ux^ + a? — I = o . (i6) 

To calculate 5 we use equation (9). Let x and u be the 
values of x and u corresponding to the value 5 (i + 5 ) of the 
stimulus. Then we have 


u = «(i + 5 ) . (17) 

From equations (9) and (ii) and the definition of x we ob- 
tain, after simple transformations, 


_ S — X 
5 -f- I 


(18) 


Then by defining a new variable a by the equation 

5 = x + z, (19) 

we obtain 

uz^ — 2 (x + i)s — (* 4- 1)“ = o . (20) 


Eliminating a from equations (19) and (20), we shall find an 
equation connecting u, 5 , and x. Then eliminating x from 
that equation and from equation (16), we obtain a relation 
between 5 and u. Since u is connected to the intensity 5 of 
the stimulus by means of equation (15), being actually pro- 
portional to S, we thus arrive at a relation between the Weber 
ratio 5 and the stimulus intensity S. 

For further details and discussion of the foregoing equa- 
tions we must refer the reader to the original paper by A. S. 
Householder.' Figures 41-45, taken from Householder’s 
paper, show a comparison of the calculated and observed 
relations between 5 and 5 ' (or, what amounts to the same, 
between 5 and u). For convenience S is plotted not against « 






T 


5 


10 log,„ u 
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but against logic u. From equations (15), (16), and (20) it 
follows that only one parameter, namely, the quantity \ 
is involved in the final relation between 5 and S. 

If the relation (i) between Ei and S is not a linear one, 
then equation (5) is not strictly equivalent to equation (3). 
If F(S) in equation (i) is of the type discussed in chapter viii, 
and if the deviation from linearity is small, within the range 
of values of S used, then a proper correction applied to equa- 
tion (5) will result, as shown by Householder,^ in a slight up- 
turn of the S, logic w curve, as indicated on Figures 41 and 
42 by the broken line. As will be noticed, the experimental 
points also indicate such an upturn. 
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CHAPTER XI 


MATHEMATICAL BIOPHYSICS OF PSYCHO- 
PHYSICAL DISCRIMINATION 

Hitherto we have discussed the problem that arises when 
we ask whether two intensities which are only slightly differ- 
ent from each other are definitely perceived as different or 
not. This problem, however, has another aspect. The two 
intensities, or, more generally, the two stimuli, may be so 
close to each other that a definite discrimination is impossible. 
Yet, even then, there may be a certain probability for a cor- 
rect statement with regard to whether the two stimuli are 
equal or not. The problem acquires a still slightly different 
aspect when we ask not merely for a statement whether the 
stimuli are equal or different hut also for a statement as to 
which of the two stimuli is the greater (or the stronger) one. 
The closer the two stimuli are to each other with respect to 
their intensities the less will be the probability for a correct 
judgment. We may ask for the functional relation between 
the difference of the two stimuli and the probability of a cor- 
rect judgment. 

An interesting approach to this problem has been made 
by H. D. Landahk by considering the neurological structure 
shown in Figure 46. This structure is a particular instance of 
the structure considered in MB, chapter xxii. 

Let the stimulus (Fig. 46) elicit a reaction R, through a 
rbai-n of excitatory fibers connected by synapses Si and 53, 
while the stimulus S2 elicits the reaction R^ through synapses 
$2 and 54. However, let the synapses and ^2 excite also 
inhibitory fibers which are arranged as shown on Figure 46. 
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We have essentially the well-known scheme of reciprocal 
innervation. 

For simplicity, let us consider the case where the excita- 
tory fibers are of purely excitatory type (chap, viii, p. 124), 
while the inhibitory ones are of a purely inhibitory type. 


R, 



4 

Fig. 46 


That is, the former produce only the e-factor; the latter, only 
thej-factor* This restriction is not an essential one. 

The stimuli Sx and S2 result in the production, at the 
synapse .^3, of an amount €3 of e by fiber III and an amount 
7*4 oij by the fiber /F'. Similarly, at the synapse .^4 we have 
the amounts €4 and J3 when the stimuli Si and S2 are applied. 
If the constants of all the fibers are the same, and in particu- 
lar when A — B and a — b, then we always have €3 = 7*3 and 
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64 = When Si = Si, then at the synapses and ^4 we 
have, correspondingly, 

“ *4 i ji ~ Ja j 

and hence, 

«3 -^4 = £ 4 -/s = o- ( 1 ) 

In other words, when the two stimuli Sj and Si are equal, 
they mutually inhibit each other, and neither Ri nor Ri is 
produced. 

If, however, one of the stimuli— for instance, 5 — is much 
larger than the other, S3, then the amount 63 is increased, 
resulting in an increase of e, — At the same time,y3 is also 
increased, resulting in a decrease of 64 —jy Therefore, when 
€3 — ^4 becomes so large as to exceed the threshold h! of fiber 
V, Ri wiU be elicited, but Ri will still be inhibited because 
€4 — jj < o. For a perfectly S3nnmetrical scheme, like the one 
shown in Figure 46, we have, quite generally, 

^a-3i = -i^i-jA)- ( 2 ) 

Therefore, at only one of the synapses or ^4 can e — y be 
positive. Hence, the simultaneous presentation of 5 , and Si, 
regardless of their intensities, results either in no reaction at 
all {Si = 5 j) or only in one reaction (i?i if 5,^52; Ri if 
Si«Si). 

Instead of speaking of the stimuli and 52, we may speak 
of the amounts e, and 62 of « at the synapses Si and ^2, since 
these amounts are monotonically increasing functions of 5 , 
and 53, respectively. In order to produce the response Ri, 
the quantity €3 — at the s}mapse jj must exceed the 
threshold h'. But in order that this should happen, we must 
have a sufficient excess of Ci over €2; in other words, e, — e, 
must exceed a threshold value h, so that 

Cx €3 ^ ^ • 


(3) 
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The quantity h is not to be confused with the threshold h' of 
the fibers III, IV, V, and VI or with the threshold of the 
fibers I and 11. The quantity A is a function of V , but it in- 
volves other paranaeters also. 

This situation holds only as long as all neurophysiological 
processes are perfectly regular and constant, not being sub- 
ject to any external or internal disturbances. However, as we 
have discussed in M.B, chapter xxviii, in general we must 
expect spontaneous fluctuation of excitation to occur in the 
central nervous system. Such fluctuations may be due to 
fluctuations of metabolic activity or to excitation carried to a 
given region from a number of other regions of the brain, 
which are randomly excited by the stream of oncoming 
exteroceptive as well as proprioceptive and enteroceptive 
stimuli. 

Let us consider such fluctuations at the synapses Xj and Xa. 
They result in an addition of a varying amount of < — y to 
either €i or Ca- In the absence of fluctuations, whenever we 
have €i — €3 > h, the reaction is produced. Since an in- 
crease of ei is equivalent to a reduction of € 2 , as far as the re- 
lease of Ri is concerned, we may consider fluctuations of «t at 
the synapse Xj only. 

In the presence of such fluctuations the following will gen- 
erally occur. Suppose Sj ^ ^aj then €i )£> Ca- In the absence 
of fluctuations this will result in the reaction Ri. But if, owing 
to the fluctuations, an amount (e — j)' < o is added to 
at the synapse Xi, then if the absolute value of this amount 
(e — j)' is sufficiently large, the net e — y at the synapse Xi 
will become less than h; and, therefore, either no reaction at 
all will be produced or a reaction Rt will occur because of a 
decrease of y'j at the synapse x^, in spite of the fact that Sl> S^- 
Let us use the designation “correct response” for the response 
J?i when Si > Si, or for Ri when Si < S,', and the term “wrong 
response” for the reaction Ri when < 5a or for R, when 
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Si > S2. With this terminology we see that the spontaneous 
fluctuations of excitation may result in wrong responses. The 
probability of a large fluctuation {e — jj is smaller than the 
probability of a small fluctuation. Since, when Si > Si, the 
larger the Si — Si the larger must be the absolute value of 
the additional (e —j)' < o in order to produce the wrong re- 
sponse; therefore the probability of a wrong response increases 
when the difference Si — Si decreases. 

Quantitatively we can determine the relation between the 
probability of a given response and the difference of the 
stimuli in the following way.*’* From the preceding discussion 
it follows that, in order to produce a reaction Ri, we must 
have (cf. Eq. [3]) 


[«i + (€ - j)'] - ei> h. (4) 

The wrong reaction is produced when the left-hand side of (4) 
is less than — h. When the left-hand side of (4) lies between 
— h and + h, no response is produced. 

Hence, 

I. If (e — j)' > — (€i — €2 — h), the correct response is 
made; 

II. If {e — j)' < — («i — €2 H- h), the wrong response is 
made; and 

III. If — (€i — €2 — A) ^ — j)' ^ — (ti — £2 -f h), no 

response is made at all. In other words, things occur as if 
Si = Si. We shall call this last case “equahty response.” 

Let 

f{e-jy die-jy 

denote the probability of having the fluctuation lie between 
(e - j)' and (e - j)' -f d{e - j)'. We have 

J^y{e-jydie-jy = j. 


is) 
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Since the correct response is made whenever case I 
is satisfied, the probability of the correct response is 
equal to the probability of having (e — jY lie in the interval 
[— (€i “ €2 “ ^), + 00], a probability given by 




p(x)dx ; x = (e- jY 


Similarly, we have for the probability of the wrong re- 
sponse and for the probability P, of the equality response 
from cases II and III, respectively, 


p(x)dx 

Pe = I p(x)dx . 

(«i— €a4”^) 


pix)dx 


If the function p{x) is given, then equations (6), (7), and (8) 
give us Pc, Pu, and Pc in terms of «2, h, and of any param- 
eters of the function p(x). On the other hand, the quanti- 
ties €i and 62 can be expressed in terms of intensities of con- 
stant stimuli Si and Si and the duration, t, of their applica- 
tion. In fact, denoting the intensities of excitation of fibers 
I and II (Fig. 46) by Ei and Ej, respectively, we have (chap, 
viii) 

ca = (i — er^‘) . 

a 

Using now, together with equations (9), the relations (2) and 
(3) of chapter viii, which gives, approximately, 


El = K log ■ 


& = .S' log ' 


K = lah 
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where K is the threshold of fibers I and II, we find relations 
between €t and St and 62 and S2. Thus, equations (6), (7), 
and (8) express the probabilities P*, P„, and P„ in terms of 
St, St, the time t and the different parameters entering into 
the function p{x). In the experiment to be described later, 
to which some of the present considerations are applied, 
the time i is kept constant. 

If we assume for p(x) the normal distribution function^ 

p(x) = , (ii) 

V 27r<r 

then equations (6), (7), and (8) express P., P„, and P, in terms 
of St, S2, and the two parameters h and <t. In this case, if we 
take two observed values of either P„ P„, or P* for any two 
pairs of St and S2, we can calculate h and o- and then calculate 
the values of P„, P,„, and P* for any other pair St, S2. Since 
in the case of a normal distribution (ii) the integrations in 
equations (6), (7), and (8) cannot be made in closed form, 
numerical tables of the probability integral have to be used. 

Other distribution functions, such as, for instance, 

pix) = j (j2) 

may also be considered. In equation (12) Vi is a constant and 
I X I denotes the absolute value of x. In this case the integrals 
in equations (6), (7), and (8) can be calculated in closed 
form. A comparison of the functions (ii) and (12) is shown 
in Figure 47. 

These theoretical results have been applied by H. D. 
LandahP to the cases of discrimination of weights. Table i 
shows the experimental data. This table means, for instance, 
that when a standard weight of 200 gm ( 5 t) is compared with 
another one of 195 gm (So), in 15 per cent of the judgments 
by a subject the actually greater weight was estimated as 
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being smaller, in 25 per cent both weights were estimated as 
equal, and in 60 per cent the correct estimate was made. The 
second lines, corresponding to the weights in the first column, 
are obtained by modifying the experiment somewhat, namely, 
by aUowiag only “greater-than” or “less-than” judgments, 
excluding the equality judgments. In this case, for instance, 
when a comparison of the 195-gm weight with the standard 



Fig. 47. — ^The full line represents the normal distribution function, 
equation (ii). The broken line represents the function given by equation 
(12). The areas of the two curves on this figure are equal. 

200-gm weight is made, it is found that there are 20 per cent 
wrong judgments and 80 per cent correct ones. 

The theoretical argument leading to this second mode of 
experimenting needs a little more discussion. Whenever con- 
dition III is satisfied, neither Ri nor R2 are produced — in 
other words, the subject does not state either stimulus to be 
greater than the other. If, however, the subject is instructed 
not to make any equality judgments but always to decide on 
a “greater-than” or “less-than” judgment, this will result in 
an effort to sharpen his own acuity for the case, when he other- 
wise would pronounce an equality judgment. Such an effort 
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may be interpreted biopbysically as a general increase in 
excitability, resulting m a lowering of the threshold h. When 


TABLE 1 

Experimental Data 


s 

(Grams) 

Pw 

Pe 


i8s 

0.05 

0.04 

0.91 

•OS 


( - 95 ) 

190 

. 12 

.18 

.70 


.15 


•85 

IQS 

•15 

•25 

.60 


.20 


.80 

200 

•30 

( -42) 

.28 


.52 


.48 

205 

■55 

•35 

.10 


.85 


.15 

210 

.70 

.18 

.12 


■85 


■IS 

215 

.85 

0.09 

.06 


0.93 


0.07 


TABLE 2 


Theoretical Values 


s 

(Grams) 

Pw 

Pe 

Pc 

i8s 

0.03 

0.06 

0.91 


.05 


•95 

190 

.07 

■13 

.80 


.11 


.89 

IQS 

.14 

.27 

•59 


.24 


.76 

200 

.29 

.42 

.29 


■50 

i 

1 

•SO 

20$ 

.58 

.28 

.14 


.76 


.24 

210 

.79 

.14 

.07 


.88 


.12 

215 

.89 

0.07 

.04 


0.94 


0.06 


h — o, then from equation (8) we have = o. For o, 
we find from equations (5), (6), (7), and (8) 

P, + P^ + Pe= = I • (13) 


For A = o, we have 


4 " -Pu) = ^ } 


(14) 


as we should expect from physical considerations. 
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In using data for “two-category” judgments only, we put 
A = o in our equations. In this way we can calculate <r from 
any one pair S^, Sa- 



Fig. 48. — Graphs showing the comparison between the experimental 
data in Table i and the theoretical data in Table 2, for probabilities of 
correct, wrong, and equality judgments, A perfect agreement would bring 
all points onto the straight line. From H. D. Landahl.* 

The values in parentheses in Table i have thus been used 
for computing h and crj, assuming for p{x) the distribution 
function (12). With the values thus found, all the remaining 
probabilities for other values Sx, Si were calculated; Table 2 
shows these calculated values. A comparison of the two 
tables is best made graphically and is shown in Figure 48, in 
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which the theoretically calculated values are plotted against 
the experimental values. 

From equations (6), (7), and (8) it is seen that P„ and 
Pf are functions of the difference Cj — €2. For a constant 
value t of the duration of stimuli, such as is used experi- 
mentally, we see from equation (9) and (10) that therefore 
P(,, Pw, and P, are functions of the ratio Si/Si. 


.JO 



Fig. 49. — The curve represents the theoretical relation between the 
probability of a correct or wrong judgment and that of an equality judg- 
ment (cf. Eq. [14]). The circles and crosses represent experimental data 
by J. P. Guilford.'* From H. D. Landahl.' 

Hence, for p{x), given by equation (12), 

P«=/,(|,^,<rx) ; P, = /.(|,A,<Tx). C15) 

Since for a given experimental setup h and oi are constants, 
therefore elimination of Sx/St from the two equations (14) 
gives us a relation between Pj and Pj. Similarly a relation is 
obtained between P„ and P«. The theoretical relation ob- 
tained by numerical evaluation of the integrals involved is 
compared with experimental data in Figure 49. 

For a different set of data on weight discrimination, similar 
comparison was made by using a normal distribution func- 
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tion, represented by equation (ii). The results of the com- 
parison of the theory with experimental data are shown in 
Tables 3 and 4 and in Figures 50 and 51. 

One complicating circumstance must, however, be dis- 
cussed in this case. The fimctions (ir) and (12) are 55011- 


TABLE 3 TABLE 4 

EXPERIMENTAL DATA THEORETICAL VALUES 


5 

(Grams) 

Fv , 

Fe 

p S 

® (Grams) 

Pvi 

Pe 

Pc 

84 

0.0X2 

.020 

0.027 

0.961 84 

,980 

0.0039 

.0103 

0,0207 

0.9754 

•9897 

88 

.021 

( -053) 

.082 

.897 88 

•947 

.0248 

( .0530) 

.0774 

.8978 
• 9470 

92 

.096 

.185 

( .181) 

.723 92 

.815 

.1027 

.1791 

( .1811) 

.7162 

.8209 

96 

.275 

.420 

.266 

•4S9 96 

.580 

.2845 

.4092 

.2646 

•4509 
■ 5908 

TOO 

-502 
( -683) 

.267 

.231 100 

•317 

-5512 
( .6830) 

.2502 

.1986 

.3170 

104 

.842 

.920 

.103 

.055 104 

.080 

• 79S7 

• 8797 

.1401 

.0642 

.1203 

108 

•915 

0.963 

0.065 

.020 108 

0.037 

.9320 

0.9663 

0.0545 

•0135 

0.0337 


metric with respect to ic = o. This means that a positive 
fluctuation is as likely to occur as a negative one. This occurs 
physically when there is no bias of any sort in favor of one 
of the stimuli. Frequently, however, a bias is introduced by 
the experimental setup. In the experiments of J. P. Guilford,^ 
represented in Table i, the two weights to be compared were 
presented in different orders several times. In the experi- 
ments of F. M. Urban,s represented in Table 3, the standard 
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I7S 

weight was always presented first. This results, in the case of 
Urban, in a bias of the subject in favor of the second weight, 
which appears to be larger. Such a bias has the same effect as 



Experiment — ^ 

Fig. 50. — Same as Fig. 48 but for a different set of data, as given in 
Tables 3 and 4. Computed by H. D. Landahl. 



Fig. 51. — Same as Fig. 49 but for a different set of data. Computed 
byH.D. Landahl. 
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shifting the distribution function to the right or to the left 
along the axis of the abscissae. This gives, instead of (ii), 

^(x) = - ^ ( j 6 ) 

V 2X0- 

where Xq is another parameter. A corresponding modification 
is introduced into function (12). Since now the distribution 
function is a two-parametric one, we have altogether three 
parameters: A, cr, and Xq, Therefore, three values from the 
experimental Table 3 are now used in order to calculate the 
remaining ones. 

Another difference to be noted between the data of J. P. 
Guilford and those of F. M. Urban is that, in the former, all 
observations were made on one subject, while, in the latter, 
seven subjects were used, so that Table 3 is obtained by- 
averaging for all subjects. 

Analysis of visual and auditory data have also been made 
by H. D. Landahl.^ He has also discussed more complex cases 
of more than two stimuli, as well as a relation between the 
probability of a correct response, the difference between the 
stimuli, and the time allowed for comparison of the stimuli.^ 
This leads us into the interesting problem of the relation be- 
tween the performance, difficulty of a task, and time allowed. 
For details we must refer the reader to the original papers. 
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CHAPTER XII 


PERCEPTION OF VISUAL PATTERNS 

As we have mentioned in chapter viii, when a constant 
stimulus is suddenly applied to a fiber characterized by 

BKA; Ka-, (,) 

then at first, for a short time, « exceeds y; finally, however,; 
again exceeds e and remains larger than e as long as the stimu- 
lus lasts. This is due to the fact that, because of A/a < B/h, 
the asymptotic value AE/a of e is less than the asymptotic 
value BE/h of j. Yet, because of R < A and i <a, the 
quantity t increases more rapidly at first, and also approaches 
its asymptotic value more rapidly than the quantity j. If, 
while the first stimulus So lasts and; is therefore still larger 
than €, we again apply, suddenly, an additional stimulus 5j, 
so that the total stimulus now becomes 5o + ■S'l, then, as can 
be readily seen from Figure 52 , again for a short time there 
will be excitation, followed again by a lasting inhibition {MB, 
chap. xxvi). It is, however, necessary that the additional 
stimulus be applied sufficiently suddenly {MB, chap, xxvi) 
and that it exceeds a certain threshold. 

When such a fiber forms a synapse with another fiber, then 
the latter becomes excited for a short time only when the 
intensity of excitation of the first fiber suddenly increases. If 
the first fiber is excited continuously with a constant intensity 
E, the second fiber is not excited. As we have seen in MB, 
the intensity of excitation of the second fiber increases with 
increasing A£ where AE is the amount of sudden increment 
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of E. As can be readily seen from the discussion chap- 

ter xxvi, the liminal ZVE which still produces an excitation of 
the second fiber is a linearly increasing function of E. Since 
E is, in general, a function E = J{S) of the intensity of the 
peripheral stimulus, which involves the first fiber, therefore 
this liminal AE will, in general, be a more or less complicated 
function of S, depending on the shape of f{S). 



If the second fiber leads to a region Ci, discussed in chapter 
X, then to each intensity of excitation of the second fiber 
there will correspond a definite group of excited neurons 
And, since a definite intensity E* of excitation of the second 
fiber corresponds to a definite relative sudden increment 
AE/E of the intensity of excitation of the fiarst fiber, while this 
relative increment itself corresponds to a definite relative 
increment of the peripheral stimulus, therefore, to every 
value of AS/ S there will correspond a definite individual 
group of neurons ne^ As has been shown in MB, chapter 
xxvi, a similar mechanism is obtained for a sudden decrease of 
E or of 5 — that is, for a negative iS/S. 
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We thus have a situation in which a sudden change of in- 
tensity of a constant external stimulus results in a short 
excitation of a definite group of central neurons, different 
groups corresponding to different amounts of change. 

Combining this result with that of chapter x, we may say 
that to each absolute value of intensity of a peripheral stimu- 
lus there is a corresponding excitation of a definite center 
consisting of a definite group of individual neurons ne ^ ; and to 
each value of relative change of the intensity of peripheral 
excitation there is also a corresponding excitation of a similar 
definite center. The finite, though very large number of, 
neurons in the brain causes an overlapping of some of these 
centers and thus causes the perception of continuity of the 
possible values of intensity of peripheral excitation. 

We shall now introduce some neurophysiological assump- 
tions which may appear rather questionable. We are intro- 
ducing them here, not because we consider them to be any 
more likely than other possible assumptions, but in a tenta- 
tive way, as a sort of working hypothesis, in order to show 
how, by means of such or similar assumptions, we may de- 
velop a physicomathematical theory of some phenomena of 
visual perception, the discussion of which hitherto frequently 
remained on a purely qualitative level. Later on we shall dis- 
cuss some possible modifications of these hypotheses — ^modi- 
fications which, however, may leave the fundamental formal 
results unaltered. 

When we look at a segment of a straight line, we succes- 
sively fix our attention on different points of this line. This 
results in movement of our eyes along that segment, back and 
forth. Experiments on eye movements^ show that there is no 
simple relation between the shape of a figure contemplated 
and the path of the eye movement. Actual following of a 
rectangle gives a rather irregular path for the eye movement.® 
However, inasmuch as the eye muscles are innervated by 
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several centers in a rather complicated fashion, the following 
hypothesis will not necessarily be at variance with observa- 
tions. Suppose that there is a group of brain centers which 
innervate the eye muscles in such a way as to make the eye 
follow exactly a contour which is projected on the retina. 
Each eye muscle is, however, mnervated by a number of other 
centers, which produce different movements, superimposed 
on the above-mentioned movements. Thus, the actual move- 



Fig. S3 


ment of the eye may have no apparent relation to the shape of 
the contour contemplated; yet it contains as a part those 
movements which trace an exact replica of the geometrical 
pattern looked at. 

In Figure 53, Cs represents schematically the center which 
would produce such exact movements of the eye, through the 
muscle M, while Cm stands for all other centers which may 
impart to M movements that have nothing to do with the 
shape of the contour. The character of the movement which 
would be produced by M under the influence of Cs alone de- 
pends on the character of excitation of Cs (intensity, dura- 
tion). If through a synapse 5 a fiber a is excited which leads 
to a sensory center Cp (Fig. 53), then to each definite type of 
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excitation of Cs, or, in other words, to each definite t3rpe of 
movement produced by M under the influence of Cs alone, 
there corresponds a definite intensity and duration of excita- 
tion of Cp. li M were to be innervated by Cs only, then the 
same type of excitation in Cp could be obtained through pro- 
prioceptive fibers, which would lead from M to Cp and which 
would be excited by the contraction of M. In such a case the 
proprioceptive impulses would have a definite relation to the 
shape of the contemplated pattern. Actually, because of the 
other centers Cm, the proprioceptive impulses coming from 
M have no relation to the shape of the pattern. But the im- 
pulses coming through s over a (Fig. 53) do have such a rela- 
tion. In other words, according to the scheme represented by 
Figure 53, things will occur at Cp as if the movements pro- 
duced by M would exactly correspond to the shape of the con- 
tour contemplated and as if these movements would send off 
proprioceptive impulses directly to Cp. 

The considerations above therefore enable us to speak, for 
simplicity, of the eye as actually following the contemplated 
contour, and of proprioceptive impulses corresponding to such 
movements being sent to some sensory centers. Although 
actually things are much more complex, by means of this 
scheme we may interpret our schematized statements m 
terms of a more exact neurophysiological picture. Thus, when 
in the following we say; “The contraction of the eye muscle 
sets up proprioceptive impulses which behave so and so,” we 
do not mean it literally, but actually imagine a scheme like 
the one represented by Figure 53. 

Neurological mechanisms that would produce such a move- 
ment can be readily suggested. For the time being, let us 
make, for simplicity and as a special, purely theoretical case, 
the assumption that this movement of the eyes goes on with 
constant velocity v. Let a line segment AB be placed in a 
vertical plane, and let its angle with the vertical be 6 (Fig. 
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54). Neglecting the possible role of the oblique eye muscles, 
the movement of the eyes in the direction AB involves a con- 
traction of the rectus superior of both eyes — the right rectus 
lateralis externus and the left rectus lateralis internus — and 
a relaxation of both recti inferiors — the right 

V lateralis internus and left lateralis externus/ 
Again assuming for simplicity that the co- 
ordinated movements of the left lateralis in- 
ternus and the right lateralis externus are 

A controlled from one common center, and that 
the same holds for the right lateralis internus 
and the left lateralis externus, we may greatly 
A o A' simplify our considerations by considering 
Fig. 54 movements of one eye — for instance, 

the right one. Then the movement in the 
direction AB (Fig. 54) involves a contraction of rectus su- 
perior and of lateralis externus and a relaxation of rectus in- 
ferior and lateralis internus which may be written in an ab- 
breviated form, thus: 

AB:rs A- Ze + li— (2) 

The movement in the direction BA involves 

BA :rs — le — ri + li + (3) 

The changes of muscular contraction of the four muscles pro- 
duce proprioceptive impulses of different intensities; and if 
the proprioceptive fibers lead to a region Ci, discussed in 
chapter x, then any given intensity of contraction will produce 
excitation in a definite individual group of neurons ne^ 
Again, for simplicity, we shall consider that only positive con- 
tractions produce proprioceptive impulses, a relaxation being 
characterized by a lack of such impulses. Let the intensity 
of these impulses be proportional to the velocity of the con- 
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traction of the muscles. Then, if v is the velocity of the eye 
movement, the velocities of contraction of rs and /e in (2) are, 
respectively, 

V cos d and v smd ; 

and the corresponding intensities of the proprioceptive im- 
pulses are 

£„ = av cos 0 ] Eu = av sind , (4) 

where a is a coefficient of proportionality. 

The excitation excites a definite individual group of 
neurons ne^ in a center which we shall call V (for vertical) ; 
while Eie excites a group of neurons in a center (horizontal). 

The movement of the eye in the direction BA produces the 
following excitation: 

Eri = av cos 0 ] Eii = avsm0 , (s) 

and causes the excitation of two other groups of neurons 
ne^, one in V and one in E. If the eye moves up and down 
indefinitely, four groups of neurons are excited. If we now 
consider another segment of a straight line, which is char- 
acterized by a different 0 — 0’, again four groups of neurons 
ne^ will be excited; but these four groups will be different from 
the former four, because E„, Eu, Eu, and Eh are different. 
Some of the new groups may partly overlap the old ones only 
when 0 — 0' = A0 is very small. Thus, every position of a 
straight-line segment in a vertical plane corresponds to the 
excitation of four distinct individual centers, two in H and 
two in V. 

Considering, now, the segment A'B' (Fig. .54), symmetrical 
to AB with respect to the vertical 00', we shall find, by a simi- 
lar argument, that the scannmg of that segment with the 
eye in the direction A'B' gives 

£ ', = av cos 0 = Er, ; E'h = ot sin 0 = Eu , (6) 



i84 mathematical BIOLOGY 

and in the direction gives 

EU = cos d - Eri] E'le == av sin 6 = Eu, (7) 

That is, AB and jL'B', scanned in both directions, will excite 
identical groups of neurons ne3 in H and V. This difSculty 
can be avoided in several ways, only one of which we shall 
consider here. 

Let, as suggested by Gale Young, the neurons of higher 
order, corresponding to li, le, and n, have such high thresh- 
olds that they remain unexcited no matter how strong the 
excitation of the peripheral fibers. However, let those higher- 
order neurons also be excited by collaterals from rs^ so that 
while, for instance, excitation of li alone does not result in 
any central excitation, yet simultaneous excitation of rs and 
li results in such a central excitation. 

In that case a contraction of the muscle rs, that is, a move- 
ment of the eye upward, if accompanied by a lateral move- 
ment, produces an excitation both in V and H. But a down- 
ward movement of the eye, produced by a contraction of ri 
and a relaxation of rs, does not result in any excitation in 
either V or H. Of course, other brain centers may be excited 
by such movement by means of other fiber connections. But 
at present we are interested only in V and H. We now have, 
for the movement along AB (Fig. 54), 

ABirs “f“ ; le A' \ ri — ; li — * (8) 

and for BA: 

BA :rs — ; le ; ri *— ; li ~ . (9) 

For A'B^ we have 

A'B^irs + ; W + ; ri — ; le (10) 



PERCEPTION OF VISUAL PATTERNS 185 
and 

B'A'irs - ; U - ; ri - h - . (11) 

From (8) and (9) it foUows that scanning of AB in both direc- 
tions results in 


Ers = av cos 6 ; Eu = av sin 8 , (12) 

while scanning A'B' in both directions results in 

E'rs = av cos 8 ; EU = av sin 8 . (13) 

Both AB and A'B' involve rs in the same amount and there- 
fore excite the same group of neurons iu V. But, while AB 
involves le, A'B' involves li. Hence, the two groups on 
neurons m the H-center, which correspond to the s cannin g 
of the two lines in both directions, wiU be different. Every 
segment of the straight line excites one group of neurons in E 
and one in V, the groups differing for different 6 ’s. Lines 
symmetrical with respect to either a vertical or a horizontal 
line have a common group in V. Parallel segments produce 
identical excitation. 

In order to obtain excitation in H center when a horizon- 
tal line is scanned, not involving any upward movement of 
the eye, we must introduce an additional assumption, name- 
ly, that only an actual downward movement of the eye re- 
sults in a relaxation of rs, while holding the eye on the same 
level still requires a tonic contraction of rs. 

The proprioceptive excitations, owing to the contempla- 
tion of any segment of a straight line, may thus, through 7 
and E, be conditioned to a response R. A different segment 
will, in general, not produce R, because it involves different 
groups of neurons in 7 and E. A segment S)anmetrical to the 
original one with respect to a vertical line may, however, pro- 
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duce a weaker R through the group of neurons in V, which it 
has in common with the original segment. 

Now consider two segments, AB and BC (Fig. 55), form- 
ing an angle at B. Following this figure with the eye results 
in the excitation of two groups of neurons in V and of two in 
H. If, however, the proprioceptive fibers, leading from rs, li, 
and le to V and E, send off collaterals, which are character- 
ized by relations (i) and each of which leads to a center A*, 
A* will become excited every time the eye passes the angle B, 
because of a sudden change in the intensity of excitation E„ 
and Eie (or of Bra and Eh for a differently oriented angle). 



Fig. sS 

To each value of the angle B there is a corresponding dif- 
ferent intensity of excitation in A*; and if A* is connected to 
a center A, of the kind discussed in chapter x, then to each 
value of the angle B there corresponds a definite group of 
neurons ne^ in A. We thus have a center for the perception of 
angles. 

Making different hypotheses concerning the physical 
characteristics and constants of all centers involved, we may 
derive various quantitative relations which may suggest ex- 
perimental studies. As an iEustration, an expression has been 
derived for the relation between the minimum perceptible 
angle formed by the two straight lines AB and BO and the 
angle which one of them — say AB — ^forms with the vertical.'’ 

Finally, we may consider not only that the speeds of the 
contraction of rs, li, and le control the excitation of some cen- 
ters but that special fibers are provided in which the inten- 
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sity of excitation is proportional to the amplitude of the 
muscular movements. In this fashion, by considerations 
similar to those above, we shall arrive at a picture in which 
a particular group of neurons ne^ is excited in a center L, the 
group being characteristic for the length of the segment of 
straight line considered. 

If a polygonal contour, consisting of n sides and n angles, 
is presented to the subject, then it follows from the foregoing 
that, as this contour is followed by the eye, in general, n 
groups of neurons will be excited in F, « groups m E, n 
groups in A, and n groups in L. The intensities of excitation 
of the different groups, even belonging to the same center, 
such as H or F, will generally be different. Two groups of 
neurons corresponding to two different distinct intensi- 
ties E, will not only be distinct but will, in general, contain 
different numbers of neurons. If the polygon possesses some 
symmetry properties, which result, for instance, in some of 
the angles being equal, this reduces the number of distinct 
groups in ^ . If there are m equal sides, then they all excite 
the same group of neurons ne^ in L, etc. As the contour is 
followed by the eye, the excitation of each group comes and 
goes, discontinuously, occurring only when the attention is 
concentrated on the particular element (segment, angle, or 
length). If, out of n elements, m are identical, the group ne^ 
corresponding to the m identical elements k, will be excited, 
on the average per unit time, m times more often than those 
groups we, which correspond to elements that have no other 
identical ones. Hence, the average intensity of excitation of the 
group corresponding to m identical elements k will be m times 
the excitation due to a single element h. 

We may now consider a scheme, discussed in MB, chapter 
xxii, and in chapter x of this present work, in which an excita- 
tory fiber leads from neurons of each center E, V, A, and 
L to neurons of higher corresponding centers, E', V', A', 
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and L'. This excitatory fiber branches off into inhibitory 
fibers, leading to the neurons of all other centers. We may 
ask for the total excitation corresponding to a given polyg- 
onal contour. If this total excitation is transmitted to a cen- 
ter whose excitation results in a sensation of pleasure {MB, 
chap, xxx), then its intensity may be considered as a measure 
of the pleasantness or of the aesthetic value of a given con- 
tour. The problem is perfectly definite and can be treated 

1 iO 

□ o 


so 60 

h 0 

Fig. 56. — These ten polygons, selected from one hundred given in 
G. Birkhoff’s book,*' were used in the experiments of jR. C. Davis' to de- 
termine the relative aesthetic values of different geometric patterns. 

provided the intensities of excitation of all the groups we, 
in E, V, A, and L are given — in other words, provided the 
distribution function N{h), discussed in chapter x, page 149, 
for different centers, as well as the function given by equation 
(i) of chapter x, is specified. For different choices of these 
functions the problem wiU be of different degree of mathe- 
matical difBculty. A particularly simple case is obtained 
when we consider a very special case where the intervals 
(ki, hx) (cf. chap. X, Fig. 39) always contain the same number 
of neurons and where the intensity of excitation of any such 
group, falling in any interval (fh, h^, is approximately con- 
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stant. Such a case is rather unlikely to occur exactly, but we 
shall consider it here as a possible approximation. In that 
case two different elements — for instance, angles — are char- 
acterized by two distinct groups ne^, having, however, the 
same total intensity of excitation when the element is per- 
ceived. To m identical elements there corresponds one group 
with an average intensity mEt,, since the group is stimulated m 
times more frequently. 

We shall apply the foregoing considerations to a group of 
polygons for which G. Birkhoff'* has 
calculated quantitative aesthetic val- 
ues on the basis of some general con- 
siderations of complexity and order 
and for which those values have been 
measured by R. C. Davis by the 
rank-order method® (Fig. 56). Let 
the intensity Ea of excitation of a 
group ne^ be taken as o.i. The in- 
tensity of inliibition, which is approx- 
imately proportional to £0 (chaps, ix 
and xi) is taken to be o.oi per unit 
of Ea. Then for No. 50 of Figure 56, 
represented separately on Figure 57, we have the following: 

There are, altogether, eight segments of horizontal lines 
which do not give any excitation at V: 3-4, 5-6, 8-9, lo-ii, 
15-16, 17-18, 20-21, and 22-23; they aU cause excitation of 
one group ne^ in H, with an intensity of 0.8. The sides 1-2, 

6- 7, 13-14, and 18-19 are all parallel and therefore identical 
in their H and V groups. Hence, in the H center they give one 
group with Ea = 0.4. Similarly, 1-24, 19-20, 12-13, and 

7- 8 give, in the H center, one group with Ea = 0.4. In the 
V center the sides 2-3, 4-5, 9-10, 11-12, 14-15, 16-17, 21-22, 
and 23-24 give one group with Ea = 0.8; while 1-2, 6-7, 
7-8, 12-13, 13-14, 18-19, 19-20, and 24-1 give another group 
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in V with Eo = 0.8. In the A center we have one group with 
Eo = 1.2 corresponding to the right angles 3, 4, 5, 9, 10, ii, 
15, 16, 17, 21, 22, and 23; a group with Eo = 0.2 for i and 13; 
another group with Eo = 0.2 for 7 and 19; one group with 
Eo = 0.4 for 2, 12, 14, and 24; and one group with Eo = 0.4 
for 6, 8, 18, and 20. In the L center we have one group with 
Eo = 0.8 for the eight equal sides 2-3, 4-5, 9-10, 11-12, 
14-15, 16-17, 21-22, and 23-24; another group with Eo = 0.8 
for the eight equal horizontal sides; and a third group with 
Eo = 0.8 for the eight equal inclined sides. Altogether we 
have the following scheme: 


H I group of Eo = o . 8 

2 groups of Eo = 0.4 

V 2 groups of Eo = 0.8 

A I group of Eo = 1.2 

2 groups of Eo = 0.2 
2 groups of Eo = 0.4 
L 3 groups of Eo = 0.8 


To each group in the centers H, V, A, and L there is a cor- 
responding group in the centers H', V', A', and The first 
group of H' receives inhibitory fibers from all other groups, 
that is, from twelve groups altogether. The corresponding 
intensities of inhibition are 0.004, 0.004, 0.008, 0.008, 
0.012, 0.002, 0.002, 0.004, 0.004, 0.008, 0.008, and 0.008 — 
altogether 0.072. Hence, the net excitation of the first group 
mE' is 0.8 — 0.072 = 0.728. Each of the two other groups 
in E' again receives inhibitory fibers from twelve other groups, 
and the net excitation of the group is thus found, in a similar 
way, to be 0.324. Calculating, thus, the net excitation for 
each group and adding them all up, we obtain, for the total 
value of excitation, 7.04. In a similar way the total intensity 
of excitation of the groups were computed for the other ten 
polygons of Figure 56, and the results are plotted on Figure 58 
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(full light line). The broken and alternate lines of Figure 58 
represent the experimental rank-order values of R. C. Davis 
for two different groups of students. The dotted line repre- 
sents values calculated by Birkhoff. Inasmuch as our values 



Fig. 58. — The two full lines represent the theoretical values of the 
intensity of aesthetic sensation in arbitrary units for the ten polygons 
represented in Fig. $6. The broken and the alternate lines represent the 
experimental data by R. C. Davis.* The dotted line represents theoretical 
values calculated by G. Birkhoff.< See text, pp. 189 ff. 

give the actual intensities of central excitation, measured in 
arbitrary units, while Davis’ values give only the rank order, 
in which a certain arbitrariness enters in the assignments of 
weights, we cannot expect a close numerical agreement be- 
tween our values and those of Davis. However, the general 
trend agrees, with the exception of No. 20. 
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We notice that for No. 20 and No. 60 our theory gives 
values that are too low. The fact that No. 60 has actually a 
much higher value than No. 10, although both are septagons 
with approximately similar symmetry properties, suggests 
that the F-groups, which correspond to smaller e, are more 
strongly excited than those corresponding to larger d’s. The 
assumption of unequal excitation of the different groups 
would make our calculations rather cumbersome. Leaving 
this more rigorous study for the future, we shall approximate- 
ly take into account a possible preference for vertical direc- 
tions by multiplying the values of total intensity of excita- 
tion, as calculated, by the ratio of the heights to widths of 
the corresponding polygons. This is obviously a crude meth- 
od, holding only for ratios which are not too large. The re- 
sults are shown in Figure 58, heavy line. The general agree- 
ment is somewhat improved, at the expense of No. 80, which 
is too much lowered. In view of the very crude approxima- 
tions made, these results should command our attention. 
Further refinement and generalization of the theory, along 
lines indicated above, is likely to bring further improvements. 

If we do not make the rather crude assumption that all 
the centers are excited with equal intensity, then, obviously, 
the final results will depend on the assumption about the 
intensity of excitation of the different centers. These intensi- 
ties of excitation may vary from individual to individual, 
thus accounting for the differences in tastes (MB, chap, 
xxix). On the average, these intensities may be approximate- 
ly equal, which would account for the fair agreement of the 
present theory with observations made on large groups of 
people. 

The influence of various symmetries, that is, the appear- 
ance of several identical elements, may be easily analyzed by 
considerations similar to those given in chapter xxii of MB. 
If we have, for instance, » equally externally excited but 
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different elements, then each of the corresponding groups ne^ 
receives an inhibition from n — 1 other groups. If m of these 
elements are identical, then we have only n — m + 1 groups. 
The total applied excitation remains the same. In the first 
case it was nE; now it is (w — m)E + mE. But the group 
corresponding to the m identical elements receives inhibitory 
fibers from n — m other groups. Thus the total inhibition is 
decreased, and the net excitation increased. If we take, as a 
measure of complexity, the total number of elements in the 
pattern, then various types of sjntnmetry have the effect of 
reducing the “effective” complexity; but this reduction oc- 
curs by a process of subtraction rather than by a process of 
division, as in Birkhoff’s theory. It must also be remarked 
that, even when all elements are different, the total excitation 
in our case does not monotonically decrease with the complex- 
ity n. Equation (r) of chapter xxix of MB gives us the in- 
tensity of excitatory process at one of the n synapses, cor- 
responding to n different elements with equal intensity of 
excitation. The total intensity of excitation is obtained by 
multiplying that equation by n, which gives 

Etot = {K^S - K,)n - {K,S - KMn - i) • 

This expression is zero for » = o and has a maximum for 

_ (E. + K,)S - {K, + E4) 

” 2 {K^S - K,) 

after which it decreases to zero and becomes negative. We 
should therefore expect that in a rather simple pattern, con- 
sisting of a very few elements, too much symmetry will be 
rather a disadvantage from the aesthetic point of view, for 
in that case the total excitation increases with increasing com- 
plexity. But for very complex patterns, where the complexity 
is greater than the optimal »*, some symmetry is pleasant, as 
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reducing the complexity and bringing it closer to the optimal. 
In the general case of unequally excited elements, these rela- 
tions become much more complicated. 

If a subject or animal is conditioned to a certain polygonal 
figure and then he is presented with another figure, which has 
some elements in common with the first one, the intensity of 
the response to the second figure can be calculated by the 
method used in chapter xxix of MB. Thus, a quantitative 
theory of discrimination of patterns may be developed. Two 
patterns will be the more similar the more they have common 
groups ne^ in various centers. If aU the neurons in A send off 
fibers to a center A a, so that Ao is always excited, whenever 
any angle, no matter of what size is presented, then A o is a 
center that responds to “angularity” in general. An animal 
trained to choose a square, when a square and a circle are 
presented, will choose the more angular figure when two other, 
but different, figures are presented. This provides a basis 
for the understanding of some interesting results of H. 
Kluver^ on the equivalence or nonequivalence of certain 
pairs of geometrical figures in the behavior of monkeys. 

The mechanism discussed here suggests a definite pro- 
cedure for the generalization of the theory, including the 
aesthetic measure to curvilinear figures. If we again consider 
a fiber of the kind discussed on page 177, but if we consider, 
instead of a sudden change of E, a gradually increasing E, 
then we shall find, by using the fundamental equations de- 
veloped in MB, chapter xxii, that to any given variation of 
E with respect to time there will correspond a definite in- 
tensity of excitation of a definite duration. When a curve is 
followed by the eye with a definite velocity, then Er,, Eu, 
and Eu are not constants but are known functions of time, 
determined by the shape of the curve. Therefore, we may 
calculate the corresponding intensity of the excitation in a 
center C^, to which appropriate fibers are leading. We may 
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then calculate the minimum perceptible curvature and other 
interesting problems. We may also ask what particular shape of 
curve gives the maximum excitation at Cc (cf. MB, chap, xxii) . 

Experiments by Edmund Jacobson^ indicate that thinking 
of a vertical object, while not causing actual muscular con- 
tractions, sets up very weak action currents in the eye 
muscles, such as would correspond to a vertical movement of 
the eye. Thinking of horizontal objects gives currents cor- 
responding to horizontal movements. This suggests also that 
proprioceptive impulses of the kind discussed on page 180 
may be set forth without the muscles performing the actual 
movements. It has been observed that complete paralysis 
of the eye muscles does not affect recognition of shapes. This 
may be due to the taking-up of the functions of the eye 
muscles by other muscles of the body, which may imper- 
ceptibly move as a whole; or it may involve a mechanism of 
Gestalt recognition, discussed previously (MB, chap, xxvii). 

Finally, it may be shown* that a mechanism based on our 
fundamental equations developed in. MB can be conceived in 
which, when a chain of neurons is excited along a contour line 
in a brain center, a wave of excitation is set up, moving along 
that chain with constant velocity. The vertical and hori- 
zontal component of the velocity of that wave may then 
formally play the role of the vertical and horizontal move- 
ments of the eye. If such a mechanism is at work, then the 
central wave may continue to travel even for some time after 
the visual object is removed from the field of vision. 

The mam thing is that we have here a picture which pro- 
vides for the analysis of complex patterns and also gives a 
basis for a quantitative, exact theory. Most likely this 
mechanism is operating simultaneously with the one dis- 
cussed before (MB, chap, xxvii). The latter provides for 
recognition of gross features and perception of “wholeness”; 
the former, for analysis of details. 
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Linear algebraic equations, 13 
Linear change, 33 
Linear equation, 98 
Linear functions, 149, 153, 155 
Linear increasing function, 178 
Lines, 189 
segrnents of, 181 
straight, $8, 183, 186 
Liquid, 96 
Liquid medium, 95 
Liquids, 37 
Living cells, loi 
Loeb, J,, 102 
Logarithm, 75, 99 


Logarithm scale, 58 

Lorente de N6, R., 125, 129, 130, 

133 

Love, A. E. H., 68 
Lungs, 148 

Macdoneld, P. A., and Robertson, 
D. M,, 161 

Magnitude, orders of, 5, 7, 14, 28, 
43, 67, 96, 97, 106, 142 
Malignant, 80 
Mass, 78 
of a cell, 69 

Mathematical biology, 3 
Mathematical biophysics, 103, 127 
Mathematical Biophysics^ Bulletin 
of, X, 17, 36, 68, 81, 92, 98, 146, 
176, 196 

Mathematical neurology, 127 
Maxipia, 152, 154, 155 
Maximum excitation, 195 
Maxwell, J. C., 102 
Mean values, 81 
Measure of pleasantness, 188 
Mechanical effects of forces, 40 
Mechanical forces, 54, 55, 68 
Mechanically stable, 4 
Mechanism, 60, 178, 194, 195 
of division, 4 
physicochemical, 128 
Mechanisms, excitatory, 127 
Medium, 37, 38, 69, 78 
colloidal, 56 
external, 6, 8, 14 
internal, 14 
liquid, 95 
moving, 95 
resistance of, 3 

Membrane, 7, 8, 14, 17, 28, 37, 40, 
41, 78 
cell, 77, 104 

physical constitution of, 30 
resistance of the, 79 
surface of the, 30 
Membrane permeabilities, 7 
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Metabolic activity, loi, 166 
Metabolic flow, 77 
Metabolic forces, 88, 90 
Metabolic rate, 90, loi 
Metabolism, 2, 67, 92 
cellular, 92 
cessation of, 90 
effective rates of, 57 
respiratory, 75 

Metabolite, 14, 30, 40, 54, 56, 94 
respiratory, 78 

Metabolizing cell, 67 
Micella, colloidal, 56 
Minimum difference, 151 
Minimum perceptible angle, 186 
Minimum perceptible curvature, 

19s 

Mitosis, 60 
Mobility of ions, 104 
Molar rate, 55 
Molecular collision, 38 
Molecular structure, 37 
Molecular weight, 31, 38, 54, 55 
Molecule, 2, 34, 35, 38, 54, 56, 72, 
126, 128 
complex, 70 
free, 56 

Molecules, movement of, 2 
Mollendorf, W. von, 59 
Moment, dipole, 56 
Monkeys, behavior of, 194 
Monotonic, 193 

Monotonically increasing func- 
tions, 165 

Motion 
circulatory, 95 
steady, 95 

Motoneurons, 130 
Motor reactions, 148 
Movement 
of chromosomes, 60 
eye, 183, 184 
internal, 93 


of molecules, 2 
muscle, 187 

Movements, superimposed, 180 
Multicellular slice, 28 
Multiple response, 147 
Multiplication of organisms, i 
Muscles, 183, 189 
eye, 179, 180, r8i, 182, 184 
Muscular contraction, 182, 195 
Muscular movement, 187 
Mutual inhibition, 165 
Myelin sheath, resistance of, 118 

Natural sciences, i 
Neck, 51, 52 
“Necking,’’ 52 

Negative concentrations, 109 
Negative roots, X5 
Nemst, W., 103, 118 
Nerve, 102, 103, 105, 107, 124, 147 
peripheral, 129, 132 
Nerve cells, loi, 102 
Nerve centers, 151 
Nerve excitation, 103 
Nerve fiber, loi, 108, 118 
acoustic, 102 
afferent, 122 

electrical properties of the, 117 
Nerve fibers, bundle of, 147 
taste, 102 

Nerve-muscle preparation, 107 
Nerve trunk, 147 
Nervous impulse, 134 
conduction of the, 117 
Nervous system, central, 120, 126, 
127, 130, 166 
Neuroelements, 120, 123 
Neurological diagrams, 134 
Neurological mechanism, 148, 181 
Neurological structure, 163 
Neurology, mathematical, 127 
Neuron, 123, 126, 127, 129, 130, 
131, 135 ) 138, m , 140, 141, 
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Neuron — continued 

144, 149, 150, 152, 178, 179, 
182, 183, 184, 188, 194 
central, 179 
excitatory, 145 
of frog, 120 
inhibitory, 146 
of man, 120 

Neuronic structure, 139, 145 
Neurons 
chain of, 137 

geometrical arrangement of, 134 
group of, 1 8s, 186 
of higher order, 184 
interaction of central, 125 
Neurophysiological assumptions, 
179 

Neurophysiological picture, 181 
Neurophysiological processes, 166 
Newton’s third law, 2 
Niiella syncarpa, 99 
Nonequivalence of geometrical 
figures, 194 

Nonexcited regions, 117 
Nonlinear equations, 83 
Nonmetabolized substances, 78 
Nonspherical shape, 90 
Nonstationary case, 13, 15, 83 
Normal distribution function, r69, 

173 

Nucleus, 62, 93 

chemical and physiological role 
of the, 93 
dividing, 59, 60 

Oblong cell, 7, 12, 54, 58 
Oblong shape, 44, 88, 91 
Observations, 7, 18, 58, 59 
Offner, F., in, 119 
Oppenheimer, C., 99 
Optical phenomena, 102 
Optics, 102 
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and complexity, 189 
of magnitude, $, 7, 14, 28, 43, 
67 j 96, 97, 106, 108, 142 


Ordinates, 74, 154 
Organisms, 2 
multiplication of, i 

Organs, end-, 139, 141 
Origins, 141 

Osmotic pressure, 40, 78 
Output of lactic acid, 29, 34 
maximum, 34 

Over-all size, 7 
Oxidation, 30, 37 
of lactic acid, 30, 34 
rate of, 69 

Oxidized, 20, 25, 35 
Oxygen, i8, 26, 30, 32, 35, 37, 56 
concentration of, 22, 24-27, 34 
consumption of, 22, 24-27, 30, 

34, S4, 55; 58, 59 
rate of, 26, 58 
relative, 22, 23 
external supply of, 31 
Oxygen pressure, 30, 34 

Parabole, inverted, 155 
Parameters, 28, 29, 58, 88, 130, 
152, 162, 166, 168, 169, 176 
gas-kinetic, 126 

Partial differential equations, 13, 
103 

Particles, 38, 39, $6 
colloidal, 39 

Perception 
of angles, 186 
of continuity, 1 79 
of ^'wholeness,” 195 

Performance, 176 
Periodical contractions, 92 
Periodical solutions, 16 
Peripheral excitation, 131 
Peripheral fiber, 134, 147, 149 
excitation of, 184 
Peripheral nerves, 129, 132 
Peripheral stimulus, 121, 132, 136, 

151, 179 

Periphery, 5 
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Permeabilities, 21, 22 
membrane, 7 

Permeability, 7, 37, 48, 78, 82, 94, 
96, lOI 

experimental determination of, 
30 

infinite, 17 
Pharynx, 148 
Phase, 61 
Phosphate, 31 

Physical characteristics, 186 
Physical constants, 3, 14, 13 1 
Physical factors, loi, 102 
Physical interpretation, 127 
Physicochemical constitution, 58 
Physicochemical interpretation, 
127 

Physicochemical mechanism, 128 
Physicochemical reactions, 69 
Physicochemical state, loi, 123 
Physicomathematical theory of 
visual perception, 179 

Physics, 102 

Physiological conditions, 12 1, 147 
Physiological interpretation, 127 
Physiological properties of the 
neuron, 120 

Physiological role of the nucleus, 
93 

Pi6ron, H., 138, 139 
Pincussen, L., 99 
Pitch, 14S 
Plane, 85, 87, 181 
Plant cells, ^3 
plasmolysis of, 30 

Plastic flow, laws of, 41 
Plate, equatorial, 60 
Pleasantness, measure of, 188 
Polar molecules, 56, 72 
Poles, 47, 60, 61 
Polygonal contour, 187, 188 
Polygonal figure, 194 


Polygons, 189, 190, 192 
Positive ions, 104, 107 
Positive values, no 
Postulates, fundamental, 121 
Potato fragments, respiration of, 27 
Potential, 42, 95 
action, 107 

Preparatory interval, 142 
Preparatory stimulus, 141, 142 
Pressure, 40, 41, 47, 77, 78, 128 
capillary, 52 
components of, 42 
osmotic, 40, 78 
oxygen, 30, 34 

Pressure changes, mechanical, 102 
Presynaptic impulse, 129 
Probability, 168, 169, 172, 176 
for correct statement, 163 

Probability integral, 169 
Problem 
cell-kinetic, 29 
diffusion, 13, 103 

Processes, breakdown, 18 
Produced substance, 14, 45, 54 
Production 

anaerobic lactic-acid, 34 
of carbon dioxide, 54 
of excitatory factor, 164 
of inhibitory factors, 164 
of lactic acid, 54 
of one substance, 94 
rates of, 3, 5, 9, 13, 30, 51, 53, 
94, 97, 106 

Products, intermediate, i8 - 
Propagation, velocity of, 137 
Properties of excitation, loi 
Proportionality 
coefficient of, 20, 31, 104, 183 
constant of, 153 
inverse, 9 

Proprioceptive excitation, 185 
Proprioceptive fibers, i8r, 186 
Proprioceptive impulses, 181, 182, 

183, 19s 
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Proprioceptive stimuli, i 66 
Protein breakdown, 54 

proteins, 78 
complex, 70 

Protoplasm, 4I; 102, 104, 107 
Protoplasmic streaming, 95, 9 ^, 99 
Protozoa, 92 

Psychometric Corporation, x 
Psychophysical discrimination, 1 21 

Quadrant, 87 

Quasi-homogeneous structure, 120 
Qtiasi-solid meshwork, 39 
Quasi-spherical cell, 70 
Quasi-stationary state, 25, 82, 89 

Radii of curvature, 47 
Radius, 67 

Rank-order method, 189 
Rank-order values, 19 1 
Rashevsky, Mrs. Emily, x 
Rashevsky, N., 16, 17, 38, 60, 67, 
68, 81, 92, 107, 1 1 2, 1x8, 1 19, 
X2I, 127, 129, 1335 140, 190 
Rate 

of change, 9, 73 
average, relative, 41 
relative, 52 
constant, 49^ 
of consumption, 3, 15, 7° 
of decomposition, 26 
of decrease, 104 
of deformation, 41 > 43? 
of elongation, 47, 48, 61, 83 
average, 50 

of formation, 16, 33, 72 
of growth, 57, 69, 70, 71, 74, 75, 
77, 79, 80 
of increase, 79, iix 
metabolic, 90 

of metabolism, efiective, 57 

molar, 55 

of oxidation, 19 

of oxygen consumption, 26, 58 

of production, 3, 5, 9, iS, 

SI. SS, 94 . 97 . 106 


proportional 

to average concentration, 19 
to concentration, 15 
to frequency, X31 
of reactions, 72, 94 
constant, 93 

of sugar consumption, 52 
Ratio of length to width, 4, 192 

Reaction, 3X, 33, 37, 93, 9^, 141, 
142, 144, 147, 151, 163, 165, 
x66_ 

reversible, 31, 78 
Reaction rates, 93 
Reaction times, 121, 134 ff- 
Reactions 
cellular, 19 

chemical, 29, 33, 37, 54 , io 4 , ^00 
constant, 29 
rate of, 93 ^ 
coupled chemical, 18 
physicochemical, 69 
rate of, 72, 94 
respiratory, 54, SS, 73 
simultaneous, 18 
wrong, 167 

Real characteristic roots, iii 
Real roots, 15 

Reciprocal innervation, 164 
Recognition 
Gestalt, 195 
of gross features, 19S 
Recti inferiors, 182 
Recti superiors, 182 
Rectus lateralis internus, 182 
Rectus lateralis, externus, right, 
182 

Red blood corpuscles, 88 
Reflexes, conditioned, 145 
Refractiveness, relative, 130 
Refractoriness, 130 
Refractory period, 122, 131 
Refractory state, 122, 130, 131 
Reiner, J. M., 56, 68 
Relations, energy, 4 
Relative intensity, 152 
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Relative oxygen consumption, 22, 

23 

Relative rate of change, 52 
average, 41 
Relaxation, 182, 1S5 
Repulsion, 51 
Reserve of sugar, 25, 26 
Resistance, 14 

of the cell, total diffusion, 14 
diffusion, 26, 28, 29, 30, 32, 33, 
36 

force of, 2 
of medium, 3 
of the membrane, 79 
of the myelin sheath, 118 
per unit area, 98 
viscous, 60, 98 

Resistance forces, viscous, 41 
Respiration, cell, 30 
Respiratory metabolism, 25, 75 
Respiratory metabolites, 78 
Respiratory reactions, 54, 55, 73 
Respirometer, 27 
Response, 185, 194 
correct, 166, 167, 168, 176 
equality, 167, 168 
multiple, 147 
wrong, 166, 167, 168 

Resultant, 54, 56 
Resultant force, 40 
Resynthesis, 19, 30 
Retina, 180 
Reversible reaction, 78 
Rheobase, 105, 114, 118 
Rhythmic excitation, 145 
Riesz, R. R.j 159, 160 
Robertson, D. M., 161 
Root, inverse cubic, 58 

Roots, IS 
negative, 15 
real, 15 

real characteristic, iii 
Rounding up of a cell, 90, 92 
Rushton, W. A. H., 118, 119 


Sac, 37 

Saline solution, 25 
Scale, logarithmic, 58 
Scanning, 183, 184, 185 
Sciences, natural, i 
Segment, 183, 185, 186, 189 
Segments, parallel, 185 
Semi-logarithmic paper, 99 
Sensation of pleasure, 188 
Sense, tactile, 102 
Sense organ, 145 
Sensory center, 180, 181 
Sensory discrimination, 121 
Septagons, 192 

. Shape, 3, s, 14, 28, 43, 48, 63, 90, 

92, 179 

of cell, 4, 23, 82, 92 
of the contour, 181 
of the curve, 142, 194 
detailed, 4 
dumbbell, 51 
of a figure, 1 79 
flattened, 88 
nonspherical, 90 
oblong, 22, 44, 88, 91 
rounded-up, ii 
size and, 2 
spherical, 4 

Shapiro, H., 68 

Sheath, resistance of myelin, 118 
Shoup, 0., 24, 36 

Sides, 6, 7, 43 , 44, 47, 97, 187, 189, 
190 

area of the, 9 
of the cell, 41 

Simultaneous equations, 112 
Simultaneous presentation, 165 
Simultaneous reactions, 18 
Sinusoidal current, 114 
Size, 14, 28, 45, 49, S6, 71, 73, 75, 
7S, 79 
average, 58 
of a cell, 69 
over-all, 7 
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Size — continued 
and shape, 2 

statistical study of cell, 80 
Size distributions, 81 
Slice, multicellular, 28 
Slices of tissues, 28 
Slope, 58, 100 
Snail, 60 
Sol, 41 

Solandt, D., x, 114, ii 5 , iiQ 
Solid cell constituents, 39 
Solid particles, disconnected, 40 
Solid strands, 39 
Solid structures, 40, 95, 96, 98 
Solutes, 37, 38 
Solution, saline, 25 
Solutions, periodical, 16 
Solvent, 37, 38 
Spermatocytes, 60 
Sphere, 3, 15, 51, 87, 88, 89 
perfect, ii, 40 
volume of, ii 

Spherical particle, 39 
Spherical cell, quasi-, 70 
Spherical cells, 4, 22, 45, 46, 84, 103 
Spherical particle, 39 
Spherical shape, 4 
Spheroidal cell, 48 
Spheroidal cell, 49, 50, 54 
Spheroidal shape, 88, 89 
Spindle, 62 
Spindle fibers, 60 
Spontaneous division, i, 4 
Square, 194 
Stability, 45 
Stable equilibrium, 88 
mechanically, 4 
Standard deviations, 81 
Standard weight, 169, 174 
Stationary case 
non-, 83 
quasi-, 82, 89 


Stationary state, 13, 14, 16, 21, 43, 
71, 106, 148 
Stationary water, 40 
Statistical significance of excita- 
tion, 128, 129 

Statistical study of cell sizes, 80 
Stimuli, 125, 130, 13s, 140, 163, 
164, 165, 167, 174 
comparison of, 176 
continuous, 129 
difference of, 163 
enteroceptive, 166 
equal, 163 
external, 145 
exteroceptive, 166 
greater (stronger), 163 
proprioceptive, 166 
Stimulus, 124, 140, 141, 142, 147, 
148, 153, 136, 170 
constant, 148, 177 
continuous, 125 
intensity of, 135, 138, 168 
external, 179 

intensity of, 122, 129, 138, 155, 
178 

natural, 123 

peripheral, 121, 132, 137, 151, 
178, 179 

preparatory, 141, 142 
warning, 141, 142, 143 

Stimulus pattern, 151 
Straight line, 58, 183, 186 
segment of, 179, 181, 185, 187 
Streaming, 96, 98 
protoplasmic, 93 
velocity of, 96, 97, 98, 99 
Structural variations, local, 60 
Structure, 37, 39, 40 
anatomical, 128 
colloidal, 37 
complex, 28, 120, 121 
immobile, rigid, 95 
inhomogeneities of, 93 
internal, 5 
molecular, 37 
neuronic, 139, 145 
rigid, 90, 95 
solid, 95, 96, 98 
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Subliminal intensityj 129 
Subscripts, 32 
Substance, 5, 13, 69, 70, 72 
chemical, 102 

concentration of a, 5, 14, 15 
diffusing, 37, 38 
flow of, 2 

nonmetabolized, 78 
produced, 14, 45, 54 
production of one, 94 

Substances, factors as, 8, 123 

Sugar, 18, 19, 20, 25, 30, 35, 37 
consumption, rate of, 32 
external concentration of, 32 
reserve of, 25, 26 

Sugars, concentration of different, 
34 

Superliminal intensity, 129 
Surface, 41, 42, 47, 52, 98 
concentration at the, 43 
Surface integral, 42 
Surface membrane, 30 
Surface tension, 45, 46, 47, 48, 50, 
58, 60, 63, 83 
lateral, 52 
Suspension, 27 
Symmetric functions, 174 
Symmetric scheme, 165 
Symmetry, 40, 187, 192, 193 
axial, 42, 103 

Synapse, 123, 124, 126, 130, 131, 

134, 13s, 138, 140, X41, 3:49, 

ISO, 151, 152, IS3, 163, 165, 
177, 180, 193 

Synaptic delay, 129, 130, 135, 136, 
138 

Synaptic impulse, pre-, 129 
Synaptic transmission, 144 

Tactile sense, 102 
Tang, P. S., 24, 25, 27, 28, 36 
Taste nerve fibers, 102 
Tastes, differences in, 192 


Temperature, 28, 98, 99 
absolute, 38 
gradient of, 126 
variation with, 28 

Terminology, 167 
Theory 

of discrimination, 152, 194 
of excitation, 105 
of heat flow, 127 
kinetic, 127 
two-factor, no 

of visual perception, physico- 
mathematical, 179 

Thermal conductivity, 126 
Thermal motion, 126 
Thinking 

of horizontal objects, 195 
of a vertical object, 195 

Threshold, loi, 103, 105, 116, 122, 
129, 130, 131, 135, 1447 1457 
147, 149, 150, 153, 165, 166, 
169, 171, 184 

Threshold concentrations, 102 

Time, it, 43, 60, 64, 69, 71, 74 , 
75, 105, 106, 131, 141, 14s, 
169, 176, 194 
conduction, 141 
excitation, 116, 129 
reaction, 121, 139, 141, 142 
of transmission, 138, 139 
variation with, 135 

Time relations, no, 125 

Tissue, 23, 70, 75 
heart, 58 

Tissues 

contractility of, 92 
normal, 80 
slices of, 28 
tumor, 75 

Tongue, 148 
Tonic contraction, 185 
Total diffusion resistance of the 
cell, 14 

Transcendental equation, 16 
Transformations, 156 
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Transmission 
direction of, 134, 135 
synaptic, 144 
time of, 138, 139 
Transport, 2, 20, 21 
convection, 97 
diffusion, 97 
heat, 126 

Tschermak, A., 196 
Tumor tissues, 75 
Tumors, 56, 69 
benign, 80 
rapidly growing, 80 
Two-factor theory, no, in, 112, 
118, 129, 131, 132 

Two factors, 123, 135 
Two-parametric curve, 143, 176 

Umrath, K., 119 

Unicellular organism, variation of 
shape with metabolism, 91 
Units, 152, 153, 191 
absolute, 7 

Unstable equilibrium, 88 
configuration of, 89 
Upper limit, 96, 98 
Urban, F. M., 174, 175, 176 

Variables, 104, 156 
Variations, S, 6, 58 
individual, 12 
local structural, 60 
with temperature, 28 
with time, 135 
Velocity 

of conduction, 118 

of contraction, 182, 183, 186 

of flow, 20 

gradient of, 98 

of propagation, 137 

of streaming, 96, 97, 98, 99 


Viscosity, 41, 95, 98 
Viscous drag, 38 
Viscous incompressible body, 53 
Viscous resistance, 60, 98 
Viscous resistance forces, 41 
Viscous strands, 39 
Vision, field of, 195 
Visual data, 139, 142, 176 
Visual perception, physicomathe- 
matical theoiy of, 1 79 

Volume, II, 38, 39, 41, 42, 43 j 49 , 
54, 58, 59, 64, 73 y 9 < 5 , 97 y 128 
of body, 41 
of the cell, 9 
element of, 40, 50 
of the sphere, 4 
Volume force, 42 
Volume integral, 42 

Walls of the cell, 95 
Warburg, 0 ., 56, 57, 68, 69, 81 
Warning stimulus, 141, 142 
Water, i8, 20, 31, 37, 38, 39, 40, 
„ 79 . 98 
flowing, 5S 

Weber ratio, 154, 156 ff. 

Weight 

discrimination of, 169, 173 
molecular, 31, 38, 54, 55 
standard, 169, 174 
Weinberg, Alvin, ix, x, 16, 17, 1 12, 
118, 119 

Width, ratio of length to, 4 
Williamson, Robert R., ix, x, 68 
Woodrow, H., 143 
Wrong response, 166, 168 

Young, Gale, ix, x, 12, 13, 16, 38, 
42, 46, SI, 52, 63, 68, 96, 100, 
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